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ABSTRACT 


An  experiment  has  been  conducted  to  determine  the  sensitivity  of  an 
HF  FSK  communications  system  to  time-  and  frequency- sel ecti ve  lading  phe¬ 
nomena.  Binary  error  rate  wa;  measured  on  AN/FGC-29  teletype  terminal 
equipment  over  an  HF  path  from  Fort  Monmouth,  New  Jersey  to  Palo  Alto, 
California  in  the  fall  of  1964.  Simultaneous  measurements  were  made  of 
channel  si gn al - to- noi se  ratios,  channel  time-delay  profiles,  channel 
Dopp 1 er- sh i f t  profiles,  and  signal  crrrelation  between  spaced  antennas. 
Instrumentation  included  a  phase-stable  t r an smi t ter- recei ver ,  an  oblique- 
incidence  ionospheric  sounder,  and  a  special  correlation  meter.  The  mea¬ 
sured  error  rates  were  compared  to  those  predicted  by  a  mathematical  model. 
The  model  for  the  channel  considers  randomly  time-varying  and  frequency- 
selective  effects  appropriate  for  the  HF  propagation  mechanism.  The  system 
model  approximates  the  AN/FGC-29  system.  Theoretical  and  measured  error- 
rate  curves  are  presented  and  discussed.  Measured  system  performance 
showed  a  well-defined  sensitivity  to  channel  time-  and  frequency-selective 
effects.  The  asymptotic  error  rates  measured  at  high  s i gna  1  - to-noise 
ratios  are  in  reasonable  agreement  with  those  predicted  from  the  model. 
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I  INTRODUCTION 


A.  HISTORICAL  BACKGROUND 

Since  1956,  the  Communication  Laboratory  at  Stanford  Research 
Institute  has  been  ective  in  various  fields  related  to  the  analysis  of 
strategic  military  communications  networks.  This  activity  has  resulted 
in  a  series  of  computer  programs  that  predict  the  performance  and  vulner¬ 
ability  of  such  networks  under  a  wide  variety  of  propagation  and  electro¬ 
magnetic  environments.  1,2*3’4  Recently,  this  work  has  been  expanded  to 
include  consideration  of  time-  and  frequency-dispersive  effects  on  the 
performance  of  selec  ed  types  of  systems.  Since  the  backbone  of  strategic 
military  communications  networks  is  HF  radio — even  though  there  has  been 
recent  inclusion  of  t roposcatter ,  cable,  and  satellite  systems — the 
emphasis  has  been  on  HF  systems  and  propagation.  The  computational  models 
used  in  the  propagation  predictions  are  based  on  predictable  geophysical 
factors.  As  a  result,  predictions  of  signal  -  to-noi  se  ratio  (S/.N)  on  a 
radio  path  can  be  computed  months  and  years  in  advance.  Even  though  these 
predictions  are  not  exact,  they  are  valuable  in  day-to-day  station  operation 
and  in  the  planning  of  future  networks. 

To  transform  the  results  of  the  propagation  predictions  into  estimates 
of  communications  link  performance,  it  was  necessary  to  establish  a  com¬ 
putational  •’’odel  for  each  communications  system  used  in  the  particular 
terminals  of  the  network.  It  was  also  necessary  to  establish  a  quantitative 
definition  of  the  performance  of  these  systems  and  criteria  for  the  utility 
of  the  link  based  on  this  per  foi  .nance  rating.  The  measure  of  performance 
adopted  was  binary  error  rate.  This  parameter  has  the  following  advantages: 

(1)  It  can  often  be  evaluated  analytically  from  basic  system 
models  by  making  certain  assumptions  on  the  signal  and 
no i se  statistics. 

(2)  It  can  usually  be  translated  further  to  other  user-defined 
performance  criteria,  -uch  as  character  error  rate, 
reliability,  or  message  intelligibility. 

(.11  i*  can  be  directly  estimated  from  nonsub  i  <•«  1 1  ve 
m**asu rement  s 
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With  the  increasing  use  of  data  and  digi tal - voi ce  transmission,  binary 
error  rate  will  undoubtedly  become  the  standard  measure  of  communications 
circuit  performance. 

Recently,  the  Communication  Laboratory  has  been  active  in  deriving 
more  detailed  error-rate  models  of  specific  systems  designed  for  HF.  These 
models  are  unique  in  that  they  consider  dispersive  channel  effects  but  re¬ 
quire  the  specification  of  additional  channel  parameters  not  predicted  by 
the  present  propagation  programs.  It  is  hoped  that  these  propagation 
programs  can  be  updated  by  including  the  additional  parameters  so  that 
eventually  a  computation  capability  for  evaluating  entire  networks  by  using 
these  new  system  and  channel  models  can  be  obtained. 

To  gain  confidence  in  error-rate  computational  models,  a  direct  com¬ 
parison  with  actual  measured  results  was  needed.  Stanford  Research 
Institute  was  asked  to  conduct  a  controlled  experiment  to  determine  the 
accuracy  of  the  error- rate  computational  model  for  a  particular  system 
when  ali  of  the  required  channel  characteristics  are  known.  This  experi¬ 
ment  was  performed  in  the  fall  of  1961;  results  of  this  experiment  are 
presented  in  this  report. 


B.  REVIEW  OF  ERROR-RATE  MODELS 


Error-rate  computational  models  are  based  on  statistical  principles. 

As  such,  they  do  not  attempt  to  predict  the  exact  number  of  errors  that 
will  occur  but  only  to  answer  the  question:  What  is  the  probability  that 
a  receiver  will  decide  that  State  A  was  transmitted  luring  a  particular 
signaling  interval  wher  State  B  was  actually  transmitted?  To  answer  such 
a  question,  probabilistic  concepts  must  be  introduced  and  with  them  certain 
random  variables  that  can  he  described  only  in  a  statistical  sense. 

The  first  attempts  to  derive  an  error- rate  model  for  communications 
channels  were  based  on  tier  following  assumot  i  o'is : 


(1)  An  ideal  mat  « -heel  -ft  1 1  ci  receiver 
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interval.  Coherent  matched- fi 1  ter  detection  requires  that  the  receiver 
have  exact  knowledge  of  all  possible  phase  states  of  the  transmitted  wave¬ 
forms  when  they  arrive  at  the  receiver.  Incoherent  matched- f i 1  ter  detec¬ 
tion  does  not  utilize  a  priori  phase  information  of  the  received  signal. 

The  assumption  of  an  ideal  matched- fi 1  ter  receiver  implies  exact  time  syn¬ 
chronization  with  the  transmitted  data  stream  and  ideal  integration  and 
sampling  of  the  received  waveform.  This  is  followed  by  complete  energy 
dumping  or  quenching  of  all  integrating  filters. 

The  next  evolutionary  stage  in  the  development  of  error-rate  models 
was  to  allow  the  received  signal  envelope  to  vary  slowly  with  respect  to 
a  signaling  element  ( i.e to  fade)  and  to  assume  a  distribution  of  this 
variation  in  calculating  the  effect  on  the  probability  of  error.  It  was 
also  assumed  in  these  models  that  fading  over  the  transmission  bandwidth 
of  the  signal  was  correlated.  Such  a  channel  has  been  termed  flat-flat 
by  Bello  and  Nelin,5  because  the  time  autocorrelation  function  of  the 
channel  is  constant  over  a  time  difference  on  the  order  of  a  signaling- 
element  time  interval,  and  the  channel  frequency  -  cor rel at i on  function  is 
constant  over  a  frequency  difference  comparable  to  the  bandwidth  of  the 
transmitted  signal.  The  frequency  correlation  is  directly  related  to  the 
time-delayed  “ mul tipath  ”  propagation  modes  which  are  observed  at  HF  and 
can  be  predicteu  by  propagation  theory.  The  time  autocorrelation  is 
directly  related  to  the  difference  in  Doppler  shift  associated  with  these 
modes.  Diversity  improvements  were  modeled  by  assuming  the  availability  at 
the  receiver  of  independent  fading  signals  and  specifying  various  combining 
rules  for  the  receiver.  The  flat-flat  models  show  a  decreasing  error  rate 
with  increasing  S/N. 

More  recent  channel  models  include  the  effects  of  Doppler  and  delay 
spreading  of  the  transmitted  signals.  The  use  of  these  models,  however, 
requires  additional  statistical  parameters  to  describe  the  spreading  of 
the  signal  by  the  channel  in  time  and  frequency.  These  models  show 
asymptotically  constant  error  rates  with  increasing  SN.  The  asymptotic 
values  depend  critically  on  interaction  of  certain  system  parameters  and 
on  Use  statistical  channel  spread  measures. 

\.  HF.  the  simpler  flat -flat  models  have  generally  been  used  for 
present  binary  FSR  *•  stems.  These  are  adequate  however,  only  when  gro«s 
estimates  of  pr  r  formal!' e  at  moderate  fi  ror  lates  are  needed.  For  normal 
100-wpm  teletype  channels  predominately  handling  messages  with  F.ngltsh- 
langiiage  redundancy,  i;  binary  error  rale  less  than  10  5  has  usually  been 
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considered  acceptable.  With  the  present  increase  in  less  redundant  data 
traffic  and  the  corresponding  increase  in  reliability  requirements,  the 
simpler  models  are  becoming  inadequate.  As  newer  modulation  systems  and 
higher  data- transmi ssion  rates  are  introduced,  time  and  frequency  spreading 
of  the  channel,  even  at  HF,  cannot  be  ignored  in  a  performance  model.  But 
estimating  or  predicting  these  channel  scattering  functions  from  geophysical 
parameters  presents  a  problem.  Although  the  present  programs  can  predict 
individual  signal  strengths  for  each  propagating  mode,  together  with 
propagat ion- time  differences,  they  cannot  adequately  predict  Doppler  effects. 

C.  OUTLINE  OF  THE  VALIDATION  EXPERIMENT 

The  basic  plan  for  the  validation  experiment  was  to  take  a  typical  HF 
FSK  data  modem,  a  single  channel  from  a  standard  AN/FGC-29  system,  and 
adapt  an  appropriate  computational  model  for  it.  from  this  model,  a  series 
of  binary  error-rate  characteristics  were  calculated  with  various  values 
of  channel  parameters.  The  next  step  was  to  run  an  experiment  on  an  actual 
HF  link  using  a  system  of  this  type,  to  measure  continuously  all  of  the 
pertinent  channel  parameters,  together  with  obtaining  a  measured  estimate 
of  the  binary  error  rate. 

Chapter  II  discusses  the  computational  system  and  channel  models  used 
and  the  system  and  channel  descriptors  that  must  be  specified.  Chapter  III 
describes  the  actual  experiment  performed,  including  a  description  of  the 
instrumentation  and  of  the  data- recording  procedure.  Chapter  IV  outlines 
the  data- pi ocessing  operations — the  criteria  and  procedure  used  to  obtain 
confident  measures  of  the  channel  parameters  and  the  binary  error  rate. 
Chapter  V  presents  the  results  in  the  form  of  computed  and  measured  error- 
rate  curves,  together  with  a  discussion  of  the  results. 

Chapter  VI  snows  the  results  of  an  allied  experiment  that  demonstrates 
quite  dramatically  the  correlation  between  f requeue y - se l ec t i ve  effects  and 
detected  errors  on  an  HF  communications  system 

Chapter  VII  presents  some  results  on  the  measurement  of  the  correla¬ 
tion  of  signal  lading  as  received  on  two  space  - sepai at ed  antennas.  Com¬ 
parison  of  these  results  with  simultaneously  measured  *rrors  on  a  dual* 
diversity  system  and  with  measured  time -da  I  ay  'lur>u  trust  ns  on  tin* 
channel  are  also  shown. 
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II  DESCRIPTION  OF  THE  COMPUTATIONAL  MODEL 


A.  I  INTRODUCTION 

The  computational  model  used  to  compute  the  theoretical  error-rate 
curves  was  developed  in  a  prior  study.6  In  this  study,  the  effects  of 
time-  and  f requency - se 1 ec t i ve  fading  on  an  idealized  KSK  system  were 
analyzed  for  a  simple  HF  channel  model.  For  the  purposes  of  the  present 
effort,  the  computational  model  and  associated  computer  programs  were 
modified  to  include  the  important  effects  >n  a  non-zero  transition  time 
when  the  FSK  transmitter  switches  from  one  frequency  to  another.  The 
model  is  more  realistic  than  the  previous  one  which  assumed  an  instanta¬ 
neous  frequency  change.  It  was  found  that  the  inclusion  of  a  small  transi¬ 
tion  time  (approximately  2  percent  of  the  si gna 1 i ng- e 1 emen t  duration) 
caused  the  predicted  performance  of  the  FSK  system  to  display  a  moderate 
deterioration  (as  opposed  to  a  phase-continuous  instantaneous  frequency 
change)  at  high  S/Ns  for  typical  ionospheric  time-delay  and  Doppler 
condi lions. 

B.  SYSTEM  AND  CHANNEL  MODELS 

The  system  model  assumes  a  s i ngl e - channe 1  matched- f i 1 te r  receiver 
with  the  filters  matched  to  ideal  FSK  transmissions.  A  quadratic  de¬ 
tector  obtains  the  squared  envelope  of  ihe  output  of  each  matched  filter; 
the  receiver  then  decides  which  signal  was  transmitted,  on  the  basis  of 
the  larger  squared  envelope. 

With  the  exception  of  a  linear  variation  of  frequency  in  the  transi¬ 
tion  intervals,  the  FSK  transmission  is  assumed  to  be  ideal.  If  i  is 

-  • 

the  m.irk  frequency  and  f  is  the  space  frequency,  then  Fig.  I  displays 
the  linear  frequency  transition  when  the  sequence  SPACE -MARK' SPACE  IS 
transmitted.  In  this  figure,  8  is  the  s  i  gna  l  t  ng  -  e  I  >'nn n  t  duration,  and 
2 L  is  the  transition  lime.  Observe  that  in  the  transition  i nt er v a l s  t hr 
transmission  i a  chirp  signal  that  spans  all  frequent  les  between  the* 
two  signaling  frequencies.  The'  receiver  filters  (which  are  matched  to 
ideal  \  .‘>K  transmissions)  are  not  matched  t  •*  l  !.«■  above  transmissions 
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FIG.  1  LINEAR  FREQUENCY  VARIATION  AT  TRANSITION  INTERVALS 

As  the  transition  time,  2 L,  is  increased,  the  performance  of  the  FSK 
system,  operating  through  a  time-varying  and  f requency- sel ec t i ve  radio 
channel,  deteriorates.  For  typical  ionospheric  Doppler  and  time-delay 
conditions  and  small  transition  times,  the  deterioration  appears  as  an 
increase  in  the  binary  error  probability  for  high  S'Ns.  The  randomly 
selective  radio  channel  transforms  the  chirp  transmissions  at  the  tran¬ 
sit  i  on  t imes to  noise-like  signals  which  impair  the  detectability  of  the 
FSK  transmissions. 

The  channel  model  has  been  adequately  described  in  previous  Stanford 
Research  Institute  reports.6'7  This  model  assume''  that  the  randomly  time- 
varying  linear  channel  is  composed  of  a  discrete  number  of  independent 
propagation  paths.  Hath  path  is  described  by  a  time  delay,  a  Doppler 
shift,  and  a  random  complex  gain.  The  model  tan  be  considered  as  an 
idealized  version  of  the  HF  channel  when  the  time-delav  spread  and  Doppler 
spread  of  the  individual  paths  are  small  compared  to  the  time-delay  dif¬ 
ferences  and  Doppler- shi ft  differences  between  the  paths  Alternatively, 
as  the  number  cf  paths  bet  omes  large,  this  channel  approaches  one  with  a 
continuum  of  time  tie  lay*  ami  Doppler  shifts.  The  channel  model  accounts 
for  the  important  signal  distort  mg  effett*  caused  bv  time-  and  frequency- 
selective  failing,  su'cti  as  intersymbol  mt  cr  ferrh*  c  and  last  lading  o  f  the 
signal  amplitude  within  a  s ign t I t ng - e i eneni  duration 
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C.  COMPARISON  OF  THE  COMPUTATIONAL  SYSTEM  MODEL 
AND  THE  EXPERIMENTAL  SYSTEM 

The  AN/FGC-29  FSK  receiver  used  in  the  experiment  and  the  incoherent 
matched- fi 1  ter  FSK  receiver  of  the  computational  model  display  notable 
differences  in  detection  procedures.  The  AN/FGC-29  employs  discriminator 
detection  without  hard  limiting  and,  in  contrast  to  the  incoherent  matched- 
filter  receiver,  does  not  quench  the  energy  stored  from  previous  signaling 
elements  before  processing  any  given  received  signaling  element.  In  the 
absence  of  hard  limiting,  the  discriminator  can  be  modeled  with  two  linear 
narrow-band  filters  with  the  appropriate  center  frequencies  and  bandwidths 
for  the  two  transmitted  FSK  signals.  Hence,  one  of  the  major  differences 
between  the  computational  model  and  the  experimental  system  lies  in  the 
energy-quenching  property  of  the  computational  model.  In  the  experimental 
system,  the  signal  energy  due  to  previous  signaling  elements  decays  at  a 
rate  determined  by  the  filter  time  constant.  In  the  computational  model, 
the  energy  in  the  filters  is  quenched  instantaneously  e.ery  I F  seconds,  and 
detection  is  based  only  on  the  signal  received  in  a  time  interval  of  length 
H  seconds.  If  the  filters  of  the  experimental  system  do  not  decay  stored 
energy  rapidly  enough,  it  is  possible  that  f requency - se 1 ec t i ve  fading  could 
cause  the  stored  energy  in  one  filter  to  differ  markedly  from  the  stored 
energy  in  the  other  filter.  This  anomalous  situation  would  cause  an  ini¬ 
tial  bias  error  which  could  adversely  affect  the  system’s  detection 
capability.  The  computational  model  dees  not  account  for  an  unequal  biasing 
of  the  filters  before  the  processing  of  a  signaling  element;  however,  it 
does  account  for  the  effects  of  time-  and  f requency - se 1 ec t i ve  fading  oc¬ 
curring  during  a  si gnal i ng- e 1 emen t  duration.  In  filtering  the  received 
pulses,  the  experimental  system  smooths  and  shapes  the  time  and  frequency 
characteristics  of  the  received  signal  instead  of  essentially  narrow-band 
integrating  over  a  s i gna 1 t ng- e 1 emen t  duration  as  does  the  computational 
mode  1 . 

1).  THE  ERROR  RATE  PROGRAM 

The  error- rate  program  utilises  the  simple  HF  channel  model  and  the 
idealised  FSK  system  model  (including  the  t runs i t i on - 1 t me  modification) 
to  compute  a  binary  error  probability.  Input  parameters  to  the  program 
describe  the  particular  channel  scattering  function'  desired  in  the  error- 
rate  analysis.  These  parameters  include  the  number  of  propagation  paths 
and  a  relative  strength,  time  delay,  the  Doppler  shift  for  each  path. 
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In  the  computation  oi  the  binary  error  probability,  it  is  assumed  that  the 
channel  transfer  function  is  a  homogeneous  Gaussian  random  field7  and  that 
the  output  of  the  channel  is  further  perturbed  by  an  additive  white  Gaussian 
noise.  The  details  of  this  computation  can  be  found  in  Ref.  6,  where  the 
cross- ambi  gui  ty  function  [Eq.  (33),  Ref.  t>]  has  been  modified  to  include  the 
transition-time  effect.  It  is  assumed  that  the  time  delays  (measured  from 
the  average  time  delay  determined  by  the  time-delay  profile  of  the  channel) 
of  the  various  propagation  paths  are  all  less  than  a  signaling-element 
duration.  Hence  intersymbol  interference  is  limited  to  adjacent  signaling 
elements.  The  program  computes  the  probability  that  the  MARK  in  the  trans¬ 
mitted  sequence  SPACE -MARK -SPACE  is  detected  improperly  as  a  SPACE.  This 
sequence  corresponds  to  the  actual  alternating  mark  and  space  sequence  that 
was  employed  in  the  experiment,  which  was  chosen  to  maximize  the  effects  of 
frequency- selecti ve  fading  (intersymbol  interference  and  unequal  fading  of 
the  mark  and  space  frequencies). 

The  output  of  the  program  is  a  computer-plotted  curve  of  binary  error 
probability  versus  S/N  for  the  desired  channel  scattering  function.  Each 
plot  displays  two  curves,  one  for  no  diversity  (labeled  FSKl )  and  one  for 
independent  dual  diversity  (labeled  FSK2 ) .  Each  plot  displays  a  channel¬ 
scattering-function  diagram  in  the  lower  left-hand  corner.  In  this  diagram 
each  path  is  represented  by  an  X  with  size  proportional  to  the  strength  of 
the  path  and  with  position  determined  by  the  time  delay  and  Doppler  shift 
of  the  path.  Time  delays  relative  to  the  average  time  delav  of  the  channel 
are  read  along  the  logarithmic  horizontal  scale  and  are  normalized  to  the 
signaling-element  duration  (T  -  the  time  delay  of  the  path  divided  by  the 
s ignal i ng- e 1 ement  duration).  Doppler  shifts  relative  to  the  average 
Dopph  (  shift  of  the  channel  are  read  along  the  logarithmic  vertical  scale 
and  are  also  normalized  to  the  signaling-element  duration  ( /)  -  Doppler  shift 
of  the  path  multiplied  by  the  s i goal i ng - el ement  duration).  For  example. 

Fig.  2(a)  displays  the  error-rate  curves  for  an  FSk  system  operating  through 
a  propagation  medium  composed  of  two  equal  - st rengt h  paths  spaced  at  a  time- 
delny  difference  equal  to  0.2  of  the  s i gnn 1 i ng - e 1 ement  duration  and  with  no 
Doppler  spreading  in  frequency. 

F.  CM AN NEC  -  SCATTER  I  Mi  -  Kl'NCT  l  ON  APPROX  I  MAT  I  ON 

Validation  of  the  theoretical  error- rate  analysis  required  simulta¬ 
neous  measurement s  of  the  channel  scattering  func  tion.  >  N.  and  binary  error 
probability.  The  scope  of  the  experiment  (the  a%  a  t  1  .»!>  t  >*  time  and  rate  oi 
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effort)  made  it  neceisary  to  limit  the  measurement  of  the  channel  scat¬ 
tering  function  to  the  time-delay  and  Doppl er-sbi ft  profiles.7  It  was  also 
necessary  to  concentrate  the  major  portion  of  the  effort  on  obtaining  data 
for  the  no-diversity  case,  since  the  dual -di versi ty  case  would  add  another 
dimension  to  the  classification  of  data:  spacial  correlation.  The  amount 
of  data  required  to  obtain  statistically  reliable  results  for  dual - di versi ty 
operation  at  all  possible  states  of  S/N,  time  delay,  Doppler  shift,  and 
spacial  correlation  far  exceeded  our  si ngl e- channel  capability  for  collecting 
dft^a  and  available  rate  of  effort  for  data  processing.  However,  a  limited 
amount  of  data  was  obtained  for  dual-diversity  operation;  some  interesting 
results  concerning  spacial  correlation  are  discussed  in  Chapter  VII. 

Although  time-delay  and  Doppler-shi  ft  profiles  were  measured  on  v,he 
channel  throughout  the  experiment,  it  was  not  feasible  or  necessary  to  use 
all  of  the  information  that  these  functions  contain  for  correlating  the  ob¬ 
served  error  rate  with  the  state  of  the  dispersive  channel.  Various  com¬ 
puter  programs  were  run,  investigating  the  sensitivity  of  the  error  rate 
displayed  by  the  theoretical  model  to  the  shape  of  the  scattering  function. 

It  was  found  that  if  the  scattering  function  is  approximately  symmetrical 
and  of  moderate  width  in  the  time-delay  and  Doppl er- shi ft  directions,  then 
the  error  rate  is  strongly  dependent  on  the  time-delay  and  Doppl er- shi ft 
second  moments  and  relatively  insensitive  tc  the  scattering- function  shape. 
Figure  2  illustrates  this  phenomenon  in  the  time-delay  direction. 

Figure  2(b)  and  (c)  corresponds  to  error  rates  computed  for  a  time-delay 
profile  consisting  of  twenty  propagation  paths  wi*h  energy  distributed 
uniformly  and  triangularly,  respectively.  In  Fig.  2(a)  the  error  rate  is 
computed  r >r  two  equal - st rength  paths.  In  all  three  of  the  above  scattering 
functions,  the  standard  deviation  of  the  time  delay  profile  is  equal  to  0.1 
of  the  signal i ng-element  duration,  or  ul!  three  shapes  possess  an  rmu  time- 
delay  spread  equal  to  0.2  of  the  s igna l i ng-e I ement  duration.  For  most 
practical  applications,  the  error-rate  curves  of  Ftg  2(a)  are  reasonable 
approximations  to  the  curves  of  Ftg.  2(b)  and  (<  ),  with  ,i  tendency  to 
underestimate*  the  error  rate,  lienee,  symmetrical  -nattering  functions 
with  moderate  lime-delav  and  floppier  spreads  i  iin  l>  *  approximated  with  two 
equa  1  -  si  rength  paths  spared  tn  t  ime  delay  at  the  t  tme-de|;»y  spread  of  the 
s< altering  function  and  spat  ed  tn  Doppler  «t  the  fn.pplet  -plead  of  the 
scattering  function.  This  approximation  introduce-  a  major  s imp  I  I f it  at  ton 
in  the  analysis  d  the  da l a ;  when  valid,  it  alio*-  one  to  um-  only  the 
second  moments  associated  with  the  time-delay  and  Dopplef  profile,. 
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For  scattering  functions  possessing  large  spreads  with  strongly 
asymmetrical  distributions,  the  tendency  of  this  approximation  to  under¬ 
estimate  the  error  rate  becomes  excessive.  For  example,  Fig.  3(a)  dis¬ 
plays  the  error-rate  curves  for  a  channel  with  two  unejual  paths  spaced 
at  0.5  of  the  signal ing- element  duration,  »i th  one  path  five  times  the 
strength  of  the  other  path.  The  rms  time-delay  spread  of  this  channel  is 
0.35  of  the  signaling-element  duration.  Application  of  the  approximation 
of  two  equal  -  strength  paths  spaced  at  0.35  of  the  signal i ng- element  dura¬ 
tion  yields  the  curves  of  Fig.  3(b)  which  excessively  underestimate  the 
error  rate.  In  analyzing  phase  systems,  Gaarder  of  Stanford  Research 
Institute  found  that  asymmetrical  scattering  functions  of  the  type  illus¬ 
trated  in  Fig.  3(a)  can  be  effectively  approximated  with  two  equal -strength 
paths  spaced  at  the  identical  spacing  of  the  original  unequal  -  strength 
paths  and  that  this  approximation  remains  valid  for  a  wide  range  of  relative 
path  strengths. 8  Application  of  this  approximation  to  Fig.  3(a)  yields  the 
curves  of  Fig.  3(c)  which  were  computed  for  two  equal  -  strength  paths  spaced 
at  0.5  of  the  signaling-element  duration.  For  this  case,  the  curves  of 
Fig.  3(a)  and  (c)  are  essentially  identical. 

F.  TRANSITION-TIME  SENSITIVITY  AND  IRREDUCIBLE 
ERROR  PROBABILITY 

The  sensitivity  of  system  performance  as  a  function  of  transition 
time  is  illustrated  by  Fig.  4.  In  thi3  figure,  error-rate  curves  are 
displayed  for  a  fixed  scattering  function  that  describes  a  channel  com¬ 
posed  of  two  equal -strength  paths  spaced  in  time  delay  at  0.1  of  the 
signaling-element  duration  and  with  zero  Doppler  shift.  Deterioration  of 
system  performance  with  increasing  transition  time  is  clearly  indicated; 
the  transition  times  are  0.005,  0.01.  and  0.02  of  the  signaling-element 
duration  in  fig.  4(a).  (b).  and  (c).  respectively'. 

A  significant  property  of  the  error- rate  curves  for  an  incoherent 
matched*  f » 1  ter  system,  opt  rating  through  a  dispersive  channel,  is  the 
existence  of  an  irreducible  error  probability  at  infinite  ,S  \s.  Hence 
the  irreducible  error  probability  serves  as  a  meaningful  measure  of  the 
limitations  imposed  on  system  performance  I,  the  dispersive  properties 
of  the  propagation  medium  trreducible  error  probability  as  a  function 
of  the  time-delay  spread  of  two  equa I - st rengt n  paths  with  xero  IVippler 
is  plotted  in  fig  5.  Similarly,  irreducible  error  probability  os  a 
function  of  the  Doppler  spread  of  two  equal • st rength  paths  with  a 
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APPROXIMATION  OF  A  STRONGLY  ASYMMETRICAL  TIME-DELAY  PROFILE 


oot  OK)  !00  lOOO 

TIME-DELAY  SPREAD  (  NORMALIZED  TO  THE  SI6NAUN&-EI.EMENT  OURATION) 


FIG.  5  IRREDUCIBLE  ERROR  PROBABILITY  AS  A  FUNCTION  OF  TIME-DELAY 
SPREAD,  2<’a  -  0,  21  -  0.C2 
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time-delay  spread  equal  to  the  transition  time  is  plotted  in  Fig.  6. 

Both  the  time-delay  and  Dopp  1  e  r  -  sh  i  f  t.  i  r  reduc  i  bl  e- er  rn  r- p  robabi  •  1  ty 
curves  were  plotted  for  a  transition  time  equal  to  0.02  of  the  signaling- 
element  duration. * 

The  irreducible-error-probability  curves  exhibit  an  interesting 
transi ti on- time  phenomenon.  For  small  time-delay  or  Dopp 1 e r- spreads  and h  i  gh 
S/Ns,  the  system  performance  is  essentially  limited  by  transition  time; 
this  limitation  is  illustrated  by  the  flat  portion  of  the  irreducible- 
error-probability  curves  at  low  spread  values.  At  lower  spread  values, 
the  curves  maintain  a  constant  error  probability;  at  higher  spread  values, 
the  irreducible  error  probability  increases  as  the  spread  raised  to  some 
positive  power. 


*  fli*  ruri'i  4o  not  tltap'av  irredaetbl*  arr<  r  probability  fur  tiaw-4tlay  *pr»a4a  l*aa  than  (Ha  ir<a>iti°a 
laaaa.  A*  tba  UM-falay  -praad  aa  fanaixf  bol««  tKa  liaaaiiitia  t:aa,  lh»  < r rtdut iblt  arror  prob¬ 
ability  bagiaa  to  4*cr«aaa  again. 
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FIG.  6  IRREDUCIBLE  ERROR  PROBABILITY  AS  A  FUNCTION 
OF  DOPPLER  SPREAD,  2<ir  -  0.02,  2L  -  0.02 
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Ill  DESCRIPTION  OF  THE  EXPERIMENT 


A,  INTRODUCTION 

The  purpose  of  the  experimental  program  was  to  compare  measured 
error-rate  S/N  characteristics  with  predicted  characteristics  derived  from 
a  theoretical  error-rate  model  for  an  HF  FSK  communications  system.  The 
design  of  the  experin,  nt  was  so  formulated  that  the  experimental  program 
could  logically  be  divided  into  several  distinct  subexperiments.  All  the 
subexperiments  were  performed  simultaneously,  and  the  results  were  used 
to  compare  measured  error  rates  with  predicted  error  rates  for  the  mea¬ 
sured  c  hannel  condi  t  ions  .  The  subexperiments  will  be  catagorized  as 
fol lows : 

(1)  Error-rate  and  S/N  measurement 

(2)  Doppler  profile  measurement 

(3)  Time-delay  profile  measurement. 

t 

Each  measurement  with  its  associated  equipment  is  described  in  detail  in 
this  chapter. 

The  measurements  were  performed  on  a  single  binary  data  channel  using 
two-stage  modulation.  The  input  binary  data  frequency-modulated  a  sub- 
carrier.  The  resulting  FSK  tones  single- si deband-modul ated  the  HF  carrier, 
with  a  substantial  portion  of  the  carrier  being  retained.  The  data  rate 
was  100  baud,  which  is  slightly  higher  than  that  of  a  single  100-wpm  tele¬ 
type  channel.  The  operating  frequencies  were  It. 360  and  7.366  Me,  used 
for  approximately  the  same  amount  of  time.  These  frequencies  were  chosen 
to  maximize  the  amount  of  data  collected  in  the  scheduled  experiment  time 
and  to  obtain  data  under  a  variety  of  propagation  conditions.  The  trans¬ 
mitter  was  located  at  the  Earle  Test  Site,  Fort  Monmouth,  New  Jersey,  the 
receiver  site  was  at  an  SRI  field  site  near  Palo  Alto,  California.  The 
great - c i re  I e  distance  between  the  sites  is  1(00  km. 
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B.  ERROR- RATH  AND  S/N  MEASUREMENT 


1.  Discussion  of  thf.  Experiment 

The  purpose  of  this  experiment  was  to  measure  binary  error  rate  and 
S/N  for  an  FSK  data  channel  using  HF  ionospheric  propagation.  The  data 
modem  consisted  of  a  single-channel  AN/FGC-29  transmitter  unit  (FSK 
modulator)  and  the  corresponding  AN/FGC-29  receiver  unit.  The  RF  trans¬ 
mitter  was  a  modified  AN/FRT-51  operated  into  a  horizontal  log-periodic 
antenna.  Dipole  antennas  were  used  for  reception  with  R-390  receivers 
and  CV- 157  sideband  converters.  The  receiving  antennas  were  so  positioned 
that  space-diversity  reception  was  also  available.  The  data  were  simulated 
by  an  alternating  mark-space  sequence  with  a  bit  or  element  length  of 
10  msec  The  error  rate  was  estimated  by  measuring  the  number  of  errors 
detected  at  the  receiver  relative  to  the  number  o 1  binary  digits  transmitted. 
The  S/N  was  measured  by  sequentially  sampling  a  linearly  detected  and 
filtered  version  of  the  signal  within  the  data  channel  with  signal  plus 
noise  (S/N)  present  (FSK  data  tones  on)  and  noise  only  present  (FSK  data 
tone  off).  The  choice  of  obtaining  signal  and  noise  samples  at  adjacent 
times  but  within  the  same  bandwidth— rather  than  the  alternative  of  ob¬ 
taining  signal  and  noise  samples  at  adjacent  bandwidths  but  at  the  same 
time — was  based  on  the  assumption  that  noise,  particularly  narrow-band 
man-made  interference,  would  be  more  correlated  over  a  short  time  interval 
than  over  a  small  bandwidth  separation. 

The  data-gather ing  technique  was  synchronous  in  the  sense  that  per¬ 
tinent  ope  rat i ons  were  pe r f o rmed  within  a  time  accuracy  ol  less  than  a 
millisecond  at  both  the  transmitting  and  receiving  sites  A  30-second 
interval  was  chosen  as  the  basic  data- gather i ng  interval.  During  each 
30  second  interval,  s.x  data  samples  were  taken,  seconds  apart.  Each 
data  sample  consisted  of  a  three-digit,  counter  readout  and  a  digital  volt¬ 
meter  readout.  The  counter  reudout  contained  the  number  of  errors  that 
had  been  detected  in  the  binary  data  stream  since  the  start  of  the 
30-second  interval.  The  voltmeter  reading  represented  a  linearly  detected 
and  filtered  version  of  the  received  signal  contained  within  i  lie  170-eps 
bandpass  of  the  data  channel.  During  the  first  five  samples,  the  data 
tones  were  on  so  that  received  S  N  was  detect.nl.  During  the  sixth  sample, 
the  data  tones  were  off,  and  noise  only  was  detect  *d 

The  data  tones  were  gated  off  at  the  transmit  lei  1 1  mu  the  23th 
second  unt i  l  the  end  of  the  30 -second  period.  Hit.  period  was  long 


18 


enough  lor  the  various  detector  filters  at  the  receiver  to  stabilize  and 
the  noise  level  to  be  sampled  and  recorded.  A  loss  of  the  data  tones 
would  normally  affect  the  receiver  gain  through  the  AGC.  To  avoid  such 
an  effect,  the  AGC  control  voltage  was  derived  exclusively  from  the  carrier 
signal,  which  remained  on  continuously  through  all  da t a - ga ther i ng  periods. 

The  error- total i zing  counter  was  reset  to  zero  at  the  beginning  of 
each  30-second  data  interval.  The  fifth  data  sample  contained  the  fifth 
S/N  voltage  reading  and  the  cumulative  error  count  for  approximately 
23  seconds  of  data  transmission  This  cumulative  error  count  when  divided 
by  the  total  number  of  binary  digits  transmitted  during  that  interval  repre¬ 
sents  the  binary  error  rate.  The  error  count  for  the  sixth  data  sample  was 
discarded  because  the  data  tones  had  been  turned  off  for  the  noise  measure¬ 
ment.  A  binary  error  in  the  received  signal  was  detected  by  comparing  the 
loop  keyer  output  of  the  AN/FGC-29  “combiner”  unit  with  a  locally  generated 
replica  of  the  transmitted  signal  (with  path  and  equipment  delays  taken 
into  consideration). 

2.  Data  Timing  Cvcle 

The  automated  collection  of  data  resulted  in  data  collection  cycles. 

The  time  relationships  within  the  30-second  data  period  are  shown  in  Fig.  ", 
This  basic  cycle  started  at  7  minutes  after  the  hour  and  was  repeated 
eighteen  times  for  a  total  of  9  minutes.  The  tenth  minute  was  used  for  a 
sounder  transmission  at  the  receiver  site,  and  data  transmission  was  not 
attempted  then.  This  minute  provided  an  opportunity  for  the  call  sign  to 


FIG.  7  TIME  RELATIONSHIPS  WITHIN  THE  BASIC  30-SECOND  DATA-CATHERING 
INTERVAL 
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FIG.  8  TIME  RELATIONSHIPS  WITHIN  A  10-MINUTE  DATA-COLLECTION  PERIOD 

be  transmitted,  minor  adjustments  to  be  made  on  the  equipment ,  and  a  time 
reference  to  be  printed  on  the  data  tape.  The  Lime  relationships  within 
a  10-minute  period  are  shown  in  Fig.  8 . 

Because  timing  was  of  utmost  importance  in  this  experiment,  both  the 
transmitting  and  the  receiving  site  used  frequency  standards  that  were 
maintained  in  synchronization  with  receptions  from  WVW  .  These  standards 
drove  digital  clocks  that  delivered  timing  pulses  at  programmed  times  ..o 
generate  the  desired  dat a - gathe r i ng  cycle. 

3.  Description  of  the  Kxperimentm  Setpp 

a.  TRANSMITTING  TERMINAL  EQUIPMENT 

\  block  diagram  of  the  equipment  located  at  the  transmitting 
site  is  shown  in  Fig.  9.  The  frequency  standard  supplied  100-kc  and  i-Mc 
signals.  The  frequency  synthesizer  generated  a  seite>  "I  phase  st  ab  I  e 
frequencies  for  the  mixers  so  that  the  transmittei  iould  operate  at  either 
7.30b  M<  or  1 4 .  3f»0  Me.  The  requirement  t  «>  t  high  phase  stability  was 
associated  with  the  Ikippler  profile  measurement.  wln<b  is  described  in 
more  detail  in  See.  ('.  of  this  chapter 

The  (irenger  Vssoriatrs  t(J  \>  Model  *10  p  >  «  amm.e  r  *  was  also 

driven  b\  the  JOO-kc  frrqurmv  standatd  This  orug i ammo r  is  an  m  curate 
digital  clock  with  a  series  of  programmable  timing  outputs  The  .  lock 
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was  set  to  an  accuracy  of  1  msec  by  using  WWV.  The  programmer  provided 
the  following  timing  signals: 

(1)  Clock  pulses,  occurring  every  10  msec,  which  drove  the 
digital  word  generator 

(2)  A  blocking  signal  which,  when  applied  to  the  AN/FGC-29 
transmitter,  turned  off  the  modulation  so  that  periodic 
in-channel  noise  measurements  could  be  made  at  the  re¬ 
ceiving  sites 

(3)  A  timing  signal  that  initiated  the  call  sign  transmission 
each  time  the  transmitter  was  turned  on 

(4)  A  timing  signal  that  turned  off  the  transmitter  for 
15  seconds  once  every  10  minutes. 


The  exact  time  relationship  between  these  operations  has  been  indicated 
in  Figs.  7  and  8. 


FIG.  9  BLOCK  DIAGRAM  OF  ERROR-RATE  EXPERIMENT  TRANSMITTER  SITE 


The  Mixon  Model  1 12  digital  word  generator 10  is  capable  of  pro¬ 
viding  a  selected  repetitive  binary  sequence  >f  any  I  engt  h  up  l  o  1 2  b  I  n  a  r  V 
digits.  In  this  experiment,  a  simple  alternating  sequence  was  used, 
because  this  was  the  type  of  signal  assumed  for  the  theoretical  error- 
late  rurves.  Therefore,  the  word  generator  in  this  case  was  essentially 
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a  time-8ynthroni zed  square-wave  generator.  The  only  required  input  for 
the  word  generator  was  the  10-msec  clock  pulses  derived  from  the  pro¬ 
grammer.  The  output  binary  data  stream  was  used  to  drive  the  FSK  keyer. 

The  binary  output  signal  from  the  word  generator  consists  of  the 
voltage  levels,  zero,  and  ~5  volts.  The  signal  normally  required  by  the 
AN/FGC-29  transmitter  keyer  is  a  0-  or  60-ma  current  source.  Therefore, 
minor  modifications  were  made  to  the  FSK  keyer  unit  to  make  the  units 
compatible.  The  AN/FGC-29  Channel  1  transmitter  unit  is  essentially  a 
subcarrier  modulator  or  FSK  keyer.  The  center  frequency  for  the  channel 
is  1785  cps  with  MARK  at  1827.5  cps  and  SPACL  >t  1742.5  cps.  The  FSK  out¬ 
put  of  this  unit  was  used  directly  to  modulate  a  si ngt e- si  deband  trans¬ 
mitter.  The  FSK  tone  output  was  blocked  (i.e.,  the  modulation  signal  was 
turned  off)  by  grounding  the  keyer  output  through  a  set  of  relay  contacts 
in  the  programmer  so  that  an  in-channel  noise  voltage  sample  could  be 
taken  at  the  receiver. 

The  AN/FRT- 51  consists  of  a  twi n- sideband  modulator,  a  system  of 
mixers,  an  exciter,  and  the  final  power  amplifier.  The  unit  is  capable  of 
providing  a  power  of  5  kw;  however,  under  these  operating  conditions,  the 
total  output  power  was  limited  to  2  kw,  with  approximately  half  the  power 
in  the  carrier  and  half  in  the  sideband.  The  sideband  in  this  case  was 
the  upper  sideband,  which  consisted  of  the  single  c Manuel  FSK  data  tones. 

All  of  the  signals  used  in  the  modulating  and  mixing  operations,  other  than 
the  modulation  signal  itself,  were  derived  from  the  phase-stable  frequency 
synthesizer  or  its  source,  the  frequency  standard.  A  block  diagram  of  the 
modulator-mixer  chain  (Fig.  10)  shows  how  the  upper  sideband  and  carrier 
were  generated.  Kach  mixer  was  followed  by  a  ^-multiplier  so  that  only 
the  desired  signals  were  retained. 

The  transmitter  site  was  located  at  Fort  Monmouth  (Karle),  New 
Jersey.  V  horizontal  log-periodic  antenna  dire* ted  toward  Palo  Alto, 

Gal  t  torn)  a.  was  used  on  both  7.df>t>  Me  and  l-I.JnO  M* 

h.  RE<  r  I  V  t  Mi- TEHUIN  At  tgttPMKNT 

The  ret  r  i  v  i  ng  - 1  e  rm  i  n  a  1  for  (he  i-tror-rate  experiment  was  located 
at  an  SRI  site  near  Palo  Alto.  C,a  I  »  forn  i  a .  \  block  diagram  of  the  teieivtng 

'V.*t  e»  is  shown  in  fig.  11.  The  s>  »t  >•“’  u<ipi|  two  pa  i  r  *  n!  horizontal  dipole 
antennas — one  pair  tuned  for  each  of  the  operating  fie.joetu  o.,  The  antennas 
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FIG.  10  PHASE-STABLE  MODULATOR-MIXER  CHAIN  IN  MODIFIED  AN  FRT-51 
TRANSMITTER 
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FIG.  11  BLOCK  DIAGRAM  OF  DUAL  ''‘ACE-DIVERSITY  ERROR-RATE 
EXPERIMENT  RECEIVER  SITE 
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of  each  pair  were  physically  separated  by  a  distance  of  approximately 
800  feet.  The  spaced  antennas  were  used  for  diversity  reception. 

A  set  of  spaced  antennas  was  connected  directly  to  separate 
R-390A/URR  radio  receivers 12  tuned  to  the  transmission  frequency.  The 
R-390A/URR  is  a  standard  military  receiver  wi  «-h  a  frequency  range  from 
0.5  to  32  Me.  For  the  demodulation  of  single-sideband  signals,  the  IF 
signal  from  the  R-390/URR  was  used  as  the  input  for  a  single-sideband 
converter. 

The  single-sideband  converters  used  were  standard  military  types 
designated  CV-157/URR.3  The  455-kc  IF  from  the  R-390A/URR  was  mixed  to 
100  kc  in  the  converter.  The  converter  separates  the  transmitted  carrier 
frequency  from  its  sideband  frequencies  and  detects  the  audio  in  the  side¬ 
bands.  The  audio  in  this  case  was  the  subcarrier-modulated  (FSK)  tones. 

The  availability  of  both  the  carrier  and  the  sideband  channels  provided 
several  alternatives  for  the  generation  of  the  AGC  voltage.  The  choice 
of  appropriate  AGC  for  the  experiment  was  considered  at  some  length  and  is 
discussed  in  a  following  paragraph.  The  upper  sideband  was  selected  to 
carry  the  FSK  tones,  and  the  audio  from  this  channel  was  fed  to  the 
Channel  1  input  in  the  AN/FGC-29  receiver  unit. 

Up  to  this  point  in  the  signal  processing,  two  identical  channels 
were  available  from  the  space-diversity  antennas.  The  two  distinct 
channels  were  processed  further  through  the  AN/FGC-29  receiver  unit.  In 
the  receiver  unit,  the  signals  were  appropriately  weighted  for  post¬ 
detection  diversity  combining.  Prior  to  the  weighting,  both  diversity 
channels  had  equal  gain.  Within  each  channel  in  the  AN/FGC-29  receiver, 
the  FSK  signal  v-as  first  passed  through  a  170-cps  bandpass  filter  centered 
at  1785  cps  (Channel  1).  Provisions  were  available  for  adjusting  the  relative 
gain  and  relative  delay  in  both  channels. 

After  the  two  diversity  channels  had  been  amplified,  the  FSK  tones 
were  detected  in  discriminators  in  the  AN/FGC-29  combiner  unit.  The  two  re¬ 
sulting  de  telegraph  signals  were  then  added  to  complete  the  combining 
operation,  The  combined  signal  was  hard-limited  and  drove  the  loop  keyer 
which  provided  the  binary  teletype  output. 

To  determine  t,he  range  of  system  linearity.  th«  relational^  p 
between  sideband  input  S  N  and  output  SN  vas  measured.  This  measurement 
is  significant,  because  a  signal  measurement  (with  sideband  tones  on) 
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followed  by  a  noise  measurement  (with  sideband  tones  off)  was  used  In 
estimate  S/N  with  carrier  continuously  present.  In  a  sense,  this  mea¬ 
surement  is  similar  to  a  measuiement  of  linear  dynamic  range  of  the  system 
for  different  AGC  levels.  A  fixed  noise  level  was  added  at  the  input  of 
the  R-390A/URR  because  the  noise  figure  of  the  receiver  is  a  function  of 
the  AGC  level.  By  deriving  the  AGC  from  the  carrier,  sideband  signal  and 
noise  output  levels  were  measured  by  alternately  turning  on  and  off  the 
sideband  modulation.  The  input  signal  was  fed  to  the  R-390A/URR  followed 
by  the  CV-157/URR  and  the  receiver  unit  in  the  AN/FGC-29.  The  output  S/N 
was  measured  at  a  point  following  the  170-cps  band-pass  filter  in  the 
AN/FGC-29  receiver.  The  AGC  was  developed  in  the  carrier  channel  of  the 
CV-157/URR  and  applied  in  the  R-390A/URR.  The  results  of  the  measurement, 
as  shown  in  Fig.  12,  indicate  that  the  sequential  measuring  technique  esti¬ 
mates  S/Ns  with  reasonable  accuracy  up  to  about  60  db. 


FIG.  12  INPUT-OUTPUT  S/N  CHARACTERISTICS 


The  channel  weighting  for  the  combining  operation  was  performed 
in  the  AN/FGC-29  receiver  unit,  ami  the  post-detection  addition  of  the 
signals  was  performed  in  the  AV  r‘GC-29  combiner  unit.  Before  combining  a 
series  of  voltage#  representing  S  i  N  (di vers 1 1 >  channels),  it  is  general l\ 
advantageous  to  weight  the  individual  voltage#  in  same  manner-  If  the 
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noises  are  mutually  independent,  a  maximization  of  the  output  S/N  is 
achieved  by  weighting  each  incoming  signal  in  proportion  to  its  ratio 
of  signal  voltage  to  noise  power.  On  the  assumption  that  the  noise  is 
constant  and  relatively  small  compared  to  the  signal,  the  combined  out¬ 
put  is  approximately  the  sum  of  the  squares  of  the  inputs.  This  type 
of  combining  is  called  maximal  -  rat  io  or  optimum  combining.14  The  com¬ 
bining  characteristics  of  AN/FGC- 29  receiver  that  were  used  in  this 
experiment  approximated  this  rule;  they  are  shown  in  Fig  13- 

For  the  combiner  to  perform  the  combining  operation  that  it  is 
designed  to  perform,  the  gain  through  the  two  diversity  channels  up  to 
the  point  where  the  weighting  takes  place  should  be  the  same.  To  insure 
this,  both  channels  must  be  balanced  for  equal  gain,  and  any  AGC  control 
voltage  other  than  the  combiner  weighting  must  be  the  same  for  both  channels. 
For  the  experimental  program,  the  channels  were  balanced,  and  a  special  AGC 
selector  circuit  was  used.  Within  each  channel,  an  AGC  voltage  was  developed 
by  placing  the  FUNCTION  switch  on  the  R-390A/URR  in  the  MGC  position.  In 
this  position,  no  AGC  voltage  was  developed  in  the  R-390A/UHR.  On  the 
CV-157/URR,  theAGC  SELECT  switch  was  placed  in  the  CARRIER  position  and  the 
AGC  TIME  switch  in  the  SLOW  position.  Thus,  for  each  channel,  an  AGC  voltage 
was  developed  that  was  slowly  varying  and  roughly  proportional  to  the  carrier 
level  in  that  channel.  The  separate  AGC  voltages  were  combined  in  a  diode 
network  in  which  the  larger  of  the  two  voltages  was  "selected”  and  used  to 
operate  the  AGC  control  circuits  in  both  of  the  R-390A,URRs.  In  this  manner, 
approximately  identical  channel  characteristics  were  achieved  up  to  the 
point  of  combiner  weighting.  The  AGC  voltages  developed  by  the  R-390A  URR 
CV-157/URR  combinations  are  shown  in  Fig.  14. 

The  measurement  of  the  voltages  for  the  computation  of  S/N  was 
made  in  one  of  the  diversity  channels  following  the  170-cps  bandpass  filter 
in  the  AN/FGC- 29  receiver  unit.  This  filter  was  centered  at  ,785  cps.  The 
signal  level  at  the  filter  output  was  estimated  by  passing  the  signal 
through  a  linear  envelope  detector- fi 1  ter  combination  and  monitoring  the 
dc  voltage  on  a  digital  voltmeter.  (The  detector  is  briefly  discussed  in 
Appendix  A.)  The  signal  responses  at  the  input  of  the  voltmeter  to  step 
on  and  off  envelope  changes  of  a  f i xed- f requeue »  sideband  at  the  input 
to  the  R-390A/URR  are  shown  in  Fig.  15.  The  input  output  characteristic 


26 


20  LOG 


FIG.  13  AN/FGC-29  COMBINER  WEIGHTING 
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FIG.  14  AGC  VOLTAGE  v*.  INPUT  SIGNAL  LEVEL 
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FIG.  15  RESPONSE  OF  SIGNAL-LEVEL  DETECTION  SYSTEM 
TO  ON  AND  OFF  STEP  CHANGES  IN  ENVELOPE 

of  the  detection  system  where  the  input  is  considered  to  be  the  input  to 
the  AN/FGC-29  receiver  unit  is  given  in  Fig.  16.  The  nonlinearity  is 
introduced  by  the  AN/FGC-29  receiver  unit.  liy  monitoring  the  detected 
signal  levels  on  the  digital  voltmeter,  an  attempt  was  made  to  maintain 
the  output  level  below  15  volts.  In  this  way.  most  of  the  nonlinear  ef¬ 
fects  in  the  AN/FGC-29  receiver  were  avoided. 

The  digital  voltmeter,  a  Hewlett-Packard  (ll-Pl  Model  I05A,  was 
set  to  range  automatically  and  was  triggered  by  the  G  \  programmer  once 
every  5  seconds.  When  a  voltmeter  reading  was  completed,  it  was  printed 
on  paper  tape  on  the  H-P  Model  560A  digital  recorder.  The  G  \  programmer, 
essentially  a  digital  clock  with  programable  output  timing  pulses,  is 
identical  to  the  unit  used  at  the  transmitter  site.  This  clock  was  main¬ 
tained  on  the  correct  time  ii>  comparison  with  WWV  and  was  driven  b\  a 

frequency  standard  periodically  compared  with  NHV  For  I  minute  in  each 

10  minutes,  a  negative  voltage  was  applied  to  the  digital  voltmeter  input 

through  relay  contacts  in  the  programmer  unit.  During  this  minute,  voltage 
samples  were  taken  once  every  10  seconds.  The  negative  polarity  caused  an 
asterisk  to  be  printed  by  each  voltage  sample.  The  series  ol  si\  astei  isks 
was  usrd  as  a  time  indicator  on  the  papet  data  i ape 
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FIG.  16  SIGNAL-LEVEL  DETECTION  SYSTEM  CHARACTERISTICS 
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The  binary  received  signal  from  the  AN/FGC-29  combiner  unit 
(the  output  from  the  loop  keyer)  was  fed  into  the  Mixon  Model  1032  digital 
woid  analyzer.  10  The  analyzer  was  driven  by  10-msec  clock  pulses  from  the 
programmer  and  generated  a  synchronized  copy  of  the  binary  sequence  gen¬ 
erated  at  the  transmitter.  The  clock  pulses  were  so  adjusted  in  the  pro¬ 
grammer  that  nominal  path  delay  and  any  system  delay  were  taken  into 
consideration.  After  initial  synchronization  of  the  transmitted  signal 
and  the  locally  generated  sequence,  synchronization  was  maintained  at  both 
transmitting  and  the  receiving  terminal  through  use  of  a  common  time 
standard  (WWV).  In  the  analyzer,  the  loop-keying  waveform  was  sampled  at 
the  nominal  center  of  each  received  bit.  The  value  of  each  sample  was 
compared  with  that  of  the  corresponding  bit  of  the  locally  generated 
sequence;  a  pulse  was  generated  each  time  these  values  differed.  Although 
an  error  counter  was  available  in  the  analyze  unit,  an  external  counter 
compatible  with  the  H-P  digital  printer  was  used  (H-P  523B  counter).  As 
has  been  indicated  in  the  timing  diagram  (Fig.  7),  provisions  wer;  made  to 
reset  the  counter  at  the  beginning  of  each  30-second  data- gather i  <ig  interval. 
The  three  least- signi ficant  decades  on  the  counter  were  connected  to  three 
of  the  printer  columns  on  the  H-P  560A  digital  recorder. 

The  raw  data  output  from  the  error-rate  and  S/N  measurement  con¬ 
sisted  of  the  printed  paper-tape  output  from  the  H-P  digital  recorder.  An 
example  of  a  portion  of  a  data  tape — with  notes  to  aid  in  the  interpretation 
of  the  format — is  given  in  Fig.  17.  For  the  exact  timing  of  the  data- 
gathering  operation,  one  should  refer  to  Figs.  7  and  8.  Each  data  tape  was 
dated  .'nd  appropriately  identified  by  the  operators  at  the  receiving  site. 

In  addition,  at  least  once  an  hour,  the  time  to  the  nearest  minute  was 
written  near  a  series  of  six  asterisks  on  the  tape.  In  this  manner  the 
exact  time  associated  with  any  dnta  gathering  interval  could  be  determined. 

A  brief  discussion  of  some  of  the  computational  problems  and 
possible  biases  associated  with  the  estimation  of  power  and  power  ratios  is 
given  in  Appendix  B.  The  results  indicate  that  because  of  the  statistical 
characteristics  of  the  noise  and  the  use  of  a  linear  envelope  detector,  the 
S/N  estimate  may  be  2  or  3  db  high.  The  '.  >pagation  of  error  in  the  power 
estimates  into  the  error  in  S/N  expressed  in  decibels  is  examined  in 
Appendix  C.  If  one  excludes  the  bias  that  has  been  indicated  and  admits  a 
percentage  error  of  25  percent  in  the  power  estimates,  the  resulting  S  \ 
estimate  has  an  error  bound  of  +2  db.  The  complete  error  bound  is  also 
gi ven  in  Appendi x  C. 
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FIG  17  AN  ANNOTATED  EXAMPLE  OF  THE  PAPER-TAPE  DATA  OUTPUT 
FOR  THE  ERROR-RATE  AND  $  N  MEASUREMENTS 
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C.  DOPPLER  PROFILE  MEASUREMENT 


1.  Discussion  of  tke  Experiment 

The  purpose  of  this  experiment  was  to  provide  a  measure  of  the 
f requency - spreading  characteristics  of  the  communication  channel.  One 
should  refer  to  Daly,7  Gal  lager,  15  Kai  1  ath,  16  or  Bello17  for  an  explanation 
of  some  of  the  details  motivating  the  measurement.  The  normalized  Doppler 
profile  q D( \),  and  such  measures  as  the  second  central  moment  of  qt,(k)  have 
been  chosen  to  display  the  frequency-spreading  characteristic  of  the 
channel.  The  normalized  Doppler  profile  is  defined  in  terms  of  the  channel 
scattering  function,  Sy(k,r),  by 


jSv(k,T)dT 

<*D^  ^  ffSy(k,  T)drdk 


The  channel  scattering  function  may  be  interpreted  as  the  density  of  power 
scattered  by  the  channel  as  a  function  of  Doppler  shift  and  time  delay. 

An  alternate  expression  for  q D(k)  is 

jRH(a,0)e-'27rKad<x 
q  (k)  ~  - 

jjRH(a,0)e-'2l,Kad<ldk 
where  fiw(a,0)  is  given  by 

Rg{a, 0)  -  ♦  a./)} 

The  function  is  called  the  time-variant  channel  transfer  function, 

it  is  a  complex  random  variable  assumed  to  be  wide-sense  stationary  in 
both  time  and  frequency.  The  symbol  E  is  used  to  indicate  the  mathematical 
expectation.  The  quantity  //(!,/)  may  be  defined  as  the  complex  envelope  of 
the  channel  response  to  a  C'lf  tone  relative  to  the  frequency  of  the  trans 
mitted  tone.  Thus  the  goal  of  the  experiment  *a*  to  measure  the  coaiplex 
envelope  of  the  channel  response  to  •  sinusoidal  transmission  over  a  given 
period  of  time.  Given  the  coaiplex  envelope,  the  computations  required  for 
estimating  the  Doppler  profile  have  been  considered  in  detail  by  Daly.* 
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The  measurement  required  for  the  estimation  of  the  Doppler  profile 
has  been  reduced  to  a  measurement  of  the  complex  envelope  (real  envelope 
and  phase)  of  the  signal  received  over  the  channel  when  the  transmitted 
signal  is  a  CJf  tone.  The  complex  envelope  must  be  measured  relative  to 
the  frequency  of  the  transmitted  signal;  thus  high  phase  stability  is  re¬ 
quired  at  both  the  transmitter  and  the  receiver. 

2.  Description  of  the  Experimental  Setup 

With  a  few  minor  modifications,  the  experimental  setup  for  this  mea* 
surement  is  the  same  as  that  described  by  Shepherd. 19  The  experimental  set¬ 
up  was  made  more  flexible  by  introducing  operation  at  14.360  Me,  as  well 
as  at  7.366  Me.  The  transmitter  portion  of  the  system  is  the  same  as  that 
used  in  the  error-rate  S/N  measurement  (block  diagrams  are  given  in  Figs.  9 
and  10).  As  has  been  mentioned,  half  the  power  is  retained  in  the  carrier, 
and  half  is  used  for  the  transmission  of  sideband  information.  The  carrier 
is  synthesized  from  a  phase-stable  source  and  is  in  fact,  the  transmitted 
CW  tone  used  for  the  Doppler  profile  measurement.  Since  the  carrier  is 
synthesized,  its  stability  depends  upon  the  stability  of  the  frequency 
standard.  The  standard  used  at  both  the  transmitting  u-’d  the  receiving  site 
was  an  H-P  Model  103  AR,  which  has  a  stability  of  5  parts  in  1010per  day. 

Although  the  error-rate  experiment  and  the  Doppler  profile  experiment 
shared  the  same  transmitter  unit,  separate  receivers  were  requited  at  the 
receiver  site.  However,  one  of  the  dual - di versi ty  receiving  antennas  was 
used  in  common  for  both  experiments.  A  block  diagram  of  the  phase-stable 
receiving  system  is  shown  in  Fig.  18.  The  crystal  oscillators  were  avail¬ 
able  for  calibration  of  the  tape-recorded  values  of  phase  and  amplitude 
through  the  entire  system.  The  R-390A/URR  receiver  was  modified  to  accept 
synthesized  local  oscillator  signals,  and  an  IF  of  459  kc  was  generated. 

For  operation  at  14.36  rather  than  7.366  Me.  the  second  local  oscillator 
frequency  was  changed  from  27  to  16.994  Me;  the  i, reiver  was  tuned  to  the 
higher  frequency;  and  the  antenna  was  changed  The  IF  output  from  the  re 
ceiver  was  mixed  to  13  kc  and  bandpass  filtered  to  a  passbund  of  200  cp*. 

The  output  from  tie  filter  went  to  the  phase  detector  and  the  envelope  de¬ 
tector.  The  output  of  the  envelope  detector  was  re. orded  as  the  level  of 
the  received  signal  and  was  used  as  the  \GC  control  u  llage  m  the 
R-39QA  DRR  receiver.  The  output  level  as  a  function  of  th«-  input  has  been 
given  by  Shepherd  **  and  is  shown  in  Fig  19  Ihe  phene  drieitnr  «l*<‘  1»,»* 
ber,  vl<  s<  rtbed  in  detail  b>  .Shepherd.  Riirflv  the  phase  of  the  received 
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signal,  nominally  at  13  kc,  is  compared  with  the  phase  of  a  locally  syn¬ 
thesized  phase-stable  13-kc  reference  signal.  The  relative  phase  of  the 
received  signal  is  recorded  by  recording  the  detected  phase  angle  modulo 
180  degrees  and  all  positive  and  negative  transitions  through  an  integral 
multiple  of  180  degrees.  Thus,  over  any  spec  i  f  i  ed  interval  of  time:,  both 
the  envelope  an  1  phase  of  the  received  signal  are  recorded.  The  time  as 
received  on  a  HW  receiver  and  any  operator  notes  are  also  recorded  on  the 
magnetic  tape. 

Data  were  collected  for  the  Doppler  profile  measurement  from  1800  to 
2300  GMT  on  14.36  Me  and  from  0000  to  0500  GMT  on  7.366  Me.  During  these 
intervals,  data  were  collected  on  chart  paper  continuously,  and  magnetic 
tape  recordings  were  made  starting  at  6,  26,  and  46  minutes  after  the  hour. 
Each  of  these  recordings  lasted  for  a  period  of  6  minutes. 

3.  Output  Data  and  Pbepbocess ing 

The  output,  data  for  the  Doppler  profile  measurement  consist  of  a  six- 
channel  tape  recording  representing  phase,  envelope,  time,  and  a  voice  re¬ 
cording  of  the  operator’s  comments.  An  additional  supplementary  two- 
channel  Sanborn  chart-paper  recording  of  envelope  and  phase  modulo 
180  degrees  was  also  available.  This  recording  is  identical  to  the  magnetic 
tape  recording  of  two  of  the  channels.  An  example  of  this  chart  is  given  in 
Fig  20.  The  scale  for  the  signal  envelope  record  is  roughly  logarithmic. 

A  Doppler  shift  in  addition  to  the  spread  is  indicated  by  the  phase  having 
a  nominal  rate  of  change  in  one  direction  (from  the  bottom  of  the  record  to 
the  top).  The  chart  recordings  were  monitored  at  the  site  to  ensure  proper 
equipment  operation  and  to  provide  data  for  a  preliminary  evaluation. 

The  six-channel  magnetic- tape  data  required  preprocessing.  The  phase 
and  envelope  signals  were  passed  through  constant- time-delay  low-pass 
filters  wit!  a  cutoff  frequency  of  approximately  50  tps.  finally,  the  data 
were  converted  into  digital  lorm  and  recorded  in  the  appropriate  format  on 
digital  tape. 

D.  TIME -DELAY  PROFILE  MEASUREMENT 
|.  Discus  si  on  or  the  Experiment 

The  purpose  of  this  subexperiment  was  to  pro  idr  a  measure  of  the 
t ime - spread) ng  character i st i is  of  the  channel.  The  significant  effects 
of  time  spreading  of  the  commun i t at i on  channel  at  HE  have  been  recogntied 
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for  many  years.  Multipath  differential  delays  have  been  measured  by 
obi ique- incidence  ionospheric  sounders  for  some  time.  Recent  developments 
in  channel  modeling  suggest  alternate  measures  which  include  consideration 
of  the  relative  amplitude  of  various  modes.  7,15,16,17  In  a  manner  similar  to 
that  indicated  for  the  Doppler  profile  measurement  (Sec.  C.l)  the  nor¬ 
malized  delay  profile,  qT(r)  may  defined  as 

fS¥(k,r)dk 

qr{T)  ‘  fSS,iKr)drdK  ' 

The  normalized  delay  profile  and  twice  the  second  central  moment  of  q^r) 
have  been  chosen  to  represent  the  time- spreading  characteristics  of  the 
channel . 

An  alternate  expression7  upon  which  a  measurement  technique  may  be 
developed  is 


£{|z(t)|2}  »  /T(r)£{|x(t  -  r) | 2}dr 

where  z(t)  and  x(t )  are  the  complex  envelopes  of  received  and  transmitted 
waveforms,  respectively,  and  T(r)  is  the  time-delay  profile  before  normali¬ 
zation.  This  expression  suggests  the  transmission  of  a  pulse  that  is  short 
relative  to  the  resolution  desired  in  the  estimation  of  T(r),  Signals  of 
other  types  are  also  useful  as  long  as  the  quantity  x(t  ~  r) | 2  performs 
the  desired  sifting  operation  on  T(r) .  An  estimate  of  T(t)  with  a  series 
of  short  pulses  is  given  by 


j  I  x  (  t  ~  r )  (  2  dr 

and  the  normalized  estimate  is  given  by 
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To  achieve  the  desired  resolution  in  the  real  envelope  of  the 

transmitted  signal,  |*(t)l  was  chosen  to  be  a  sounder  short  pulse  of 
100  jusec.  A  special  modification  was  incorporated  in  the  G/A  sounder 
transmitter  to  provide  singl e - frequency  operation.  At  the  chosen  frequency, 
a  series  of  fifteen  short-pulse  transmissions  were  made  to  provide  confidence 
in  the  estimate.  Modifications  were  also  made  in  the  G/A  sounder  receiver 
unit  to  permit  repetitive  reception  at  a  single  frequency.  The  video  signal 
from  the  sounder  occurring  in  the  window  of  10  msec  was  supplied  to  a  unit 
designated  Sounder-Computer  Link  (SCL)  MK  I.  This  unit  used  the  sounder 
receiver  timing  signal  and  provided  analog- to-digi tal  (A-D)  conversion  of 
the  received  signals.  The  digital  output  was  recorded  on  digital  tape, 
along  with  timing  and  other  identifying  information.  The  digital  tape 
format  was  designed  to  be  directly  compatible  with  tape  input  requirements 
for  IBM  computers.  Further  processing,  including  refined  detection  based 
on  the  series  of  samples,  was  carried  out  in  the  computer.  This  process 
is  described  in  Sec.  D  of  Chapter  IV. 

2.  Description  of  yhe  Experimental  Setup 

The  sounder  system  used  to  measure  the  channel  characteristics  con¬ 
sisted  of  a  G/A  Model  901  transmitter  unit  located  at  Earle,  New  Jersey 
and  a  G/A  Model  903  receiver  unit  located  near  Palo  Alto,  California.  For 
this  experiment,  certain  of  the  automatic  features  were  de-activated  by 
relays  to  allow  repetitive  operation  on  a  single  frequency  on  both  the 
transmitter  and  the  receiver  units.  A  block  diagram  of  the  time-delay 
profile  measurement  system  is  given  in  Fig.  21.  The  frequency  standards 
supplied  both  the  sounder  and  the  programmer;  they  had  a  frequency  stability 
of  5  parts  in  10*°per  dav.  The  frequency  accuracy  was  maintained  through 
periodic  calibration  to  NllA,  The  programmers  provided  accurate  timing  for 
both  the  transmission  and  reception  of  the  sounding  pulses. 

The  sounder  was  operated  at  a  frequency  near  the  frequencies  of  opera¬ 
tion  of  the  Dopp ler’ pro fi  le  experiment  — 7  .  366  Me  and  14.360  Me  —  that  was 
also  relatively  free  of  interference.  This  generally  resulted  in  a  choice 
of  the  G/A  sounder  Frequency  Bands  \-32  (7.25  Me)  or  V  33  (7.35  Me)  and 
B*  32  (14.5  Me)  or  B-33  (14.7  Me).  A  brief  description  of  general  sounder 
characteristics  has  been  given  by  Shepherd.  ***  The  modifications  made  de¬ 
activated  the  f requency • st eppi ng  features  through  the  use  of  relays.  Both 
sounder  transmitter  and  the  sounder  receiver  were  modified  to  provide 
single- frequency  operation  as  at.  optional  operational  mode. 
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FIG.  21  TIME-DELAY-PROFILE  MEASUREMENT  SYSTEM 
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The  sintole- frequency  soundings  at  14  Me  were  taken  once  every 
5  minutes,  starting  at  1802  GMT  and  continuing  until  2257  GMT.  *  On  7  Me, 
soundings  were  made  from  2357  GMT  until  0457  GMT.  Operation  thus  concided 
with  the  collection  of  data  on  all  of  the  other  experiments.  Each  single¬ 
frequency  sounding  was  affected  by  the  transmission  of  i  series  of  fifteen 
100-yusec  30-kw  pulses  with  a  2-second  spacing  between  pulses.  Thus  approxi¬ 
mately  30  seconds  was  required  for  the  complete  sounding.  At  the  receiver, 
the  detected  RF  signal  was  used  to  observe  the  response  of  the  channel  to 
each  of  the  fifteen  pulses.  The  detected  signal  from  the  sounder  was 
supplied  to  the  SCL  for  further  processing.  The  SCL  performed  the  A-D 
conversion  and  interface  operations  between  the  sounder  and  a  computer. 

This  unit  has  been  described  in  detail  by  Sifford.20  During  a  predeter¬ 
mined  time  interval  of  10.2  msec,  samples  were  taken  once  every  50  /isec 
for  a  total  of  204  samples.  Each  sample  was  quantized  into  one  of  the 
thirty-two  levels  and  encoded  into  a  five-bit  code  word.  The  digital 
samples  with  parity  checking  bits  were  recorded  on  the  digital  tape.  The 
SCL  provided  an  identification  preamble  in  digital  form,  including  the 
time  that  the  data  were  taken.  The  data  format  on  the  tape  was  designed 
to  be  used  as  direct  input  to  further  computer  processing  programs. 


IV  DATA  PROCESSING 


A.  INTRODUCTION 

Data  from  the  experiments  were  recorded  on  various  media  and  in 
unique  formats.  To  extract  meaningful  results  and  to  present  them  ef¬ 
ficiently,  high-speed  computers  were  employed  in  a  data-processing  phase. 

A  flow  chart  illustrating  the  basic  components  of  this  processing  phase 
is  shown  in  Fig.  22.  The  end  result  was  a  determination  of  confident  ex¬ 
perimental  estimates  of  the  correlation  of  measured  S/Ns  with  measured 
binary  error  rate  for  an  FSK  communications  system  on  an  HF  channel. 
Measured  time  delay  and  frequency  spreading  of  the  channel  and  spaced  an¬ 
tenna  correlation  were  parameters. 

Accumulated  binary  errors  detected  over  a  30-second  interval,  to¬ 
gether  with  simultaneous  samples  of  received  signal  and  noise  voltages, 
were  recorded  on  adding-machine  tape  from  a  digital  printer.  Time-delay- 
spread  d^ta  were  recorded  in  digital  form  on  magnetic  tape  directly  from 
the  output  of  an  ionospheric  sounder  receiver.  Frequency- spread  or  Doppler 
data  were  recorded  in  analog  form  on  magnetic  tape  in  the  form  of  CW  amp¬ 
litude  and  phase  measurements.  Correlation  of  signals  received  on  a  pair 
of  spaced  antennas  was  recorded  on  chart  paper  directly  from  the  output 
of  a  special  correlation  instrument.  The  antenna- correl ation  data  pro¬ 
cessing  is  further  discussed  in  Chap.  VII. 

Partial  piocessing  of  the  30-second  S/N  and  bit-error  data  was  done 
on  a  Burroughs  B5000  computer  at  the  SRI  Computer  Center.  The  time-delay- 
spread  processing  was  run  on  the  IBM  7090  computer  at  Stanford  University. 
The  CW  analog  amplitude  and  phase  measurements  wer*  converted  to  digital 
form  and  recorded  on  magnetic  tape  on  an  A-D  converter  and  COC  160A  com¬ 
puter  at  SRI.  Those  data  were  than  processed  to  extract  the  channel 
Doppler  profiles  on  the  7090  computer  at  -‘anford  University.  The  results 
of  these  three  processing  program*  provided  the  input  to  the  final  clas¬ 
sification  and  summation  program  Ibis  program,  which  tabulated  the  S/N 
and  bit-error  rates  for  various  measured  channel  spread  conditions,  was 
run  on  the  B5000  at  SRI . 
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FIG.  22  DATA-PROCESSING  FLOW  DIAGRAM 


B.  PRE-EDITING  OF  S/N  AND  BIT  ERROR  DATA 


The  S/N  and  bit-error  data  were  manually  edited  before  being  punched 
onto  IBM  cards.  This  editing  consisted  primarily  of  grouping  the  data  for 
each  3C-sccond  interval  into  the  five  signal  readings  and  the  one  noise 
reading  together  with  the  accumulated  bit-erroi  totals  during  these  samples 
In  addition,  the  legibility  and  completeness  of  each  data  block  were 
checked,  and  data  suspected  of  being  faulty  were  deleted  from  the  tape. 

The  five  signal-sample  readings  and  the  nois^  sample,  together  with 
the  bit-error  totals  at  these  sample  times,  were  punched  onto  IBM  cards 
from  this  edited  tape.  The  format  for  this  card  is  shown  in  Fig.  23. 

Each  card  contained  the  data  for  one  30-second  data  period.  Time  header 
cards  were  inserted  in  the  deck  to  indicate  the  starting  time  of  a  sequence 
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FIG.  23  BASIC  30-SECOND  SIGNAL,  NOISE.  AND  BIT-ERROR  DATA  CARD  FORMAT 
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C.  S/N  AND  HIT  ERROR  PREPROCESSING 


I.  S/N  Estimation 


Five  samples  of  the  received  S  +  N,  taken  five  seconds  apart,  were 
recorded  from  the  output  of  a  linear  recti fi er/ low-pass  f t I  ter  combina- 
t i on .  It  is  shown  in  Appendix  B  that  if  the  signal  and  noise  are  inde¬ 
pendent  zero-mean  processes,  then  the  best  estimate  of  the  total  S  +  N 
is  proportional  to  the  square  of  the  average  of  the  five  sample  readings. 
Furthermore,  this  estimate  represents  the  sum  of  the  power  in  the  signal 
and  noise  components. 

The  noise-power  estimate  was  made  from  the  average  nf  the  two  noise 
samples  immediately  before  and  after  the  corresponding  five  S  +  N  readings. 
In  the  event  that  a  valid  noise  reading  was  not  taken  immediately  before 
or  after  each  of  the  five  signal  readings,  a  search  prodecure  was  ini¬ 
tiated  to  obtain  two  usable  noise  samples.  The  search  procedure,  based 
on  cursory  observations  on  the  autocorrelation  of  the  noise,  was  as  fol¬ 
lows.  If  a  valid  noise  sample  was  not  available  immediately  before  (or 
after)  a  30-second  signal  data  block,  then  the  noise  reading  for  the  pre¬ 
ceding  (or  following)  interval  was  used.  This  reading  was  weighted  by 
80  percent,  however,  in  calculating  the  average  noise  estimate  with  the 
other  sample.  However,  if  the  preceding  (or  following)  noise  sample  was 
also  missing,  the  next  preceding  (or  following)  30-second  noise  reading 
in  time  was  used,  and  an  additional  86-percent  de-emphasis  was  applied 
to  this  reading  in  determining  the  average.  The  search  for  two  valid 
noise  samples  was  continued  in  this  manner.  If  two  noise  sam  Pi  e  s  were 
not  found  within  four  30-second  periods  removed  from  the  signal  data, 
then  all  data  for  this  30-second  period  were  rejected. 

The  formula  used  by  the  processor  for  calculating  the  noise  average 

Wtt.s 


-4.v  ( r  > 


(1.3  -  0.  2KM)\iT  -  Mf>  *  <1.0  -  0  ,'.kP)\lT  ►  KP) 


where 

T  is  an  integer  designating  a  particular  30  second  block  of 
signal  samples  followed  by  a  noise  sample 

KM  is  the  lowest  integer  from  *>ne  to  four  *hrre  the  noise 
sample  at  f  *  KM  is  valid 
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KP  is  the  lowest  integer  from  zero  to  three  where  the  noise 
sample  at  T  +  KP  is  valid, 

N(T)  is  th'-  noise  value  at  time  T. 

The  procedure  wan  used  at  least  once  every  10  minutes,  since  no 
data  were  recorded  during  the  sixth  m;nute  of  each  10-minute  data- 
gathering  period.  For  the  first  30-second  data  period  of  the  seventh 
minute,  the  procedure  bridged  the  ^ap  by  using  the  noise  reading  taken  at 
the  end  of  the  second  30-second  block  of  data  during  the  fifth  minute  to 
estimate  the  noise  for  the  first  block  f  data  during  the  seventh  minute. 

The  S/N  in  decibels  was  calculated  for  each  30-second  data  period  as 
fo 1 1 ow3 : 


m  ■ 10  u*t3 

whe  - 

_ _ 2 

Ps  Ps  +  Pn  ~  Pn  P(s  +  n)  Z(s  +  n ) 

Pn  Pn  Pn  — ~2 

Zn 

Z(s  H  n)  is  the  average  of  the  fi-.e  S  +  N  samples 

Zn  is  the  average  of  the  two  noise  samples. 

If  the  ratio,  PsPn  was  lea  th-a  one,  then  the  5  V  was  assigned  a  value 

oi  -100  db. 

2.  Validation  of  Bit  Error  Count 

It  was  theoretically  possible  to  accumulate  as  many  as  2,251  bit 
errors  during  each  30-second  data- gathering  period  because  *>nlv  the 
first  three  digits  of  the  error  counter  were  recorded.  f*’s  -aunter  oc¬ 
casionally  overflowed  during  a  data  period  In  sense  an  overflow  in  the 
counter,  the  processing  pro«r?»s;  w..»  tu*iru»ted  to  monitor  the  error  count 
during  each  , 5-second  period.  Whenever  the  accumulated  error  r«u»  t  ’  any 
second  recurding  period  was  Im  than  the  preceding  count,  an  additional 
l . 000  e, ror*  were  added  t»  the  recorded  !"t«l . 

On  numerous  occasions,  excessive  errors  writ  retarded  ;«rtng  the  ex¬ 
periment.  even  with  a  high  S  ?i  These  errors  were  principally  caused  1- 
difficult*  with  the  automatic  frequency  control  and  by  occasional  loss  oj 
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bit  synchronism.  Even  though  an  attempt  was  made  to  detect  such  condi¬ 
tions  and  reject  the  resulting  data  as  they  occurred  during  the  experi¬ 
ment,  man^  incidences  were  not  detected.  Therefore,  a  criterion  was  es¬ 
tablished  in  this  program  to  reject  data  if  the  total  errors  made  during 
a  data  period  were  judged  to  be  excessive  on  the  basis  of  the  average  S/N 
measured  during  the  data  period.  The  30- second  data  were  rejected  if  the 
product  of  the  maximum  number  of  errors  detected  in  any  five-second  period 
when  data  were  transmitted  times  the  S/N  in  decibels  was  greater  than 
2,000.  For  example,  assume  the  S/N  was  20  db.  The  number  of  errors  made 
during  any  5-second  interval  then  must,  not  be  greater  than  100  for  the 
data  to  be  accepted.  The  test  was  not  applied,  however,  if  the  measured 
S/N  was  less  than  ten  db. 

3.  Preprocessed  Output 

Averaged  S/N  in  decibels  together  with  corresponding  total  bit  er¬ 
rors  for  each  accepted  30- second  time  period  were  punched  onto  IBM  cards 
by  the  preprocessing  program.  The  format  for  this  card  is  shown  in 
Fig.  24.  Each  data  card  included  the  date  and  time  for  convenience  in 


correlating  later  with  channel- spread  measurements  and  antenna  correla¬ 
tions.  A  listing  of  these  values  for  each  30 -second  period  was  also  pro 
duced  on  the  page  printer,  together  with  notations  indicating  where  and 
why  data  were  rejected  (see  Table  I).  The  notations  of  rejected  data 
permitted  the  cause  of  the  data  rejection  to  be  investigated  and,  if  pos 
sible,  corrected.  In  many  cases,  mispunched  or  misordered  input  cards 
were  detected  by  the  various  da t a- rej ec ti on  tests. 
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Ta  b 1 e  1 

SAMPLE  OF  PREPROCESSE1)  30- SE(X)NI)  DATA 


TIME 

AVSIG* 

ERRORS 

A  VN** 

S/N 

BEK 

DATE  IS 

1106 

4530 

5, 418  +  01 

9 

2.308-03 

43.7 

4*008-03 

4600 

NU  DATA 

4630 

NU  DATA 

4700 

1.118+02 

15 

2.888-03 

45.9 

6.668-03 

4730 

ABNORMAL 

ERRORS  FUR 

S/N.  DATA 

REJECTED 

4730 

1.678+02 

465 

4.368-03 

45.8 

2,078-01 

480  0 

DATA  OUT 

OF  SYNCH 

4830 

DATA  UUT 

UF  SYNCH 

4900 

DATA  UUT 

UF  SYNCH 

4930 

DATA  UUT 

UF  SYNCH 

5000 

DATA  UUT 

UF  SYNCH 

5030 

DATA  UUT 

UF  SYNCH 

5100  INSUFFICIENT  DA  1 A  TU  CUMPuTE  AVERAGE  NOISE 


5130 

6.518+01 

23 

1.338-03 

46.9 

1.028-02 

5200 

1.238+02 

4 

2.708-03 

46.6 

1.788-03 

5230 

3.348+01 

5 

8.378-03 

36.0 

2.228-03 

5300 

5.668+01 

2 

5.558-03 

40.1 

8.868-04 

5330 

7.348+01 

1 

2.358-03 

44.9 

4.448* 04 

5400 

8.548+01 

0 

3.028-03 

44.5 

0.008+00 

5430 

1.068+02 

4 

3.848-03 

44.4 

1.788-03 

5500 

6.208+01 

43 

3.608-03 

42,4 

1.918-02 

5530 

4.598+01 

19 

3.428-03 

41.3 

8.448-03 

5600 

NU  DATA 

5630 

NU  DATA 

5700 

1.338+02 

1 

4*478-03 

44.7 

4.448-04 

5730 

1.168+02 

5 

3.728-03 

44.9 

2.228-03 

5800 

1.048+02 

8 

2.558-03 

46.1 

3.558-03 

5830 

1.048+02 

1 

2.558-03 

46.1 

4.448-04 

5900 

1.908+02 

0 

6.088-03 

44.9 

0.008+00 

5930 

7.228+01 

9 

6.978-03 

40.2 

4.008-03 

10000 

1.208+02 

2 

2.408-03 

47.0 

8.888-04 

10030 

1.038+02 

0 

1.378-03 

48.8 

0.008>00 

10100 

6.508+01 

1 

2.268-03 

44.6 

4.448-04 

10130 

7.438+01 

2 

2.868-03 

44.1 

8.888-04 

10200 

7,448+01 

8 

2.268-03 

45.2 

3-558*03 

10230 

1.438+02 

2 

3.258-03 

46.4 

8.888-04 

10300 

9,278+01 

3 

3.728-03 

44,0 

1.338-03 

10330 

8.418+01 

11 

2.768-03 

44.8 

4.898-03 

10400 

1.048+02 

6 

4,428-03 

43.7 

2.678-03 

10430 

1,648+02 

0 

4,498-03 

45.6 

0. 008+00 

10500 

1.908+02 

6 

3.428-03 

47,4 

2.678-03 

10530 

4.158+01 

0 

2,558-03 

42.1 

0.008+00 

10600 

NU  DATA 

10630 

NU  DMA 

10700 

8.538+01 

2 

2.948-03 

44.6 

8.868-04 

107  30 

1.428+02 

11 

2.928-03 

46.9 

4.898-03 

10800 

1.938+02 

0 

1.988-03 

49,9 

0.008+00 

10830 

2.068+02 

5 

3.368-03 

47,9 

2.228-03 

10900 

1,368+04 

0 

4.768-03 

64.6 

0.008+00 

*  AVSIG 

*  Average  of 

five  signal  readings. 

AVN 

-  Average  of 

two  noise  readings. 
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TIME 

AVSIfi 

ERRORS 

AVN 

S/N 

BE* 

DATE  IS  1106 

10930 

1 *640+02 

2 

4.166-03 

46.0 

8 .886*04 

11000 

1*466402 

0 

3.3?>*-OJ 

46.4 

0. 006  +  00 

11030 

9.966401 

4 

2.836-03 

45.5 

1.786-03 

11100 

1.336402 

0 

3.318-03 

46,1 

0. 006+00 

11130 

1.666402 

0 

7.066-03 

43,8 

0.0U6+00 

11200 

9.726401 

) 

6.246-03 

41 ,9 

0.006+00 

11230 

5.466401 

0 

3.976-03 

41.4 

0.006+00 

11300 

6.956401 

1 

3.976-03 

42.4 

4.446*04 

11330 

5.096401 

6 

3.546-03 

41.6 

2.676-03 

11400 

3.766401 

2 

4.566-03 

39.2 

8.886-04 

11430 

1.666401 

5 

5.118-03 

35.2 

2.226-03 

11500 

5,546401 

2 

8.936403 

37.9 

8.886-04 

11530 

8.746401 

2 

7.406-03 

40,7 

8.886-04 

11600 

NU  DATA 

11630 

NU  data 

11  TOO 

1.166402 

0 

8.216-03 

41.5 

0.006+00 

11730 

4,206401 

11 

5.046-03 

39,2 

4.896*03 

11600 

4.916401 

1 

1,768-03 

44,4 

4.446*04 

11630 

6.296401 

1 

3.546-03 

42,5 

4.446-04 

11000 

2.666401 

0 

4.166-03 

38.1 

0.0U6+Q0 

11930 

4,646401 

2 

2.556-03 

42.6 

8.886-04 

12000 

4,636401 

0 

5.636-03 

39.2 

0.006+00 

12030 

9, *96+01 

3 

9.606-03 

40.2 

1.336-03 

12100 

6.356401 

11 

8.286-03 

40.0 

4.896-03 

12130 

5,736401 

5 

6.486-03 

39.5 

2.226-03 

12200 

4.736401 

0 

5.636-03 

39.2 

0. 006+00 

12230 

1.096402 

10 

7.576-03 

41.6 

4.446-03 

12300 

1.046402 

0 

1.056-02 

40.0 

0.006+00 

12330 

1.296402 

1 

1.106-02 

40.7 

4.446*04 

12400 

3.576401 

2 

3.426-02 

30.2 

8.886-04 

12430 

1.756401 

0 

3.594-02 

26.9 

0.006+00 

12500 

9.888401 

2 

2.876-02 

35.4 

8.886-04 

12530 

1,346402 

7 

2.346-02 

37.6 

3.116-03 

12600 

NO  DATA 

12630 

NU  DATA 

12700 

1.656402 

0 

2.446-02 

38.8 

0. 008+00 

12730 

9.096401 

0 

2.214-02 

36.2 

0.006+00 

12600 

6,946401 

3 

7.836-03 

40.6 

1.336-03 

12630 

5.206401 

3 

8.376-03 

37.9 

1.338-03 

12900 

7.946401 

3 

7,146-03 

40.5 

1.336-03 

12930 

1.926402 

0 

4.766-03 

46.0 

0.006+00 

13000 

1.226402 

19 

4.226-03 

44.6 

8.446-03 

13030 

9.996401 

4 

1.236-02 

39,1 

1.786-03 

13100 

5.296401 

3 

1.536-02 

35.4 

1.336-03 

13130 

2,026402 

6 

1.216-02 

42,2 

2.876-03 

13200 

1.216402 

5 

1.078-02 

40.5 

2.226-03 

13230 

6.486+01 

1 

4.766-03 

41.3 

4.446-04 

13300 

6. 126  +  01 

0 

4.626-03 

41.2 

0. 006+00 

» 
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TIME 

AVSIG 

ERRORS 

A  VN 

S/N 

BEK 

DATE  IS  : 
13330 

1106 

9,146*01 

1 

7.926-03 

40,6 

4.446-04 

13400 

6,156*01 

0 

2.046-02 

34.6 

0.006*00 

13430 

4.656*01 

0 

1.406-02 

35.2 

0.008*00 

13900 

6 • 506*01 

0 

5.556-03 

40.7 

ft. 006*00 

13530 

7.918*01 

0 

8.016-03 

39.9 

0.006*00 

13600 
13630 
13  TOO 

NO  DATA 
NU  DATA 
4,256*02 

1 

4.358-02 

39.9 

4.446-04 

13730 

3,006*02 

0 

3.266-02 

39.6 

0.006*00 

13600 

6,506*01 

0 

3.768-03 

42.3 

0.006*00 

13630 

8.066*01 

0 

4.766-03 

42.3 

0.006*00 

13900 

1.226*02 

0 

1.246-02 

39.9 

0.006*00 

13930 

7.286*01 

0 

1.196-02 

37.9 

0.008*00 

14U00 

7.696*01 

5 

5.708-03 

41,4 

2.226-03 

14030 

7.686*01 

0 

1.296-02 

37.9 

0.006*00 

14100 

5.596*01 

3 

1.468-02 

35.8 

1.336-03 

14130 

3.476*01 

16 

5.166-03 

36.3 

8. 006*03 

14200 

4,066*01 

6 

3.026-03 

41.3 

2.678-03 

14230 

7.206*01 

0 

2.406-03 

44,6 

0.006*00 

14300 

7.546*01 

6 

6.726-03 

40.5 

3.558-03 

14330 

6.446*01 

33 

1.568-02 

36,2 

1.476-02 

1440ft 

8.586*01 

0 

1.068-02 

39.1 

0.006*00 

14430 

6.676*01 

18 

6.166-03 

40.3 

8.0U8-03 

14500 

7.046*01 

0 

6.728-03 

40.2 

0.008*00 

14530 

4.736*01 

5 

7.316-03 

36.1 

2.226-03 

14600 

14630 

14700 

NU  DATA 
NU  data 
6.006*01 

6 

9.906-03 

37.8 

2.678-03 

14730 

6.646*01 

15 

1.506-02 

36.6 

6.668-03 

14600 

7.256*01 

6 

1.566-02 

36,7 

2.678-03 

14030 

abnormal 

ERRORS  FUR 

S/N#  data  hejected 

14630 

1,118*02 

92 

1.466-02 

38.8 

4.096-02 

14900 

4.426*01 

50 

4.496-02 

29.9 

2.226-02 

14930 

2.728*01 

0 

3.7o8-02 

26.6 

0.006*00 

150C0 

7.688*01 

10 

8.266-03 

39.7 

4.446-03 

15030 

3.618*01 

0 

6.976-0 3 

37.1 

0.006*00 

15100 

4.688*01 

0 

5.116-03 

39,6 

0.008*00 

15130 

ABNURMAL 

ERRORS  FuH 

S/N#  OATA  1 

REJECTED 

1513ft 

6.328*01 

615 

2.506-03 

44.0 

2.738-01 

15200 

4 . 946*0 1 

23 

7.666-03 

36.1 

1.026-02 

15230 

1.138*02 

0 

7.466-03 

41,8 

0.006*0o 

15300 

7,726*01 

l 

7.226-03 

40,3 

4,446-04 

15330 

6,638*01 

6 

1.006-02 

36.2 

2.676-03 

15400 

7.658*01 

4 

2.616-03 

44.3 

1.786-03 

15430 

4,648*01 

5 

1.526-03 

44.8 

2.226*03 

15500 

6.688*01 

0 

2.976-03 

43.5 

0.006*00 

15530 

7.448*01 

0 

4.036-03 

42.7 

0,006*00 

15600 

15630 

15700 

NU  OAIA 
NU  OATA 
7.206*01 

7 

2.868-03 

44,0 

3.116-03 

51 
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TIME 

A9SIQ 

ERRORS 

A  VN 

S/N 

BEK 

0A1E  IS 
15/30 

1106 

1. 268*02 

9 

2.769-03 

46.6 

4.008-03 

15600 

3.649401 

2 

2.219-03 

42*4 

8.869*04 

15830 

1.759402 

0 

2.219-03 

4*  ,0 

0.008*00 

15900 

9.349401 

5 

2.979-03 

45*0 

2.228*03 

15930 

5,209401 

0 

3.609-03 

41.6 

0.008+00 

20000 

1.039402 

3 

3*369-03 

44.9 

1.338*03 

20030 

9.769401 

0 

3.549-03 

44,4 

0.009+00 

20100 

6.069401 

0 

6.069-03 

40.0 

0.009+00 

20130 

9.409401 

6 

4.909-03 

42.6 

2.678*03 

20200 

5.259401 

0 

2.609-03 

43.1 

0. 009+00 

20230 

1.199402 

0 

2.609-03 

46,6 

0.008+00 

20300 

3.929401 

11 

1.989-03 

43,0 

4.898*03 

20330 

5.539401 

0 

3.369-03 

42.2 

0. 009+00 

20400 

1.169402 

0 

5.169-03 

43,6 

0.009+00 

20430 

1.479402 

3 

3.149-03 

46.7 

1.339*03 

20500 

5.569401 

0 

2.709-03 

43.1 

0.009+00 

20530 

1.179402 

0 

3.609-03 

45.1 

0.009+00 

20600 

20630 

20/00 

NO  DA! A 
NO  DATA 
1,679402 

4 

5.069-03 

45.7 

1.789*03 

2073O 

ABNORMAL 

ERRORS  FUR 

S/N.  DATA  REJFCTED 

20730 

1.169402 

245 

3.029-03 

45.6 

1.099*01 

20600 

7.259401 

6 

3.609-03 

43.0 

2.679*03 

20830 

ABNORMAL 

ERRORS  FUR 

S/N*  OATA  REJECTED 

20630 

1.049402 

968 

6. 4*9-03 

42.1 

4.309*01 

20900 

ABNORMAL 

ERRORS  fur 

S/N*  OATA 

REJECTED 

20900 

7 . 749401 

82 

2.338-02 

35.2 

3.649*02 

20V30 

3.609401 

0 

2.338-02 

31.9 

0.009+00 

21000 

6.079401 

16 

7.838-03 

40.  \ 

8. 009*03 

21030 

1.319402 

13 

?.  838-03 

42.2 

5.789*03 

21100 

4,689401 

9 

4.429-03 

40.2 

4 .009*03 

21130 

8.199401 

0 

5.118-03 

42.0 

0.009*00 

21200 

8,049401 

18 

5.408-03 

40,5 

8.008-03 

21230 

7.749401 

9 

4.629-03 

42.2 

4.008-03 

21300 

5,439401 

1 

3. 764-03 

41.6 

4.449*04 

21330 

9.699401 

0 

2.609-03 

45,6 

0.009+00 

21400 

1.179402 

0 

2.929-03 

46.0 

0.009+00 

21430 

4,649401 

5 

3.069-03 

42.0 

2.229-03 

21500 

1.359402 

0 

3.428-03 

46.0 

0 .009+00 

21330 

3.699401 

0 

3.764-03 

39.9 

0.009+00 

21800 

21630 

21700 

NO  OATA 
NO  DATA 
1,109402 

0 

3.668-03 

44.6 

0 .009+00 

21730 

4.769401 

6 

2.409-03 

43.0 

2.678-03 

21600 

2.789401 

6 

2.029-03 

41.4 

2.679-03 

21630 

A,  289401 

4 

4.699-03 

39.6 

1,788-03 

21900 

1,109402 

7 

1 .559-02 

38.5 

3.119-03 

219)0 

2.US9402 

0 

1 .0*8-02 

42.7 

0 .009+00 

22000 

abnormal 

ERRORS  FuH 

S/N*  DATA 

REJECTED 

22000 

6,369401 

126 

2.309-03 

45.7 

5.609-02 

22030 

1,199402 

11 

1.618-03 

48.2 

4.8V9-03 

22100 

9.369401 

0 

3.258-03 

44.6 

0.009*00 

TIME 
DA1E  IS 

A  V  $  1 1» 
1106 

erhuks 

AVN 

S/N 

BE* 

22130 

3« 169401 

21 

6.816*03 

36,7 

9.336-03 

22200 

1.406402 

0 

6.166-03 

44.6 

0.01)6*00 

22230 

5.236401 

0 

6.016-03 

39.* 

0.006*00 

22300 

3.306401 

4 

1.466-02 

34.1 

1.706-03 

22330 

1.496402 

2 

1.616-02 

39.7 

6.866-04 

22400 

9,636401 

2 

9.416-03 

40.2 

8,366-04 

22430 

1.546402 

0 

7.066-03 

43.4 

0.006*00 

22500 

3.726401 

2 

2.H66-03 

41.1 

8.686*04 

22530 

22600 

22630 

6.986401 

No  DATA 

NU  DATA 

0 

1.726-03 

47.2 

0.006*00 

22/00 

1.306402 

0 

2.506-03 

46.6 

0.006*00 

22/30 

5. *46401 

6 

2.026-03 

44.0 

2.676*03 

22600 

5.036401 

3 

6.016-03 

39.2 

1.336*03 

22630 

1.226402 

5 

7.316-03 

42.2 

2.226-03 

22900 

7.696401 

0 

3.256-03 

43.9 

0.006*00 

22930 

4.716401 

0 

2.216-03 

43.3 

0.006*00 

23000 

6.696401 

0 

9.616-04 

46,6 

0.006*00 

23030 

5.346401 

0 

6.506-04 

49.1 

0.006*00 

23100 

7.706401 

0 

1.446-03 

47.3 

0.006*00 

23130 

6.166401 

1 

1,446-03 

47.5 

4.4*6*04 

23200 

1.436402 

0 

1.336-03 

50.3 

0.006*00 

23230 

1.116402 

0 

4.906-03 

43.5 

0.006*00 

23300 

3.146401 

0 

6.066-03 

37.1 

0.006*00 

23330 

6.116401 

5 

2.666-03 

43.3 

2.226-03 

23400 

1.166402 

0 

3.346-03 

44,9 

0.006*00 

23430 

6.356401 

0 

5.656-03 

40,4 

0.006*00 

23500 

7.223401 

6 

5.706-03 

41.0 

3.556-03 

23530 

23600 

23630 

1.136402 

NU  DATA 

NU  OATA 

0 

6.466-03 

42.4 

0.006*00 

23/00 

9.026401 

0 

9.936-03 

39,6 

0.006*00 

23730 

1,426402 

9 

1.246-02 

40,6 

4,006-03 

23600 

7,656401 

0 

7.7*6-03 

39,9 

0 . 006*00 

23630 

1.966402 

0 

2.706-03 

46,6 

0. 006*00 

23900 

7  r  *16*01 

0 

1,856-03 

46,0 

0.006*00 

23930 

3.546401 

23 

2.866-03 

40.9 

1.026-02 

24000 

1.266402 

9 

2.266-03 

47.5 

4.006-03 

24030 

1.006402 

0 

9.006-04 

50,5 

0.006*00 

24100 

8.606401 

0 

6.506-04 

51.3 

0.006*00 

24130 

2.636401 

0 

1.566-03 

42.3 

0.006*00 

24200 

1.926403 

13 

2.506-03 

56.6 

5.766-03 

24230 

1.336402 

0 

2.976-03 

46.5 

0.006*00 

24300 

1.356402 

0 

2,356-03 

47.6 

0.006*00 

24330 

6.956401 

4 

1,066-03 

48,2 

1.796-03 

24400 

7,536401 

0 

9.616-04 

46,9 

0.006*00 

24430 

4.626401 

3 

1.416-03 

45.2 

1.336-03 

24500 

1.256402 

0 

1.266-03 

50.0 

0.006*00 

time 

4VS16 

ERRORS 

DATE  IS  1106 

34530 

3.296401 

2 

24600 

NU  0414 

24610 

NU  04T4 

24700 

5.516401 

0 

34710 

6.536401 

0 

34600 

4.066401 

3 

24610 

2.056402 

6 

24900 

7.456401 

0 

24910 

1.236402 

0 

25001, 

6.526401 

4 

25010 

5.VA6401 

5 

25100 

6,556401 

15 

25130 

6.V16401 

7 

25200 

7.126401 

0 

25210 

7.106401 

3 

2530 0 

1.906404 

0 

25330 

4.366401 

0 

25400 

6.016401 

0 

25430 

2.076401 

13 

25500 

5.606401 

7 

25530 

3.626401 

31 

25600 

NU  04  f 4 

25630 

NU  0414 

25700 

1.456402 

0 

25730 

6.276401 

0 

25600 

3.156402 

0 

25630 

1,296*0? 

3 

25900 

3.576401 

6 

25930 

3.706401 

0 

30000 

1,296402 

0 

30030 

4.626401 

0 

30100 

4.456401 

0 

30130 

7,066401 

0 

30200 

1,076402 

0 

30230 

4,036401 

6 

30300 

7,676401 

3 

30330 

1.166402 

7 

30400 

1.356402 

1 

30430 

1,466402 

0 

30500 

1,046402 

0 

30530 

6,506401 

0 

30400  NU  04 f 4 

30630  NU  0414 

30700 

1.406402 

0 

30730 

3.9064U2 

0 

30600 

1.616402 

0 

30630 

6.666401 

0 

30900 

3,269402 

0 

AVN 

S/N 

BC<< 

1.026-03 

45.1 

6,669-04 

7.436-04 

46*7 

0.009*00 

1.726-03 

45.6 

0.006*00 

2.616-03 

41.6 

1.339-03 

2.026-03 

50.1 

2.676-03 

3.426-03 

43*4 

0.006*00 

4.169-03 

44,7 

0.006*00 

1.619-03 

46.7 

1.766*03 

2.669-03 

43.2 

2.226-03 

3.069-03 

44.4 

6.666*03 

3.199-03 

43,4 

3.116-03 

6.569-03 

40.4 

0.006*00 

2.359-03 

44,6 

1.336-03 

2.359-03 

69.1 

0.006*00 

3.769-03 

40.6 

0.0U6*00 

3.609-03 

43.5 

0.006*00 

3.319-03 

36.0 

5.766-03 

9. 009-04 

47.9 

3.116-03 

1,169-03 

45.2 

1.386-02 

1.666-03 

48,9 

0.6u6*00 

2.306-03 

44.3 

0,006*00 

7.836-03 

46.0 

0.006*00 

9.316-03 

41.4 

1.336-03 

5.116-03 

38.4 

3.556-03 

3,146-03 

40.7 

0.006*00 

8,416-04 

51.8 

0.006*00 

1.376-03 

45.5 

0.006*00 

1.566-03 

44.5 

0.006*00 

9.306-04 

48.8 

0.006*00 

2.166-03 

46.9 

0.006*00 

1.766-03 

43.6 

2.676*03 

4.846-04 

52.1 

1.336-03 

4,416-04 

54.2 

3.116-03 

2.356-03 

47.6 

4 .446*04 

2.816-03 

47.2 

0.006*00 

1.166-03 

49.5 

0.006*00 

1.266-03 

48.3 

0.006*00 

1.966-03 

48.5 

0. 006*00 

1.686-03 

53.7 

0.006*00 

6.506-04 

53.9 

0.006*00 

1.196-03 

46.7 

0.006*00 

3.786-03 

49.4 

0.006*00 

54 


A  binary  error-rate  estimate  (HUH)  was  also  ea  I  rii  I  a  toil  and  printed 
oil  this  listing  for  each  30-seeond  data  period.  However,  beennse  of  the 
limited  number  of  data  bits  transmitted  during  the  period  (i.c.,  2,251) 
the  confidence  in  the  estimate  is  quite  low.  No  further  use  was  mude  of 
this  figure. 

D.  CHANNEL  TIME- DELAY- SPREAD  PROCESSING 

1.  Raw  Data  Format 

Time-delay-spread  estimations  ofthe  channel  were  made  from  digitized 
recordings  of  the  signal  received  from  a  remote  pulsed  transmitter.  The 
raw  data  consisted  of  a  time  series  of  204  amplitude  samples,  spaced 
50  fi sec  apart,  of  a  receiver  envelope  detector  where  each  sample  was 
quantized  in  to  32  ampl i tude  levels  and  recorded  digi tal ly  on  magneti c  tape. 
The  timing  for  the  series  of  samples  was  synchronized  with  the  transmis¬ 
sion  of  a  lC0-/zsec  pulse  at  the  remote  transmitter.  A  time-delay-spread 
estimate  was  made  from  an  ensemble  of  fifteen  such  recordings  spaced  two 
seconds  apart. 

2.  Standard  Devi ation- to-Mean  Test 

An  ensemble  mean  and  standard  deviation  was  calculated  for  each  of 
the  204  time  samples  over  the  fifteen  recordings.  To  eliminate  extraneous 
impulse  noise  and  interference  from  the  data,  a  limit  was  placed  on  the 
ratio  of  the  standard  deviation  to  ensemble  mean.  If  this  ratio  was 
greater  than  2.0,  the  eniemble  mean  for  the  time  sample  was  set  to  zero. 
The  standard  deviation- to-jnean  threshold  of  2.0  was  chosen  so  that  normal 
Rayleigh  fading  signals  would  not  be  rejected.  (See  Appendix  A.) 

3.  Noise-Based  Threshold  Test 

A  noise  power  estimate  was  made  from  the  average  of  the  square  of  the 
first  twelve  and  last  twelve  ensemble  mean  time  samples.  The  synchronized 
timing  of  samples  with  respect  to  the  transmissions  insured  that  no  actual 
signal  from  the  remote  transmitter  waa  received  at  these  times.  This  noise 
average  was  multiplied  by  a  constant  and  used  as  a  threshold  ov  which  the 
received  power,  represented  by  the  ensemble  mean  squared  (EMS),  in  each 
time  sample  was  compared.  The  multiplying  factor  used  was  3.0.  If  the 
EMS  for  a  time  sample  did  not  exceed  or  equal  the  noise-based  threshold, 
then  the  EMS  for  that  sample  time  waa  aet  to  aero.  !f  it  equaled  or  ex¬ 
ceeded  this  threshold,  it  was  unaltered. 
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4.  Adjacent  Time-Cell  Test 

To  further  discriminate  against  impulse  interference,  it  was  required 
that  at  least  two  successive  time  samples  be  above  the  noise-based  thresh¬ 
old  in  order  to  be  classified  as  signal.  For  each  EMS  found  above  the 
noise-based  threshold,  a  test  of  the  EMS  in  adjacent  time  cells  was  made. 
If  neither  rf  the  two  adjacent  EMS  samples  was  above  the  noise-based 
threshold,  then  the  EMS  for  this  time  sample  was  set  to  zero.  If  either 
of  the  two  adjacent  EMS  exceeded  the  noia;  threshold,  it  was  accepted  as 
true  signal  and  totalized. 

5.  S/N  Test 

A  final  test  on  the  confidence  of  the  lime-delay-spread  data  was  made 
<  the  basis  of  the  total  received  S/N.  This  test  required  the  total  S/N 
to  be  10  db  or  greater. . 

At  this  stage  in  the  process,  a  plot  of  the  EMS  samples  — normal  i  *ed 
to  the  total  power— as  a  function  of  the  time  samples  was  sometimes  pro¬ 
duced.  These  time-delay  profiles  were  used  for  detailed  examination  of 
the  channel  time-delay  dispersion.  Examples  of  these  plots  are  shown  and 
discussed  in  Chapter  V. 

6.  Time-Delay  Spread  Calculations 

From  the  series  of  EMS  at  the  various  time  samples,  a  mean  time,  de¬ 
fined  as  the  average  time  delay  of  the  path,  was  found.  Finally,  a  sec¬ 
ond  moment  about  this  mean-time  sample  was  determined.  Twice  the  square 
root  of  this  moment,  converted  to  milliseconds,  was  used  as  a  measure  of 
the  time-delay  spread  of  the  channel.  The  exact  computational  form  was 


where  r  u  the  second-moment  time-delay  spread  of  the  channel  in  milli¬ 
seconds 


P,  is  the  EMS  at  time  sample  i 
f  is  the  mean  time-delay  sample. 


Tlii  s  spread  calculation  together  with  the  time  of  the  recording  in  GMT 
was  punched  on  IBM  cards.  The  format  for  this  card  is  shown  in  Fi®.  25. 
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FIG.  25  TIME-DELAY -SPREAD  DATA  CARD  FORMAT 

These  cards  were  used  as  input  to  the  final  classification  processing  of 
S/N  and  bit-error  data  by  channel  spread  conditions. 

E.  CHANNEL  FREQUENCY- SPREAD  PROCESSING 

Doppler- frequency-spread  estimates  for  the  channel  were  made  from 
recordings  of  the  phase  and  amplitude  of  a  CW  transmission  over  the  chan¬ 
nel.  The  instrumentation  for  this  experiment,  previously  developed  by 
another  SRI  project,6 has  been  described  in  Chapter  III.  Calibrated  volt¬ 
ages  representing  phase  and  amplitude  of  the  received  CW  signal  were  re¬ 
corded  in  analog  form  for  6-minute  intervals  every  20  minutes  on  magnetic 
tape.  These  amplitude  and  phase  records  were  then  digitised  through 
an  Ad com  A-D  converter  and  CDC  160A  computer  combination  and  recorded  in 
digital  form  on  magnetic  tape.  Amplitude  and  pbaae  samples  were  taken 
every  10-mset.  «nd  each  sample  was  quantised  into,  one  ef  the  256  discrete 
levels. 

T,e  digitised  phase  and  amplitude  recordinga  were  further  proceeeed 
on  an  IHM  7090  computer  to  extrect  •  Doppler- frequency  profile.  Daly  has 
described  thie  program.1  The  procedure  eonaieta  of  the  following  beaic 
step*: 

tl)  ^termination  of  the  quadrature  component*  of  the  received 
signal  from  the  phase  end  amplitude  samples 

(2)  Calculation  of  the  autocorrel n ion  end  cross-correlation 
functions  associated  with  tht  *  components 


«- 

U  * 


(3)  Estimation  of  the  Doppler  profile  by  a  Fourier  transform  cf 
these  autocorrelations. 

Each  Doppler-profi le  calculation  was  based  on  1,200  samples  of  am* 
plitude  and  phase  spaced  at  100-msec  intervals.  This  yielded  a  10-cps 
bandwidth  spectrum  with  a  power  estimate  every  l/p0  cps.  The  results  at 
this  stage  were  often  plotted;  examples  of  these  plots  are  shown  in  Chap¬ 
ter  V. 

A  single-parameter  estimate  of  the  energy  spreading  in  frequency  by 
the  channel  was  based  on  the  variance  of  the  frequency  spread  by  a  tech¬ 
nique  similar  t»  that  used  in  the  time-delay-spread  calculation:  a  mean 
Doppler-shi ft  frequency  was  found  and  a  second  moment  about  this  mean  was 
determined.  The  square  root  of  this  moment,  multiplied  by  two,  was  used 
as  the  frequency- spread  parameter  for  this  time. 

The  frequency- spread  calculations  together  with  the  date  and  time 
were  punched  onto  IBM  cards  in  a  format  identical  to  that  shown  in  Fig.  25 
for  the  time-delay-spread  calculations. 

F.  TABLES  FOR  CHANNEL- SPREAD  CLASSIFICATIONS 

The  classification  of  channel-spread  conditions  was  based  on  two 
considerations: 

(1)  Theoretical  differences  in  error  rates  at  high  S/N  would 
be  detectable  for  different  channel-spread  conditions. 

(2)  Sufficient  data  would  be  collected  in  each  channel  group¬ 
ing  to  form  confident  estimates  of  S/Ns  and  error  rates, 
particularly  at  high  S/Ns. 

To  provide  flexibility  in  these  classification  decisions,  the  group¬ 
ings  were  defined  to  the  processing  program  in  the  form  of  a  two-dimensional 
matrix  as  shown  in  Fig.  2b.  Each  row  was  subjectively  defined  by  a  time- 
delay-apread  range  for  the  channel,  and  each  column  was  assigned  a 
frequency-spread  range.  Next,  each  element  in  the  matrix  was  assigned  a 
table  number  defining  classifications  into  which  data  with  these  channel- 
spread  conditions  would  be  grouped.  As  the  amount  of  data  in  each  classi¬ 
fication  became  apparent,  the  classifications  were  easily  altered  by  simply 
redefining  the  column  or  row  limits  or  changing  the  elements  *he  matrix. 

Each  table  in  the  proeeaaor  consisted  of  two  80-elei^nt  arrays,  thie 
of  the  arrays  contained  the  total  number  of  30-second  data  periods  found 
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at  each  S/N  from  0  to  80  db  in  1-db  in¬ 
tervals.  All  S/Ns  less  than  0  db  were 
grouped  with  the  0-db  entries,  while 
S/Ns  greater  than  80  db  were  added  to 
those  in  the  80-db  interval.  The  other 
array  contained  accumulated  bit  errors 
at  the  corresponding  S/N  intervals. 

G.  CLASSIFICATION  OF  S/N 
AND  BIT- ERROR  DATA 

Channel  time-delay-  and  frequency- 
spread  data  as  a  function  of  date  and 
time  were  generated  by  the  time-delay- 
and  frequency- spread  processing  pro¬ 
grams.  These  data,  on  punched  cards, 
were  arrange  d  chronologically  into  two 
decks.,  one  for  time-delay-spread  data 
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FIG.  26  DEFINITION  MATRIX 
FOR  CLASSIFICATION 
OF  TIME-DELAY-  AND 
DOPPLER-SPREAD  DATA 


and  the  other  for  frequency- spread  data.  Each  deck  was  entered  into  the 
cl  ssification  program  where  chronological  time  intervals  were  established 
during  which  a  particular  channel- spread  classification  was  assigned.  The 
time  interval  was  determined  from  the  recorded  time-delay-  and  frequency- 
spread  values.  The  start  of  the  time  interval  was  determined  by  calcu¬ 
lating  one-half  the  difference  between  recorded  time  and  the  recorded 
time  for  the  preceding  spread  measurement.  The  end  cf  the  time  interval 
was  determined  by  calculating  one-half  the  difference  between  that  re¬ 
corded  time  and  the  recorded  time  for  the  following  spread  measurement. 

In  this  manner,  any  missing  data  were  simply  bridged  over  by  the  bracket¬ 
ing  spread  measurements.  Each  time  interval  was  then  associated  with  a 
particular  classification  table  b>  matching  the  time-  and  frequency- 
spread  measurements  tor  the  time  interval  with  the  channel  -  spread  classi¬ 
fication  matrix  A  listing  of  these  time  intervals  and  the  table  numbers 
assigned  to  each  interval  is  shown  in  Appendix  E  according  to  the  defini¬ 
tions  of  tables  given  in  Fig.  26. 


H.  S/N  ERROR- RATE  DATA 

All  of  in*-  30-second  S/N  and  error  data  were  then  grouped  according 
to  channel  -  spread  condition*  observed  at  the  time  the  data  were  taken. 
Finally,  the  data  : n  each  classification  or  table  were  processed  further 
to  produce  a*  many  confident  error-rate  and  correaponding  S  N  ratio  points 
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as  possible  Considerations  used  to  establish 
mates  are  discussed  in  Appendix  D.  The  result 
require  that  the  total  observed  errors  for  S/N 
cation  equal  or  exceed  100.  The  computational 
contents  of  each  table  was  as  follows: 


confi  dent 
ing  proces 
intervals 
procedure 


error-rate-  esti- 
sing  rule  was  to 
in  each  cl assi f  i  - 
for  emptying  the 


Cl) 


The  total  number  of  errors  recorded  for  an  S/N  of  0  db 
was  noted. 


{2)  If  this  number  was  greater  than  100,  a  binary-error- rate 
estimate  was  made  by  dividing  the  number  of  errors  by  the 
total  bits  received  at  this  S/N.  The  total  bits  received 
were  determined  by  the  product  of  the  number  of  30-second 
data  entries  at  this  S/N  and  2,251,  the  number  of  bits 
transmitted  during  each  30-second  data  period.  The  pro¬ 
cedure  was  then  repeated  for  the  next  higher  S/N  decibel 
interval . 


(3)  If  the  total  number  of  errors  in  the  0-db  S/N  inter¬ 
val  was  less  than  100,  the  errors  at  the  next  higher 
S/N  interval  (i.e.,1  db)  were  added  to  the  total.  If 
the  total  errors  were  then  greater  than  100,  a  weighted 
average  S/N  was  calculated,  based  on  the  number  of  en¬ 
tries  at  each  oi  the  two  S/N  intervals.  The  correspond¬ 
ing  binary  error  rote  was  determined  from  the  total  er¬ 
rors  and  the  total  number  of  entries  at  both  the  0-  and 
i-db  S/iN  intervals. 

(4)  If  the  total  number  of  errors  in  the  0  and  1-db 
intervals  was  still  less  than  100,  Step  (3)  above  was 
repeated,  taking  the  errors  and  entries  at  the  next, 
higher  S/N  until  the  error  total  exceeded  100.  At  this 
tine,  a  weighted  average  S/N  was  calculated  along  with 
a  binary  error  rate. 

15)  After  the  highest  S/N  interval  had  been  entered,  the 

average  S/N  and  binary  error  rate  were  determined,  even 
if  the  total  errors  did  not  exceed  100. 


Both  punched  cards,  for  uae  in  data  plotting,  and  a  printout  of  the  aver¬ 
age  S/N  and  corresponding  error-rate  calculations  were  produced.  Print¬ 
outs  of  all  results,  sorted  into  various  channel  -  spread  classifications, 
are  shown  in  Appendix  F. 
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V  PRESENTATION  AND  DISCUSSION  OF  RESULTS 


A.  CHANNEL  MEASUREMENTS 

1.  Ray-Path  Geometry 

The  character  of  radio  signals  propagated  via  the  ionosphere  varies 
greatly  with  time  of  day  and  season  and  with  transmission  frequency.  Two 
reflecting  layers  are  responsible  for  the  return  of  HF  transmissions  on 
the  path  from  Fort  Monmouth,  New  Jersey  to  Palo  Alto,  California.  The 
principal  layer  is  the  F  layer  nominally  located  at  a  height  of  300  km 
above  the  earth.  The  other  s  the  E  layer  at  about  100-km  height. 

On  this  path,  a  great-circle  distance  of  4100  km,  it  is  normal  for 
signals  to  arrive  at  the  receiver  in  many  combinations  of  reflections  from 
these  layers  and  the  ground.  Each  of  these  propagation  paths  or  modes  has 
a  characteristic  propagation  time  depending  on  the  total  distance  traveled 
and  the  group  velocity  of  propagation  along  the  path. 

Movements  and  turbulence  of  these  layers  cause  a  Doppler- frequency 
shift  and  spread  of  the  received  signal  from  thut  transmitted.  Each  mode 
contributes  a  Doppler  component  to  the  total  received  signal.  This  com¬ 
ponent  may  be  a  fixed  frequency  shift  of  all  transmission  frequencies  or  a 
smearing  oi  the  received  signal  over  a  band  of  frequencies. 

2.  lONOGRAMS 

Propagation  time  for  each  mode  can  change  drastically  with  transmission 
frequency.  Characteristic  patterns  of  propagation  time  with  transmission 
frequency  are  often  used  to  identify  experimentally  the  various  propagation 
modes  present  on  a  path.  Obi ique- incidence  ionosphere  sounders  measure  and 
display  this  time  del  ay- frequency  characteristic  directly;  such  a  sounder  has 
been  used  on  this  particular  path  for  many  years.  Typical  innograms  produced 
on  this  path  during  the  error- rate  experiment  are  shown  in  Fig.  27.  The  MOF 
of  each  mode  is  defined  as  the  highest  frequency  propagated  by  that  mode,  ,»s 
determined  by  the  ionogram,  while  the  LOF  i#  the  lowest  detectable  frequency 
for  the  mode. 
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FIG.  27  OBLIQUE-INCIDENCE  IONOGRAMS  —  FORT  MONMOUTH  TO  PALO  ALTO 


3.  Characteristics  at  14  Me 


From  1800  to  2300  GMT  (1200  to  1700  local  time  at  midpath),  the  error- 
rate  experiment  was  operated  at  14  Me.  The  predominant  mode  at  these 
times  is  normally  a  strong  two-hop  F-layer  reflection  (i.e.,  2 F).  The  high 
ray  from  this  mode  (2 FH)  often  propagates  also.  The  signal  from  the  2FH  is 
often  unresolvable  in  time  delay  from  the  normal  2 F,  but  can  arrive  up  to 
0.5  msec  after  the  2 F  low  ray.  At  times,  usually  when  the  MOF  of  the  2F  is 
lower  than  14  Me,  it  is  possible  to  receive  signals  from  the  one-hop  F-layer 
high  ray.  However,  the  normal  IF  low-ray  node  is  usually  not  present  on 
this  path  because  of  the  earth’s  curvature. 

From  path  geometry,  it  is  impossible  for  propagation  to  occur  by  a 
single  hop  from  the  E  layer.  However,  it  is  possible  for  propagation  to 
occur  by  two  or  three  hops  from  this  layer.  A  two-hop  E  mode  is  somewhat 
rare  during  the  daylight  hours  at  14  Me,  but  a  three-hop  E  mode  is  more 
common.  The  3 E  mode  arrives  approximately  0.2  msec  before  the  arrival  of 
the  normal  2 F  mode. 

A  very  common  and  strong  propagation  mode  from  1800  to  2300,  at  14  Me, 
is  the  so-called  N  mode.  This  is  a  two-hop  combination  mode  consisting  of 
a  one-hop  reflection  from  the  F  layer  and  a  one-hop  reflection  from  the  £ 
layer.  The  jV-mode  signal  arrives  a  small  fraction  of  a  millisecond  after 
the  3 E  mode,  if  present,  or  about  0.2  msec  before  the  2 F  mode.  A  typical 
range  of  time-delay  difference  between  the  first  and  last  arriving  mode  at 
this  frequency  and  time  of  day  is  from  0  to  1.0  msec. 

Background  noise  is  either  of  local  origin  or  is  propagated  from  great 
distances  by  the  ionosphere.  Local  noise  is  generally  man-made  interference 
Propagated  noise  is  either  atmospheric  noise  or  interference. 

The  heavy  fl-region  absorption  during  the  day  acts  as  a  substantial  at¬ 
tenuator  to  all  propagated  signals,  particularly  those  at  low^jjf.|r*quencies. 
The  signal  strength  of  each  mode  is  affected  by  its  total  path  length  and 
the  distance  traveled  in  the  abrorbing  D  region.  Variation  in  absorption 
introduces  long-term  fading.  Of  more  practical  importance,  however,  is  the 
effect  of  the  antenna  gain  pattern  on  the  strength  of  the  signals  received 
(and  transmitted)  from  the  various  modes  which  arrive  (and  leave)  it  differ¬ 
ent  vertical  angles.  Short-term  variations  of  mode  signal  strength  are 
also  caused  by  polarisation  fading  and  fading  within  resolvable  modes. 
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4.  Characteristics  at  7  Me 
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The  experiment  wes  operated  at  7  Me  after  sunset  on  the  path  (1600- 
2000,  local  midpath).  Initially,  at  this  frequency  the  path  is  charac¬ 
terized  by  the  propagation  of  higher-order  modes  (i.e.,  3 F,  4 F,  and  5 F) . 
However,  aa  the  evening  wears  on,  the  electron  density  of  the  lower  iono¬ 
sphere  diminishes,  and  the  absorption  in  the  D  region  lessens.  The  initial 
effect,  therefore,  is  a  strengthening  of  the  higher-order  modes.  But 
shortly  thereafter,  the  electron  densi  ty  of  the F region  diminishes,  and  the 
MOF  of  all  modes  begins  to  drop.  As  a  result,  the  number  of  modes  propa¬ 
gating  at  this  frequency  is  large  at  first,  may  increase  slightly,  and  then 
decreases  as  the  MOF  of  the  higher-order  modes  drop  through  7  Me.  The 
two-  and  three-hop  E  modes  are  prevalent  at  this  frequency  and  time  of  day. 

A  typical  range  of  time-delay  differences  from  the  modes  propagated  at 
this  frequency  and  time  of  day  is  0.2  to  2.0  msec. 

Loss  of  the  D  region  at  night  also  causes  the  propagated  noise  and 
interference  to  increase  significantly.  The  geographical  origin  of  this 
noise  is  more  widely  spread  than  during  the  day. 

Other  modes  are  possible  on  this  path,  particularly  off-path  modes  due 
to  ionospheric  anomalies.  These  modes,  if  present,  arrive  with  appreciably 
longer  time  delays. 

5.  Time-Delay  Profiles 

Accurate  error-rate  performance  estimations  of  an  HF  communications 
systems  require  a  rather  detailed  picture  of  the  scattering  of  the  signal 
energy  by  the  channel  in  time  and  frequency.  A  time-delay  scattering 
picture  was  obtained  by  actually  measuring  and  displaying  an  average  of 
several  impulse  responses  of  the  channel.  The  measurement  technique  has 
been  described  in  Sec.  D.  2  of  Chapter  III.  Selected  examples  of  time- 
delay  profiles  at  14  and  7  Me  produced  directly  from  the  data  processor 
output  by  a  CalComp  piotter  are  shown  in  Figs.  28  through  33.  The  time- 
delay  scale  has  been  shifted  so  that  zero  time  delay  corresponds  to  the 
average  delay  of  the  path.  The  power  density,  measured  at  50-yisec  inter¬ 
vals,  is  normalised  to  the  total  received  power. 

Twice  the  second-moment  spread,  2^T.  calculated  for  each  of  these 
profiles,  is  elso  shown.  This  spread  measure  tends  to  be  less  than  the 
time  difference  »;om  the  first  to  last  arriving  mode,  as  measured  from 
ionograms .  because  i  t  wei ght s  the  amp l i tude  found  at  each  di f  ferrnt i al  t  tme  dc 1  ay . 
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FI  ytt  TIME -DELAY  PROFILES 


The  time-delay  profile  is  a  transform  of  the  f requency - se 1 e» ti ve- 
fading  characteristics  of  the  channel.  In  general,  the  wider  the  time- 

delay  spread,  the  more  severe  is  the  f requency - sel ec t i ve  fading.  It  should 
be  noted  that  a  single  mode  with  a  slowly  varying  time  delay  will  not  cause 
frequency-selective  fading  but  only  a  shift  in  the  arrival  time  of  the  trans 
mitted  pulse.  Slowly  varying  mean  time  shifts  are  normally  corrected  by 
resynchronization  in  baud-synchronous  communications  systems  or  by  automatic 
adjustment  in  nonsynchronous  systems. 

Discrete  peaks  in  the  time-delay  profile  are  indicative  of  the  propa¬ 
gation  of  discrete  modes  with  different  time  delays.  The  reciprocal  of  the 
time-delay  difference  between  these  peaks  corresponds  to  a  strong  frequency- 
selective  fading  component  in  the  leceived  signal.  If  the  profile  indicates 
only  two,  equal  -  strength  modes,  deep  frequency- selective  nulls  equally 
spaced  can  be  predicted. 

6.  Diurnal  Time-Delay  Spreads 

Diurnal  variations  in  the  time-delay  spread  are  shown  in  Figs.  34 
through  41.  Diurnal  trends  in  time-delay  spread  can  be  discerned;  however, 
there  is  a  wide  variation  in  the  data  from  day  to  day,  particularly  at  7  Me. 
The  high-multipath  time-delay  spread  generally  occurred  in  the  early  evening 
at  7  Me. 

7.  Tiwe-Delay-Sread  Histograms 

Histograms  for  all  observed  time-delay  spreads  calculated  at  14  and 
7  Me  for  all  data  from  3  to  14  November  are  shown  in  Fig.  42.  The  time- 
delay  spread  histogram  at  14  Me  shows  a  double  peak.  The  first  peak  indi¬ 
cates  a  predominant  single-mode  propagation  condition  with  little  time- 
delay  spreading.  The  second  peak  indicates  a  two-mode  cond  tion  with  a 
time-delay  separation  averaging  about  0.5  msec. 

At  7  Me,  time-delay  spreads  aie  more  evenly  distributed.  The  rather 
large  tail  in  th  distribution  is  worth  noting  since  it  indicates  the 
possibility  of  hihh  values  of  the  time-delay  spreads. 

8.  Doppler -Shi ft  Profiles 

A  Doppler  scattering  picture  was  obtained  by  aituully  measuring  and 
displaying  the  frequency  dispersion  of  the  channel.  The  measurement  tech¬ 
nique  has  been  described  in  Sec.  C. 2  of  Chapter  III.  Ksamplea  of 
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FIG.  34  DIURNAL  TIME  DELAY  —  6  AND  7  NOVEMBER.  7.3  Me 
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FIG.  35  DIURNAL  TIME  DELAY  —  10  AND  11  NOVEMBER,  7.3  Me 


FIG.  37  DIURNAL  TIME  DELAY  —  14  NOVEMBER,  7.3  Me 


FIG.  43  TIME  .DELAY  HISTOGRAMS 
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Doppler- frequency  profiles  are  shown  in  Figs.  43  through  46.  The  frequency 
scale  is  referenced  to  the  frequency  of  the  transmitted  CW  tone.  The  power 
density,  measured  at  1/80-cps  intervals,  is  normalized  to  the  total  ac  power 
detected.  Twice  the  second-moment  spread  measure,  2^,  is  also  shown  on  these 
profiles,  as  is  R,  the  ratio  of  the  dc  to  ac  power.  The  dc  power  is  the  on- 
frequency  power  often  referred  to  as  the  specular  component  of  the  received 
signal . 

The  Doppler- frequency  profile  is  a  tranufoi  i  of  the  t ime-sel jeti ve- 
fading  characteristics  of  the  channel.  In  general,  the  wider  the  Doppler 
spread,  the  more  rapid  is  the  time  fading.  It  should  be  noted  that  a  single 
mode  with  a  slowly  varying  mean  Doppler  frequency  will  not  cause  time  fading 
but  only  a  shift  in  the  frequency  of  the  arriving  signal.  In  communications 
systems,  slow  frequency  shifts  are  normally  corrected  by  AFC  circuits  in  the 
receiver  but,  if  not,  are  normally  small  enough  to  stay  within  the  bandwidth 
of  the  receiver  filters. 

Discrete  peaks  in  the  Doppler  profile  are  indicative  of  the  propaga¬ 
tion  of  discrete  modes  with  different  Doppler  shifts.  Th<  frequency  dif- 
terence  between  these  peaks  corresponds  to  a  strong  time-fbding  components 
in  the  received  signal.  If  the  profile  indicates  only  two,  equal-strength 
modes,  deep  time-fading  nulls  equally  spaced  can  be  predicted. 

B.  ERROR- RATF.  MEASUREMENTS 

1.  Introduction 

Before  embarking  on  a  discussion  of  the  results  obtained  in  the  error 
rate  experiment,  it  is  well  to  consider  once  again  that  the  experimental 
FSK  system  and  the  idealized  FSK  system  model  used  to  compute  the  theoret¬ 
ical  error-rate  curves  display  important  differences  in  detection  procedures. 
The  computational  model  is  un  energy-quenching  system;  on  the  other  hand, 
the  experimental  system  decays  received  signal  energy  at  a  rate  determined 
by  its  f  i  '•  *  e  r  time  constant.  In  filtering  the  received  pulses,  the  computa¬ 
tional  »-'ir.  merely  narrow -band  integrates  over  a  »  i  gnu  1  i  ng-e  )  rment  dura¬ 
tion,  while  the  exper tment a  I  system  smooths  and  shapes  the  lime  and  frequency 
characteristic*  of  the  reteived  signal  in  a  more  complex  manner.  In  anv 
critical  experimental  or  theoretical  analysis,  these  different**  will  be 
apparent  as  the  time-  and  f  requent  y  -  di  *prr  s  i  v  r  effect*  id  the  channel  ar*’ 
isolated  and  analyzed  in  suffirient  detail.  \  complete  analysis  of  the 
experimental  system  is  not  available  at  the  present  time,  and  1 1  was  not 
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FIG.  43  DOPPLER-FREQUENCY  PROFILES 
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the  purpose  of  this  project  to  generate  such  an  analysis.  The  purpose  of 
the  experiment  entailed  measuring  dispersive  effects  on  the  performance  of 
an  operational  FSK  system  and  comparing  these  effects  with  those  predicted 
by  the  existing  theoretical  model  of  an  idealized  FSK  system.  The  idealized 
theoretical  model  was  developed  in  the  early  stages  of  a  continuing  effort 
which  is  delving  into  the  complexities  of  modeling  real  channels  in  order 
to  more  accurately  predict  the  performance  of  various  communications  systems 
in  the  presence  of  arbitrary  dispersive  conditions.  This  model  serves  as  a 
tool  that  allows  the  communications  engineer  to  study  the  effects  of  a  dis¬ 
persive  propagation  medium  on  the  performance  of  a  particular  FSK  system,  a 
tool  which,  however  idealized  the  system,  yields  important  and  useful  infor¬ 
mation  about  the  class  of  communications  systems  represented  by  the  theo¬ 
retical  model. 

In  addition  to  the  problem  of  adequately  modeling  the  pecularities  of 
a  given  communication  system,  the  experimental  validation  of  an  error-rate 
analysis  is  fraught  with  difficulties  when  the  time-  and  f requency - sel ec t i ve 
nature  of  the  channel  is  modeled.  These  difficulties  are  of  a  practical 
nature  rather  than  conceptual  and  concern  the  validity  of  the  statistical 
assumptions  that  are  required  in  the  computation  of  the  error  probability. 

For  example,  it  is  assumed  that  the  channel  transfer  function  is  a  homogeneous 
Gaussian  .andom  field7  on  the  time- frequency  plane,  yet  this  fact  was  not 
verified  experimentally  (m.  ing  the  course  of  the  experiment.  Such  a  verifi¬ 
cation  was  completely  outside  the  scope  of  this  effort  and,  indeed,  in  its 
entirety  would  be  an  unavailing  task  in  the  present  framework  of  obtaining 
approximate  answers  to  an  exceedingly  complex  problem.  The  various  assump¬ 
tions  concerning  the  statistical  behavior  of  the  dispersive  mechanisms  of 
the  propagation  medium  are  based  on  physical  reasoning  and  experience  with 
radio  channels.  These  assumptions  are  general  enough  to  adequately  approxi¬ 
mate  the  behavior  of  a  physical  radio  channel  but  at  the  same  time  simple 
enough  to  yield  error-rate  expressions  that  can  be  effectively  solved  on 
a  large  digital  computer. 

Apart  from  their  use  lor  verifying  a  theoretical  model,  which  involves 
comparing  measured  and  predicted  results,  the  measured  data,  obtained  in  a 
ctrefully  planned  and  executed  experiment,  ore  very  valuable  for  assessing 
the  performance  of  an  operational  communications  system  in  a  dispersive 
environment.  These  data  can  be  used  to  indicate  problem  areas  in  the 
present  model  and  to  greatly  assist  in  the  development  of  in  improved  model. 
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In  addition  to  the  error-rate  data,  useful  measurements  of  propagation 
phenomena,  such  as  the  time-delay  and  Doppl er- shi f t  profiles  of  the  channel, 
were  obtained. 

It  is  hoped  that  this  preliminary  discussion  indicates  the  degree  of 
approximation  that  is  inherent  in  the  formulation  and  solution  of  an  en¬ 
gineering  problem  that  involves  developing  a  statistical  characterization 
of  the  dispersive  effects  of  a  natural  radio  channel  on  a  communication 
system. 

2.  Data  Classification 

The  basic  data  obtained  in  the  experiment  included  estimates  of  the 
binary  error  probability  and  S/N  over  a  30-second  interval.  The  data  were 
supplemented  with  measurements  of  the  channel  time-delay  and  Doppler  pro¬ 
files,  which  were  obtained  often  enough  to  follow  the  important  variations 
of  th»»se  functions.  Hence  any  binary  error  probability  measured  during  a 
30-second  interval  was  associated  with  an  S/N,  a  time-delay  profile,  and 
a  Doppler  profile.  This  experimental  procedure  made  it  possible  to  classify 
the  error-rate  data  according  to  time-delay  and  Doppler  spreads.  This 
classification  yields  an  experimental  determination  of  error  rate  as  a 
function  of  S/N  in  which  time-  and  f requency - sel ec t i ve  fading  effects  can 
be  isolated  and  observed  in  various  degrees  and  combinations.  In  addition 
to  being  subjected  to  this  classification  procedure,  the  30-second  esti¬ 
mates  of  binary  error  probability  and  S/Ns  were  also  appropriately  grouped 
and  averaged  so  that  a  minimum  of  100  observed  errors  was  obtained  in  the 
final  estimate  of  the  error  rate  at  any  given  S/N.  This  grouping  proce¬ 
dure  insured  the  statistical  reliability  of  the  resulting  estimate  of  the 
binary  error  probability;  many  of  the  final  error-rate  estimates  at  a  given 
S/N  are  based  on  more  than  100,000  transmitted  bits.  The  classification 
and  grouping  procedures  are  explained  in  more  detail  in  Chapter  IV.  and  the 
statistical  reliability  of  error-rate  measurements  is  considered  in 
Appendix  D.  The  error-rate  data  were  subjected  to  tour  different  classi¬ 
fication  procedures;  to  facilitate  reference  to  the  resulting  data  classes, 
it  is  convenient  to  refer  to  them  with  data-elass  designators  as  defined 
bv  Table  M. 
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1)A'.  ■  CLASS 
PKS  ‘  INATOR 

DESCRIPTION  OK  DATA  CLASS 

\ 

Tlie  class  composed  of  ail  data  in  which  the  measured  error  rate 
is  a  function  of  S/N  only. 

TDn 

Tlie  nth  of  nine  classes  obtained  by  classifying  the  data  accord¬ 
ing  to  both  time-delay  and  Doppler  spread.  Each  class  yields 
measured  error  rate  as  a  function  of  S/N  for  a  particular  range 
of  time-delay  and  Doppler  spread. 

Tn 

The  nth  of  three  classes  obtained  by  classifying  the  data  ac¬ 
cording  to  time-delay  spread  only.  Each  class  yields  measured 
error  rate  as  a  function  of  S/N  for  a  particular  range  of  time- 
delay  spread  independent  of  the  associated  Doppler  spread. 

Dn 

The  ntli  of  three  classes  obtained  by  classifying  the  data  accord¬ 
ing  to  Doppler  spread  only.  Each  class  yields  measured  error 
rate  as  a  function  of  S/N  for  a  particular  range  of  Doppler 
spread  independent  of  the  associated  time-delay  spread. 

The  range  of  the  seeond-monien  t  spread  measures  associated  with  each 
of  these  data  classes  appears  in  Table  III. 


Table  I  1 1 

TIME- DELAY  -  AND  DOPPLER- SPREAD  RANGK 
ASSOCIATED  WITH  EACH  DATA  CLASS 


DATA  CLASS 

— 

IT MC- DELAY  SPREAD 
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(msec) 

1 
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(2cta> 

( cp  s  ) 

A 

0  -  1.5 

0  -  1 

Till 

0  -  0.25 

0  -  0.25 

1D2 

0  -  0.25 

0.25  -  0.5 

TD3 

G  -  0  25 

0.5  -  1 

104 

0.25  -  0.75 

0  -  0.25 

11)3 

0.25  -  0  75 

0.25  -  0.5 

TDn 

0.25  -  0.75 

0.5  -  1 

nr 
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o  -  0.2  ‘ 

11  >H 
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measures  would  have  been  more  suitable  for  the  expferiment;  however,  the 
variations  of  the  observed  spread  values  were  sufficient  to  cause  a  de¬ 
tectable  sensitivity  of  system  performance  to  time-  and  f requency- sel ecti ve 
fading  phenomena. 

3.  Error-Bate  Data 

The  results  of  classifying  the  error-rate  data  according  to  the  four 
preceding  classification  procedures  are  presented  and  discussed  in  this 
section.  Measured  error  rate  is  plotted  as  an  X  on  the  corresponding 
theoretical  error-rate  curve.  The  theoretical  error-rate  curves  are  com¬ 
puted  for  a  transition  time  equal  to  0.02  of  the  signaling-element  duration, 
which  would  be  200  /xsec  for  a  10-msec  signaling  element.  This  choice  of 
transition  time  was  based  on  a  rough  measurement  performed  on  the  actual 
transmitting  equipment  with  a  special  wide-band  discriminator.  All  mea¬ 
sured  data  are  for  the  no-diversity  case  with  a  10-msec  signaling-element 
duration.  A  limited  amount  of  dual - di vers i ty  data  Was  processed;  these 
data  are  presented  in  Chapter  VII. 

a.  DATA  CLASS  A 

All  measured  error-rate  data,  without  regard  to  Lime-delay  or 
Doppler  spread,  are  plotted  as  a  function  of  S/N  in  Fig.  47.  The  theoret¬ 
ical  error-rate  curve  of  Fig.  47  was  computed  for  two  equal -strength  patis 
with  a  time-delay  spread  of  1  msec  and  a  Doppler  spread  of  1  cps,  winich  ire 
typical  of  the  spreads  observed  on  this  channel.  The  relatively  smooth 
variation  of  measured  error  rate  with  S/N  is  indicative  of  the  degree  of 
statistical  reliability  achieved  by  the  experiment.  Figure  47  i  1  lust rjites 
the  important  and  interesting  "  bottoming-out "  characteristic  generally 
displayed  by  the  measured  error-rate  data:  the  measured  error  rate  ceas  ;s 
to  decrease  with  increasing  S/N  and  essentially  maintains  a  constant 
irreducible  value  at  high  S/Ns,  as  predicted  by  the  theoretical  model. 

The  close  agreement  between  the  measured  error  rate  of  Data 
Class  A  and  the  theoretical  error  rate  indicates  that  the  theoretical 
model  yields  a  good  approximation  to  the  average  performance  of  the  ex¬ 
perimental  system  over  the  range  of  time-delay  and  Doppler  apreads  ob¬ 
served  on  the  channel.  It  >•  to  be  stressed  that  Fig.  47  is  concerned 
with  comparing  the  measured  performance  of  the  experimental  system 
averaged  over  all  dispersive  states  of  the  channel  with  a  theoretical 
curve  based  on  the  average  state  of  the  channel.  Close  agreement  between 
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measured  and  theoretical  results  in  this  case  does  not  necessarily  imply 
that  close  agreement  will  he  obtained  when  measured  and  theoretical  error 
.  rates  are  compared  for  a  particular  dispersive  state  of  the  channel.  Al¬ 
though  a  Doppler  spread  of  1  cps  was  the  maximum  observed,  the  theoretical 
error-rate  curves  of  Fig.  47  are  essentially  determined  by  the  time-delay 
spread  and  transition  time  and  would  not  display  a  discernible  difference 
when  computed  for  a  zero  Doppler  spread. 

b.  DATA  CLASSES  TD1  THROUGH  TD 9 

The  results  of  classifying  the  data  according  to  both  time-delay 
and  Doppler  spread  are  plotted  in  Figs.  48  through  55.  It  is  interesting 
to  observe  that  no  data  were  found  in  Data  Class  TD7 ,  which  corresponds 
to  the  combination  of  maximum  time-delay  spread  and  minimum  Doppler  spread. 
In  other  words,  during  the  course  of  the  experiment,  minimum  Doppler 
spreads  were  not  observed  when  the  time-delay  spread  was  maximum.  The 
theoretical  error-rate  curve  associated  with  each  data  class  was  computed 
for  two  equal -strength  paths  possessing  a  time  delay  spread  and  a  Doppler 
spread  equal  to  the  maximum  time-delay  spread  and  Doppler  spread  of  the 
data  class  plotted  on  the  curve. 

Observe  that  the  theoretical  error-rate  curves  associated  with 
Data  Classes  TDl  through  TD6  (minimum  and  median  observed  time-delay 
spreads)  are  essentially  identical  and  those  associated  with  Data  Classes 
TD8  and  TD9  (maximum  observed  time-delay  spreads)  are  essentially  identical. 
Hence,  the  theoretical  model  is  relatively  insensitive  to  variations  of 
time-delay  and  Doppler  spread  in  the  range  of  spreads  observed  on  the 
channel . 

The  measured  error  rate  is  these  classes  illustrates  the 
hot tomi ng-out  phenomenon;  what  is  more  important,  this  phenomenon  can  now 
be  observed  as  a  function  of  tim>’-delay  and  Doppler  spread.  The  most 
striking  and  important  characteristic  of  the  measured  error- rate  da t a  i s  the 
evident  sensitivity  of  the  experimental  system  to  time-  and  frequency- 
selective  fading  phenomena,  notwithstanding  the  predicted  insensitivity 
of  the  theoretical  model  within  the  range  of  spreads  observed  on  the 
eti  anne  l  , 

There  is  reasonable  agreement  between  the  measured  data  and  the 
theoretical  error  rate  curves  for  the  median  time-delay  classes  (TD4. 
and  TD6 »  The  error-rate  data  associated  *  1 1  h  data  classes  TDl. 


BINRRY  ERROR  PRQBRBIL I TY 


FIG.  50  DATA  CLASS  TD3  -  0  2-.  0.25,  0.5  2-.>  1-0 
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BINARY  ERROR  PROBABILITY 


FIG.  52  DATA  CLASS  TD5  —  0.25  <  2<>r  <  0.75;  0.25  <  2<'A  <  0.5 
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BINRRY  ERROR  PROBABILITY 


FIG.  53  DATA  CLASS  TD6  -  0.25  <  2"r  <  0.75;  0.5  <  2^A  <  1.0 


FIG.  54  DATA  CLASS  TD8  —  0.75  <  "  »r 


1.5;  0.25  <  *>x  <  0.5 
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TD8,  and  TD9  indicate  definite  departures  from  the  theoretical  model. 

Data  Class  TD1  is  the  minimum  ti me- del  ay- spread  and  minimum  Doppler-spread 
class;  TD8  is  the  maximum  time-del  ay- spread  and  median  Doppl er- spread  class; 
and  TD9  is  the  maximum  time-delay-spread  and  maximum  Doppl er- spread  class. 

It  is  interesting  to  note  that  the  measured  irreducible  error 
probability  for  Data  Class  TD1  is  lower  than  that  predicted  by  the  theoret¬ 
ical  model.  The  theoretical  model  bottoms  out  at  an  error  probability  of 
approximately  1  x  1 0 ~ 3 ,  whereas  the  measured  data  of  TD1  display  a  definite 
bottoming  out  at  an  error  probability  of  5  x  10”*.  Since  a  complete  analysis 
of  the  experimental  system  is  not  available,  we  are  not  in  a  position  to 
give  exact  technical  reasons  for  this  departure;  however,  in  view  of  the 
differences  between  the  experimental  system  and  the  theoretical  model  which 
have  been  discussed,  such  departures  are  not  surprising.  It  is  significant 
that  the  measured  irreducible  error  probability  is  lower  than  the  predicted 
irreducible  error  probability  for  the  data  class  corresponding  to  the 
minimum  observed  time-delay  and  Doppler  spreads,  which  implies  that  the  ex¬ 
perimental  system  is  more  effective  than  the  theoretical  model  in  combating 
a  small  amount  of  time-  and  f requency- sel ec t i ve  fading. 

On  the  other  hand,  the  irreducible  error  probability  indicated  by 
the  measured  data  of  Data  Class  TD9  is  approximately  4.2  x  1 0~ 3 ,  whereas 
the  theoretical  model  bottoms  out  at  approximately  1.5  x  10”3.  Hence,  the 
irreducible  error  probability  for  the  experimental  system  is  higher  than 
the  irreducible  error  probability  predicted  by  the  theoretical  model  for 
the  maximum  observed  time-delay  and  Doppler  spreads.  This  indicates  that 
the  experimental  system  is  less  effective  than  the  theoretical  model  in 
combating  the  degree  of  time-  and  f requency - sel ect i ve  fading  associated 
with  the  maximum  observed  time-delay  and  Doppler  spreads. 

The  fact  that  the  experimental  system  shapes  and  smooths  the 
time  and  frequency  characteri sti cs  of  the  received  signal  pulses  to  a 
greater  degree  than  does  the  incoherent  matched- fi 1  ter  model  is  probably 
responsible  for  the  superior  performance  of  the  experimental  system  for 
small  amounts  of  time-  and  f requency - se 1 ect i ve  fading.  Roughly  speaking, 
appropriate  shaping  of  the  time  structure  of  the  received  pulses  would  aid 
in  combating  f requency - se 1 ect i ve  (time-delay)  effects,  and  appropriate 
shaping  of  the  frequency  structure  of  the  received  pulses  would  aid  in 
combating  the  t ime- se 1 ect i ve  (Doppler)  effects. 
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To  continue  this  speculative  reasoning,  it  is  plausible  that,  as 
the  degree  of  selective  fading  on  the  channel  increases,  a  poi nt  i »  reached 
where  the  energy-quenching  property  of  the  theoretical  model  accounts  for  i  1 3 
superiority  at  the  more  severe  observed  time-  and  frequency- selecti ve  fading 
conditions.  In  other  words,  the  performance  of  the  experimental  system  is 
hindered  by  the  possib i 1 i ty of  an  unequal  bi  as  condi  tion  in  which  the  stored 
energies  in  the  mark  and  space  filters  are  different  enough  to  cause  a  de¬ 
terioration  in  detection  capability.  Observe  that  frequency- selecti ve 
fading  can  cause  one  of  the  signaling  frequencies  to  be  in  a  fade  over 
several  signaling  elements;  the  resulting  difference  in  stored  energy  in 
the  two  filters  will  intensify  the  detrimental  effects  of  frequency- sel ecti ve 
fading 

Finally,  to  complete  this  speculative  picture,  it  is  reasonable 
to  assume  that  the  two  opposing  effects  discussed  above  may  have  compensated 
for  each  other,  so  that  when  all  of  the  data  were  plotted  without  regard  to 
the  dispersive  state  of  the  channel  (Data  Class  A)  a  reasonably  good  fit 
with  the  theoretical  curve  was  obtained. 

The  measured  irreducible  error  rate  as  a  function  of  time-delay 
and  Doppler  spread  is  plotted  in  Figs.  56  and  57,  respectively.  The  three 
measured  points  of  irreducible  error  probability  as  a  function  .?f  time- 
delay  spread  were  derived  from  Data  Classes  TD2,  TD5,  and  TD8;  similarily, 
the  three  measured  points  of  irreducible  error  probability  as  a  function 
of  Doppler  spread  were  derived  from  Data  Classes  TD1 ,  TD2,  and  TD3.  The 
theoretical  irreducible  error  probability  curves  for  no  diversity  and  for 
dual  diversity  are  also  plotted  in  these  figures.  Comparison  of  the 
theoretical  and  measured  curves  of  irreducible  error  probability  emphasizes 
the  greater  sensitivity  of  the  experimental  system  to  time-  and  frequency- 
selective  fading  effects.  Observe  that  the  experimental  system  displays 
a  bottomi ng-out  effect  at  the  low  spread  values  which  is  in  agreement  with 
the  t rans i 1 1 on - 1 i me  bottoming  out  piedicted  by  the  theoretical  model 

c.  DATA  CI.ASSF*  T 1  .  T2,  AND  Tt 

The  results  of  classifying  the  error-rate  data  according  to 
minimum,  mediun,  and  muximurr  observed  time-delay  spread  (Data  Classes  Tl , 

T2.  and  T3 )  without  regard  to  Doppler  spread  are  plotted  in  Fig*.  58,  59. 
and  60  respectively. 
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FIG.  56  MEASURED  IRREDUCIBLE  ERROR  PROBABILITY  AS  A  FUNCTION 
OF  TIME-DELAY  SPREAD 
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BINARY  ERROR  PROBABILITY 
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FIG.  57  MEASURED  IRREDUCIBLE  ERROR  PROBABILITY  AS  A  FUNCTION 
OF  DOPPLER  SPREAD 
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BINARY  ERROR  PROBABILITY 
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BINARY  ERROR  PROBABILITY 
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1  lie  best  lit  o  I  measured  dat  a  with  I  lie  l  li re i  1 1  a  I  mi  vc  i  s  <d> 
tallied  for  Data  Cluss  12  (the  median  t  i  me- de  I  ay  -  ;s|>  read  r|.iss|,  dep.i  i  in  i  <•  s 
between  the  measured  and  theoretical  irreducible  error  p  i  ob.ib  i  I  i  I  y  are 
observed  in  the  low  and  high  t  i  me  -  de  1  ay  -  sp  read  classes.  The  measured  data 
display  the  same  behavior  observed  in  the  previous  r|  a.ssi  Inal  inn  proce¬ 
dure  and  emphasize  the  sensitivity  ol  tin*  per  Id  rmaiiee  ol  tin*  experimental 
system  t  f requeney- sc  1 ec t i ve  effects.  The  irreducible  error  probability 
measured  forDatu  Class  1*3  is  approximately  ten  times  that  measured  lor 
Data  Class  IT.  This  is  an  impressive  variation  of  error  rate  for  the 
relatively  small  variation  of  time-delay  spread  associated  with  the  two 
d  a  la  classes. 


d.  DATA  CLASSES  Dl.  D2.  AND  D.i 

The  results  of  classifying  the  error-rate  data  a< cording  to 
minimum,  median,  and  maximum  observed  Doppler  spread  (l)at.i  Classes  PI  , 

D2 ,  and  1)3)  without  regard  to  time-delay  spread  are  plotted  i  n  Fig.-.-  til, 

62,  and  63 

It  is  interesting  to  observe  that  the  experimental  data  display 
a  definite  detectable  sensitivity  to  Doppler  spread.  For  the  range  of 
Doppler  spread  observed  on  the  channel,  the  theoretical  model  is  not  sen¬ 
sitive  to  t i me- se 1 ec t i ve  effects,  since  the  f requeue}- - se I ec t i ve  and 
trausi tion- time  effects  are  the  major  contributing  factor-  affecting  the 
performance  of  the  therorot i ca 1  model  for  the  range  oj  observed  time-delay 
and  Doppler  spreads.  However,  the  measured  i n cducibl c  error  probabilits 
of  Data  Classes  1)1,  D2,  and  1)3,  indicates  a  small  but  definite  L  i  me- se  I  ec  t  i  \  e 
sens i t i v i ty . 

C  I'AIMTUMKNT  ItKVIEW 

The  experience  obtained  from  conducting  such  a  complex  ••xjotimenl  us 
this  naturally  suggests  improvements  that  could  lie  made  if  the  experiment 
were  ripe. ited.  tutin  suggestions  are  offered  for  possible  future  releiencc- 

ll)  Herotd  all  t  air  datu  automat  i cal  I y  and  directly  in  u 
•  •in/Mi  /  •'  t  -  com/iii  /  i  b  I  e  input  "l•,d^um.  Many  pimi’-siin 
problems  llint  caused  delnw  were  direct  I  v  irui.uhl  ■ 
to  errors  in  the  manual  rve-irdi  ng  of  dit.i,  for  ex¬ 
ample,  m  i  spun  cb  cd  data  c..i»Ij  and  on<ittci  l,.g  eni  1  i  . 

A  small  sample  of  erroneou-  ilata  can  sij.ui  I  . .  ini  I  \ 
bias  the  results  in  such  cn  expo i i men t . 
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(2)  Increase  the  data  transmission  rate.  It  was  des  i  tab  I  <• 
to  obtain  the  error  rate  for  a  number  of  classes  of 
channel  conditions.  However,  the  100-baud  maximum 
transmission  rate  an2  the  over-all  time  schedule 
limited  the  amount  of  data  that  could  be  collected 
for  each  class.  Confident  error-rate  estimations 

at  low  error  rates  require  a  large  number  of  transmitted 
bits.  More  accurate  validation  of  error-rate  character¬ 
istics  of  diversity  systems  would  require  the  measure¬ 
ment  of  much  lower  error  rates  than  attempted  here. 

With  more  data,  the  parameter  of  diversity  channel 
correlation  could  be  more  accurately  accounted  for. 

The  data  rate  for  the  AN/FGC-29  system  could  be  in¬ 
creased  simply  by  adding  f requency- di v i si  on  multi¬ 
plexed  channels.  This  would  also  permit  experiments 
to  be  conducted  on  the  effect  of  interchannel  dis¬ 
tortions. 

(3)  Transmit  other  binary  sequences .  The  alternating 
10-msec  keying  sequence  used  was  the  most  sensitive 
sequence  possible  for  observing  effects  of  time- 
delay  spread  on  error  rate.  The  error  rate  with  the 
alternating  sequence  cannot  be  more  than  four  times 
the  error  rate  with  a  random  sequence  and  is  probably 
closer  to  two  times.  However,  a  random  sequence  may 
be  more  realistic  for  evaluating  the  performance  of 

an  actual  data  system  under  various  channel  conditions. 

The  error-rate  model  uses  channel  -  spread  parameters 
that  are  normalized  to  the  transmission  rate  and  to 
the  mark-space  frequency  spacing.  Therefore  increasing 
the  element  length  affords  a  method  of  effectively 
simulating  reduced  channel  time-del  ay- spread  conditions. 
Transmission  of  a  simple  11001100  sequence  would  reduce 
the  inter-symbol  interference  effects  by  one  half.  In 
fact,  with  sufficiently  long  signaling  elements,  measure¬ 
ment  of  the  error  rate  with  a  sequence  of  all  marks  may 
greatly  reduce  frequency-selective  fading  effects  and 
may  permit  error- rate  validation  of  an  essentially  pure 
Doppl er- spread  channel. 

In  like  manner,  the  spacing  of  the  murk-space  frequencies 
could  be  widened  to  simulate  reducer  Doppl er - spread 
channel:*.  This  could  be  done  conveniently  in  the 
\N  FIX'.- 2**  system  by  signaling  with  frequencies  in  dil- 
ferent  channels. 

(4)  ('unduct  the  experiment  over  a  shorter  oath  length.  Time- 
delay  and  Doppler  effects  on  the  Fort  •noiouth  to  Palo  Mto 
path  were  minimal.  \  shorter  path  would  give  much  more 
variation  of  t i me- de I av -  spread  conditions  This  would 
permit  belter  validation  of  error- rate  models  for  higher 
time-delay  spreads. 
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(5)  Conduct  the  experiment  over  a  path  with  different  iono¬ 
spheric  his tory  .  A  path  where  auroral  effects  or  spread-F 
conditions  are  more  prevalent  would  give  more  Doppler 
variations  and  hence  a  better  verification  of  the  error- 
rate  models  under  various  Doppler- spread  conditions. 

Measurements  at  additional  transmission  frequencies  and  at  different 
times  of  day  would  also  provide  data  at  other  spread  conditions. 

Everything  considered,  it  was  felt  that  the  experiment  was  very 
sucessfui.  We  are  unaware  of  any  reported  channel  error-rate  experiment 
conducted  in  this  depth;  simultaneous  measurement  of  appropriate  channel 
parameters  and  real  system  error  rate  to  verify  assumptions  associated 
with  a  mathematical  model  for  both  the  channel  and  the  system. 
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VI  FREQUENCY -SELECTIVE  WAVEFORM  ANALYSIS 


A.  INTRODUCTION 

The  word  analyzer  used  in  the  error- rate  experiment  provided  immedi¬ 
ate  detection  of  received  data  bit  errors  With  this  capability,  i t  was 
possible  to  record,  by  oscilloscope  photography ,  received  signal  wave¬ 
forms  at  various  processing  stages  in  the  receiver  system  for  the  time 
interval  immediately  following  a  detected  error.  The  basic  technique 
utilized  a  fast  response,  dual-channel  oscilloscope  to  display  the  wave¬ 
forms  and  a  Polaroid  Land  camera  to  photograph  them.  By  triggering  one 
or  two  such  oscilloscope-cameras  with  the  error  signals,  it  was  possible 
to  record  in  detail  as  many  as  four  received  waveforms  simultaneously  in 
the  time  interval  immediately  following  the  occurrence  of  an  error.  Dur¬ 
ing  that  interval,  several  errors  usually  occurred  because  of  the  burst 
character  of  the  errors.  Some  of  these  records  and  t.hei  r  interpretation 
are  presented  in  this  chapter. 

B.  DESCRIPTION  OF  EXPERIMENT 

Detrimental  frequency -selective  effects  for  an  FSk  system  would  be 
evidenced  by  unequal  and  time-varying  channel  attenuations  at  the  mark 
and  space  signaling  frequencies.  To  observe  whether  this  condition  oc¬ 
curred,  it  was  necessary  to  compare  the  amplitude  of  the  mark  and  -ipace 
frequencies  as  a  function  of  time.  The  AN  FtiC-29  transmitter  sent  FSK 
tones  separated  by  85  cps.  Alternate  mark  and  spare  frequencies  were 
transmitted  in  10-msec  intervals.  Because  of  certain  d i ve rs i t v - combi n i n g 
techniques  used  in  the  AN  EGO -  29  receiver,  it  was  not  convenient  to  ob¬ 
tain  u  true  picture  of  c  hannel  frequency- select  ive  fading  effects  iti  the 
normal  signal  processing  by  this  receiver  Hence  a  special  mark-space 
filter  was  constructed,  des.gned  to  provide  separation  of  the  two  fre¬ 
quencies.  This  filter  was  connected  to  ihe  output  of  the  delay  amplifier 
of  the  AN  K(jC-  29  receiver. 

A  block  diagram  of  this  mark-space  filtet  is  shown  in  Iig  t-1-  The 
received  signal  passed  bv  the  mark  filter  was  envelope  detect'd  with  a 
positive  polaritv.  The  received  signal  passed  Sv  the  spare  filter  was 
detected  with  a  negative  polaritv  These  two  r  tputs  sere  combined  across 
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FIG.  64  BLOCK  DIAGRAM  Of  MARK-SPACE  FILTER 

a  grounded  common  load  resistor.  The  voltage  across  this  resistor  then 
increased  positively  when  a  mark  frequency  was  received,  in  proportion 
to  the  strength  of  the  received  signal,  and  similarly  increased  nega¬ 
tively  when  a  space  frequency  was  received.  The  3-db  band\vidfl\  of  each 
filter  was  about  40  cps. 

The  output  of  the  mark- space  filter  was  connected  to  one  input  of  a 
dual-channel  oscilloscope,  as  shown  in  Fig.  65.  Error  pulses  generated 
by  the  wo;d  analyzer  triggered  the  oscilloscope.  The  triggering  mode 
was  set  to  “ s i ngl e - sweep ”  mode  operation.  In  this  mode,  after  the  trig¬ 
gering  circuit  was  manually  “armed,”  the  first  external  trigger  pulse  to 
arrive  triggered  a  single  scope  sweep.  Additional  trigger  pulses  had  no 
effect  on  the  scope  until  the  trigger  circuit  was  again  armed. 

Transient  waveforms  were  recorded  during  detected  errors  b\  the  fol¬ 
lowing  procedure: 

(1)  Prepare  the  Polaroid  camera  on  the  oscilloscope. 

(2)  Open  the  camera  shutter. 

(3)  Arm  the  *•  ope  trigger  circuit. 

(4)  14a  it  for  a  del  et  ted  error  to  occur 

( 5  >  Close  shutter  and  process  the  photograph 

Alternating  mark  and  space  signaling  element s  of  10  msec  duration 
were  transmitted  To  observe  the  distortion  details  of  each  pulse,  the 
scope  sperd  was  adjusted  to  2'!  or  50  msec  cm.  I  reqti  no  \  -  s  e  I  ec  t  i  v  r  fad¬ 
ing  »  h  a  r  «*c  t  e  r  i  s  t  I  c  s  on  the  channel  srf"  observed  w  i  t  h  a  slower  sweep  speed. 
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The  recordings  were  made  durng  operation  at  7.3  or  14.3  Me.  Hie 
14- Me  recordings  were  made  between  the  hours  of  1800  anil  2300  GMT  and  the 
7.3-Me  recordings  between  0000  and  0200  GMT.  The  exact  date  and  time  of 
each  recording  is  indicated  on  the  figures.  The  average  S/N  when  these 
recordings  were  made  was  generally  higher  than  20  db 


C.  WAVEFORMS  WITH  FREQUENCY- SELECTI VE  FADING 

Figure  66  shows  an  example  of  almoct  identical  fading  properties  on 
the  mark  and  space  frequencies.  For  a  channel  of  this  type,  error  prob¬ 
abilities  are  predicted  by  assuming  a  channel  model  in  which  the  signal¬ 
ing  frequencies  fade  . imu l taneousl y :  produce  fading  that  is  flat  or 
completely  correlated  at  the  two  signaling  frequencies. 


FIG.  66  FREQUENCY-FLAT  FADING 


Figure  67  shows  two  examples  of  noticeable  frequency  selectivity  be¬ 
tween  the  mark  and  space  frequencies.  The  time  positions  of  actual  de¬ 
tected  errors  are  shown  immediately  below  this  waveform  as  a  series  of 
dots.  The  occurrence  of  error  bursts  was  quite  common  during  frequenrv- 
selective  fading  of  this  sort.  Furthermore,  ns  can  he  seen  from  these 
examples,  bit  errors  tended  to  occur  alternately.  lots  implies,  because 
of  the  alternating  sequence  transmitted,  that  the  I'lnos  wt*r**  of  one  tvpe 
trither  a  mark  Manspnsed  to  a  space  or  a  space  transposed  t  «•  a  mark). 

Figuie  bH  shows  two  examples  of  frequec*-  flit  fading  on  a  longer 
time  scale.  Figure  h9  show*  two  examples  *•  ’  f requenr S  se ■ ect i ve  fading 
in,  the  same  time  scale  as  that  of  Fig  Hoi*  it  *  an  b*  -cm,  that  the 

is  t  iti  ii  I  t  an  onus  channel  gams  tit  the  mark  and  'pace  lii.miu  U  »*rr  «»u  i  t  r  dt  .* 
ferent  There  are  at  least  two  occasions  on  these  fonnb  where  the  gam 
a*  one  frequency  was  in  a  null,  while  at  a  frequen«v  <>nlv  8  cps  awavthe 
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FIG.  6?  FREQUENCY  -SELECTIVE  FADING 
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FIG.  69  FREQUENCY-SELECTIVE  FADING  ON  SPACE-DIVERSITY  ANTENNAS 


the  channel  gain  was  at  a  relative  maximum.  These  t requency- se 1 ec t i ve 
effects  can  be  attributed  to  the  differential  time  delays  asroci  ated  wi  tli 
multipath  propagation  modes. 

D  rut. FORMS  ON  SPACE  ANTENNAS 

Hi  e  upper  waveform  in  tigs-  68  and  t>6  was  obtained  s  itnti  1 1  an  eou  s  1  v 
with  the  lower  waveform  from  a  second  AN  tC  -,t»  receiver  and  mark- spare 
filter  connected  to  the  second  chann.'l  <>l  the  dua !  -  <  ti  ann  e!  sc<>pe.  flic 
antenna  for  this  receiver,  however,  was  spaced  some  I . 000  feet  awav  from 
the  first  antenna  At  these  tunes,  the  waveforms  seem  to  indicate  that 
the  time  fadmg  on  these  two  antennas  was  uiicm  t  el  at  ed,  because  no  deep 
fading  incurred  s  i  mu  I  t  an  etui  s  l  v  on  both  antennas.  Sjowevet.  it  is 
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significant  that  the  frequency- sel ec ti ve  characteristics  on  the  two  an¬ 
tennas  were  similar. 

The  straight  line  at  the  beginning  of  the  sweep  in  the  waveforms  of 
Fig.  69  occurred  because  the  sweep  was  manually  triggered  while  the  tones 
were  off  at  the  transmitter.  This  occurred  every  30  seconds  in  the  trans¬ 
mission  schedule;  it  provides  a  convenient  zero  reference  level  in  these 
reco  rds . 

E.  LINE- LOOP  CURRENT  WAVEFORMS 

In  this  experiment  it  was  also  possible  to  observe  the  waveforms  of 
the  dc  line- loop  keying  signal.  This  signal,  available  as  the  output  of 
the  AN/FGC-29  combiner  unit,  is  the  limited  v  rsion  of  the  discriminator 
output  normally  delivered  to  the  user.  Distortion  on  this  signal  is  often 
measured  and  evaluated  to  determine  the  usefulness  of  the  channel.  This 
signal  is  the  input  to  the  word  analyzer  in  this  experiment.  Figure  70 
shows  two  examples  of  this  keying  waveform.  An  error  detected  by  the  word 
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FIG.  70  LINE-LOOP  CURRENT  WAVEFORMS 
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analyzer  using  this  signal  is  indicated  on  the  figure  l>y  a  dot  immediately 
below  the  waveform. 

If  a  mark  state  j  arbitrarily  defined  as  the  upper  level  on  this 
waveform,  then  the  initial  portion  of  the  top  waveform  shows  what  is  known 
as  spacing-bias  distortion.  Spacing  bias  tends  to  cause  transmitted  marks 
to  be  erroneously  interpreted  as  spaces.  Notice,  however,  that  the  dis¬ 
tortion  soon  reverses  to  marking  bias,  and  errors  of  the  opposite  type  oc¬ 
cur.  Such  behavior  can  be  explained  by  f requency * se l ec t i ve  fading  on  the 
mark  and  space  channels. 


F.  SIMULTANEOUS  WAVEFORMS  DURING  ERRORS 

Through  the  use  of  two  du  al  *  chann  el  oscilloscopes  and  Polaroid  cam¬ 
eras,  triggered  by  the  same  error  pulse,  it  was  possible  to  record  simul¬ 
taneously  four  waveforms.  figures  71  through  76  show  examples  of  these 
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FIG.  71  SIMULTANEOUS  WAVEFORMS  DURING  ERRORS 
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FIG.  72  SIMULTANEOUS  WAVEFORMS  DURING  ERRORS 
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FIG.  73  SIMULTANEOUS  WAVEFORMS  DURING  ERRORS 
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FIG.  74  SIMULTANEOUS  WAVEFORMS  DURING  ERRORS 
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FIG.  75  SIMULTANEOUS  WAVEFORMS  DURING  ERRORS 
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FIG.  76  SIMULTANEOUS  WAVEFORMS  DURING  ERRORS 
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records.  The  horizontal  time  scale  for  Fig?.  71  through  75  is  50  msec 
per  grid  line.  In  Fig.  76  the  scale  is  100  msec  per  grid  line.  The 
waveforms  from  top  to  bottom  are  as  follows: 

(1)  Mark-space  filter  output  from  space-diversity  receiver 

(2)  Mark-space  filter  output  from  normal  receiver 

(3)  Line-loop  keying  signal  from  normal  receiver 

(4)  Dots  indicating  signaling  elements  detected  in 
error  from  normal  receiver  only. 

The  gains  of  the  two  mark-space  waveforms  were  equalized  with  a  com¬ 
mon  signal  into  the  antenna  terminal.  An  attempt  was  made  to  center  the 
mark-space  filter  output  displays  vertically  so  that  zero  output  was  on 
a  hori.  ntal  grid  line.  The  line-loop  keying  signal  waveform  is  reversed 
in  polarity  with  respect  to  its  corresponding  mark-space  signal  shown 
above  it.  Hence  f requency - se 1 ec ti ve  fading  in  the  space  frequency  (upper 
level)  will  be  evident  by  the  trace  remaining  below  the  zero  threshold. 

At  the  same  time,  the  line- loop  keying  waveform  will  be  distorted  with 
marking  bias. 

7/  , 

Figure  ?4  is  a.i  interesting  example.  At  the  start  oi  the  sweep,  the 
channel  was  evidently  in  a  selective  fade  on  the  mark  channel  but  not  in 
the  space  channel.  Asa  result,  there  was  heavy  spacing  bias  as  indicated 
on  the  loop  keying  waveform,  and  errors  were  made  on  every  reception 
of  a  mark  signaling  element.  Then  there  followed  a  period  where  the 
fading  was  apparently  more  or  less  equal  on  the  mark  and  space  channels 
and  the  keying  waveform  is  not  noticeably  distorted.  However,  error.1,  were 
detected  for  tight  signaling  elements  in  a  rov  The  channel  then  went  in¬ 
to  a  condition  of  selective  fading  on  the  space  channel,  and  enors  were 
made  only  in  the  detection  of  space  signaling  elements  ()»■  possible  ex¬ 
planation  of  this  plienor.cn  is  to  Hypothesise  the  simultaneous  existence 
of  the  three  propagation  modes,  two  of  which  were  rather  strong  and  equal 
in  strength  but  with  small  differential  time  delay.  The  third  mode  would 
have  had  to  be  weaker  with  a  relatively  lung  delay.  (Hiring  the  time  that 
the  two  stronger  modes  arrived  a<  the  antenna  in  phase  opposition,  the 
weaker  mode  dominated.  However,  because  svnehrom rat  ion  was  stiil  deter¬ 
mined  bv  the  arrival  time  of  the  two  stronger  paths,  sampling  *»l  the  sig¬ 
naling  elements  with  ret « pt ion  from  the  longer  path  was  *«fiM  irrtt  h  dis¬ 
placed  to  cause  errors  in  every  eieaent- 
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G.  WAVEFORMS  WITH  DIVERSITY  OPERATION 


An  attempt  was  made  to  observe  waveforms  after  a  detected  error  dur¬ 
ing  operation  with  spaced  diversity.  However,  diversity  operation  errors 
were  so  rare  that  the  waiting  time  to  obtain  records  following  the  detec¬ 
tion  of  an  error,  was  prohibitive.  Observation  of  the  diversity  wave¬ 
forms  shown  in  the  previous  figures  indicates,  however  that  when  errors 
were  detected  on  the  normal  receiver,  very  few  errors  would  have  been 
made  had  the  system  been  operating  with  diversity.  At  the  time  these  wave¬ 
form  records  were  made,  antenna  correlation  measurements  were  al 1  under  0 . 6 . 


H.  CONCLUSIONS 


Hi gh- time- resolution  waveform  recordings  for  a 
ately  following  the  detection  of  an  error  were  made 
tions  system  over  an  HF  path.  These  have  shown  the 
frequency- selecti ve  effects  on  the  performance  of  a 
The  observations  of  these  records  can  be  summarized 


time  interval  immedi- 
on  an  FSK  communica- 
si  gni  f  i  cani  e  of 
no-diversity  system, 
as  follows: 


(1)  Even  with  the  rather  narrow  frequency  sep; ration  of 
85  cps  used  in  the  AN/FGC-29  system,  the  assumption 
that  the  two  signaling  frequencies  always  fade  to¬ 
gether  is  not  valid. 

(2)  Man-  of  the  errors  observed  occurred  in  bursts  and 
were  all  of  one  type:  either  a  mark  transposed  to 
a  space  or  a  space  transposed  to  a  mark.  In  many 
cases,  the  errors  at  the  beginning  of  the  burst  were 
of  one  type,  and  a  transition  to  errors  of  the  op¬ 
posite  type  occurred  during  the  burst.  No  cases  were 
observed  where  the  type  of  errors  during  a  burst  was 
random.  Thus  the  errors  during  a  burst  cannot  be  ex¬ 
plained  simply  by  a  low  S/N. 

(3)  Bias  distortion  of  the  line-loop  keying  waveform  is 
correlated  with  observed  f requency - se 1 ec t i ve  effects 
on  th e  jh anne 1 . 


The  number  of  records  made  was  insufficient  for  the  formation  of  con¬ 
fident  statistics  as  to  the  length  of  error  bursts  or  the  type  of  error. 
Several  cases  were  recorded  where  only  an  isolated  error  occurred.  How¬ 
ever,  bursts  of  five  to  fifteen  errors  were  also  quite  common. 
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VII  SPATIAL  CORRELATION  MEASUREMENTS 


A.  INTRODUCTION 

l .  Di versitv  Systems 

Radio  communications  systems  have  long  utilized  diversity  techniques 
to  combat  the  effects  of  signal  fading  and  thereby  improve  performance. 

In  strategic  HF  military  links,  the  use  of  dual  spaced- antenna  diversity 
is  almost  standard.  The  key  to  the  success  of  a  diversity  technique  is 
the  availability  of  uncorrelated  signals  at  the  detector  of  the  receiver. 

A  comparison  of  dual-  and  no-diversity  systems  with  respect  to  error  rates 
as  a  function  of  S/N  shows  th'  dual  diversity  consi derablv  improves  per¬ 
formance.  Futhermore,  the  margin  of  improvement  increases  with  increasing 
S/N  up  te  the  point  where  dispersive  effects  predominate.  Previous  work 
by  Staras,25  Packard,22  and  Pierce  avid  Stein23  has  attempted  to  assess  the 
effect  on  system  performance  of  partially  correlated  diversity  signals. 

The  results  have  indicated  that  for  signal  correlations  up  to  0.6  there  is 
negligible  deterioration  in  performance  over  the  zero-correlation  case. 
Little  analysis  har  been  made  fcr  time-  and  f requency - sel ec t i ve  fading  of 
the  channel  with  correlated  inputs. 


2.  Correlation  Coefficient 

Assume  that  before  detection  each reeei ved  signal  of  a  dua 1  - di vers i ty 
communications  system  can  be  represented  by  sample  functions  from  two  com¬ 
plex  random  processes,  x{t)  and  y(t).  A  convenient  measure  of  the  efter 
tiveness  of  diversity  combining  is  the  normalized  covariance  or  correla¬ 
tion  coefficient,  p  ,  of  these  functions  at  zero  time  shift, 

r  x  y 


P,, 


[(*  -  V )  ( y  -  v )  ] 


(1) 


where  x  and  y  are  the  respective  means,  and  o-2  and  a2  are  the  respective 
variances  of  x(t)  and  y(t).  The  value  of  p  is  generally  constrained  to 
the  range  from  ~1  to  +1.  A  p  of  +1  or  ~1  indicates  that  with  prob¬ 
ability  one,  a  deterministic  linear  relationship  exists  between  ,t(f)  and 
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y(t),  while  a  /o  of  zero  indicates  they  arc  unco  r  re  I  a  ted  and,  if  x(  t.  )  and 
y{t)  are  jointly  Gaussian,  independent. 

It  is  interesting  that  the  correlation  coefficient  between  the  squared 
envelopes  of  Gaussian  processes,  x(t)  and  y(t),  is  simply  the  square  of  /- 

*  y 

This  would  indicate  that  the  aquared  envelope  correlation  coefficient  should 
always  be  greater  than  zero.  It  is  also  true  that  the  envelope  correlation 
coefficient  for  the  jointly  Gaussian  process  is  greater  than  zero.  De¬ 
partures  from  this  would  indicate  that  the  sample  functions  arc  not  derived 
from  Gaussian  i andom  processes. 

3.  Purpose  of  the  Measurement 

The  original  objective  of  the  measurement  of  correlation  coefficient 
was  to  provide  information  that  would  permit  comparison  of  a  theoretical 
error-rate  model  for  dual  diversity  with  experimental  error-rate  data. 
Unfortunately,  it  was  not  possible  to  collect  enough  dual -di versi ty  error- 
rate  data  for  definite  conclusions.  Measurements  of  correlation  coeffi¬ 
cient  were  made,  however,  and  at  selected  times  the  data  receiver  was 
operated  in  dual  diversity.  From  the  data  tak  n,  certain  interesting  facts 
emerged:  (.1)  when  the  dual -di  versi  ty  system  made  errors,  the  envelope 

correlation  coefficient  was  high;  (2)  a  relationship  existed  between  long¬ 
term  trends  in  envelope  correlation  coefficient  and  time-delay  spread;  and 
(3)  previously  held  rules  for  space- di vers i ty  antenna  spacing  may  not  be 
correct.  In  this  chapter,  physical  models  for  space- selecti ve  phenomena 
are  given;  the  measurement  techniques  are  described;  and  the  results  of 
the  measurement  are  d: scuTSt d. 

B.  PHYSICAL  NATURE  OF  SPACE  -  SELECTIVE  FADING 

At  HF,  a  common  technique  for  obtaining  uncorrelated  received  signals 
is  to  use  signals  received  on  antennas  physically  spaced  by  several  wave¬ 
lengths.  Recent  work  by  Ames,  however,  has  shown  that  this  simple  rule 
may  not  be  adequate.  24  lie  ha.,  experimentally  demonstrated  the  presence  ol 
quasi-periodic  interference  patterns  along  the  ground,  which  arise  from 
differences  in  arrival  ang'e  between  ionospheric  modes.  On  this  basis, 
correlation  of  signals  received  cn  spaced  antenna  will  be  highest  when 
the  center  line  between  the  two  antennas  is  perpendicular  to  the  inter¬ 
ference  patterns  and  lowest  when  it  is  parrullel  * o  the  interference 
patterns . 
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ANTENNAS  ORIENTED 


FIG.  77  TYPICAL  MODE  INTERFERENCE  PATTERNS 


1.  Mode  Interference  Patterns 

Shown  in  Fig.  77  are  two  modes  differing  in  elevation  angle  but 
identical  in  azimuth.  This  could  represent,  for  example,  high-  and  low 
ray  F-propagating  modes,  each  with  a  different  time  delay,  Doppler  fre¬ 
quency,  and  incoming  angle  of  arrival.  If  the  transmitted  signal  is 
represented  as 

y{t)  ~  Re{e  ^°  }  ,  (2) 

the  received  signal,  Z(x,y,z,t)  ^ Z(r,t ),  is  thus 


Z(r,t) 


Re  IV 


i  [  2 j  (  f  +  r  J  )  +  k  J  '  r] 


i  [  2tt\  .  (  t  +r  „  )  +  T  „  *  r]l  t  2rr  f  t] 


+  42e  ‘  2'  2  2 


where 

k.  =  ith  mode  Doppler  frequency 
r.  =  ith  mode  time  delay 

k  =  mode  vector  Wave  number,  =  k  i  +  k  7  +  fe  £ 

r  =  xt  +  yj  +  zTi  =  general  position  vector 
A  .  =  complex  amplitude  of  ith  mode. 


(3) 
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For  illustration,  suppose 


i  2  w  A  ,  r .  i  2  A  ^  r  „  +  S 

V  1  1  -  V  2  =  -4 


(4) 


Then  Eq.  (3)  becomes 


Z( 


r,0  , 


(5) 


Z(r,t)  will  be  zero  when 


S  +  27r(\2~  A.j  )  t  +  (£2  -  tj)  •  r  =  (2n  -  1  >  n  =  0.+1.+2  ... 


(6) 


and  will  be  a  maximum  when 


S  +  2tt{K2  -  +  (t2  --  X.)  •  7  =  2n7r  n  =  0 ,  +1 ,  =  2,  ... 

(7) 

The  effect  of  time  delay  is  merely  to  put  a  fixed  phase  shift  in  the 
system  of  traveling  waves.  If  the  two  modes  have  unequal  amplitudes, 

Eq.  (7)  still  describes  the  occurrence  cf  maxima,  but  Eq.  (6)  describes 
the  occurrence  of  minima  rather  than  zeros. 

At  a  given  point  in  space  [(X2  -  Xj )  *  r  =  a  constant],  a  zero  occurs 
every  1/(A2  ~A.j)  seconds.  If  one  stood  at  a  given  point  and  measured  field 
intensity  on  a  narrow-band  receiver,  the  time  between  minima  (or  maxima) 
would  depend  only  on  the  Dopp 1 e r- f requency  difference  between  the  two  modes. 

Suppose  that  instead  of  standing  at  a  given  point  we  position  many 
antennas  along  a  straight  line.  For  simplicity  (with  reference  to  Fig.  77), 
we  make  this  line  have  the  same  azimuth  as  our  two  rays.  From  Eq.  (7)  we 
see  that  the  locus  of  maxima  is  satisfied  by  the  equation  of  the  plane 

(X2  “  X,)  *  7  =  2nrr  -  2 w(\2  -  \j)f  ,  (8) 

or,  by  expanding, 

( kx2  ~  k,2^x  +  ^ky2  "  ky\)y  +  ~  kz}')Z  =  2nJT  “  277(/V2  K\  >  *  • 

(9) 
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By  the  choice  of  the  coordinate  system,  the  two  waves  have  no  contribution 
in  the  x  direction,  kx2  =  =  0,  and  furthermore  we  are  observing  on  the 

plane  2=0.  Under  these  conditions,  only  the  middle  term  on  the  left-hand 
side  of  Eq.  (S)  survives,  and  by  solving  for  y 


2tt [n  -  (\2  -  \j)  t) 

y  =  - - - - -  n  =  0,  +1 ,  +2,  ...  (10) 

k  _  -  k  . 

y  2  y  1 


The  distance  between  successive  maxima  (or  minima)  will  be  27 7/(^>2  -  &  j)- 
This  is  the  period  of  the  interference  pattern  shown  in  Fig.  77.  This 
pattern  moves  with  time  at  a  velocity  of  2~ 7(A.,  +  2  ~  j  )  along  the 

y  axis. 


Certain  physical 
a  scalar  wave  funtion 


constraints  are  imposed  by  virtue  of  the  fact  that 


Z( r , t )  =  Re^e 

must  satisfy  the  wave  equation 


[  2  7T  ( A  +/ft  )  t  +  * 


i 


V*z 


3  2Z 

fie  - 

dt2 


(11) 


(12) 


which  leads  to  the  relationship 

(27 r)2fie(K  +  /0)2  =  |X|2  =  k2  +  k2  +  k\  .  (13) 

As  might  be  expected,  the  wave  numbers  in  the  three  principal  directions 
are  constrained  by  the  physics  of  the  medium. 

2.  Numerical  Examples 

For  our  first  example,  let  us  assume  that  our  two  waves  are  a  2 F 
high  and  low  ray  on  the  4100-km  Fort  Monmouth  to  Palo  Alto  path.  On  the 
basis  of  vi  rt.ual -height  geometry,  we  can  assume  t?  and  6^  to  be  11  and 
15  degrees,  respectively.  If  we  assume  a  frequency  of  14  Me  and  that  A ^ 
and  \2  <<  f0<  we  g-  i  ,  sn.ee  k%  -  0, 
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k2  +  k2 

y  t  « 


(2tt) 


1 


3  x  108> 


(14  x  io6) 


(14) 


and 


«  Tan'1  15°  ; 


■y  2 


Solving,  we  get 


2  1 


Tan'1  11° 


(15 


y  1 


y  2 


0.267  (meter)'1 
0 . 262  (meter)'1 


(16) 


The  distance  between  maxima  is  then 


2  7T 


k  .  ~  k  „ 
y  i  y  2 


2v 

0.^05 


1256  meters 


(17) 


for  a  frequency  of  14  Me.  (At  7  Me,  kfl  and  k^2  become  much  closer,  and 
the  distance  between  maxima  is  greatly  increased.) 

Suppose  now  that  our  two  incoming  rays  are  2 F  and  3 F  low  rays  having 
arrival  angles  of  11°  and  20°,  respectively.  Under  these  conditions,  at 
14  Me, 


and 


0.267 

0.236 


(18) 


202  meters 


The  distance  between  minima  (or  period  of  the  interference  pattern) 
for  2 F  and  3 F  is  only  about  one-sixth  that  of  the  period  of  the  pattern 
produced  by  high  and  low  2 F-  This  drastic  change  shows  the  difficulty 
of  using  a  single  antenna  spacing  for  different  propagation  conditions. 
In  practice,  spacings  on  the  order  of  ten  wavelengths,  chosen  to  satisfy 
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constraints  of  land  use,  should  probably  be  adequate,  provided  that  the 
line  between  antennas  is  approximately  in  the  azimuthal  di icction  of  the 
incoming  waves  (not  necessarily  the  great-circle  bearing). 


3.  INTFRFERENCF.  PATTERNS  AS  RANDOM  PROCESSES 

The  preceding  two  sections  treat  a  somewhat  special  case  of  only  two 
modes  at  constant  amplitude.  Generally,  it  is  necessary  to  consider  many 
modes  having  amplitudes  that  are  random  processes.  In  order  to  do  this, 
an  approach  similar  to  that  used  by  Daly  in  describing  time-  and  frequency- 
selective  fading  will  be  employed.  6 

Assume  that  the  transmitted  signal  is  continuous  of  frequency  / 

The  channel  is  assumed  to  be  varying  in  time  and  giving  rise  in  space  to 
a  number  of  discrete  plane  waves  having  different  Doppler  shifts  and 
directions.  Since  we  are  dealing  with  a  band  of  frequencies  only  a  few 
cycles  at  most  removed  from  /Q,  we  neglect  frequency- sel ecti ve  effects. 
Under  these  conditions,  the  output  of  the  channel,  analogous  to  Eq.  (3) 
ill  the  last  section,  is  for  n  plane  waves 


Z(r,t )  =  Ref  J  ^  a  .  e 


i  [  2  •  t+k  *  7  T]  i  2  v  f 


The  complex  envelope  of  Z(r,t)  denoted 


It  is  assumed  that  the  a(  are  zero-mean  complex  Gaussian  random  variables 
having  independent  real  and  imaginary  parts.  This  last  constraint  imposes 
uniformity  upon  the  phase  of  the  a..  Furthermore  we  assume  that 

E[aia.]  =  cr*  i  =  j 

=  0  i  i  j  ;  (21) 

cr ^  is  then  the  average  power  associated  with  each  plane  wave,  and  it  is 
assumed  that  any  two  plane  waves  are  uncorrelated.  Under  these  condi¬ 
tions,  H(r,t)  will  be  a  homogeneous  Gaussian  random  field. 

The  correlation  coefficient  described  in  Sec.  A. 2  has  an  interesting 
physical  interpretation:  the  correlation  of  two  complex  envelopes  at 
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different  points  in  space  and  at  some  additional  time  delay.  We  define 
the  function  /?ff(ArAt)  as  the  space- time  correlation  function 


/?H(Ar,  Af) 


E  [H(r,  t  )ll*  (  r  +  Ar,  t  +  At)] 


5 

i  =  i 


cr2e 


-  i  (  2  rt\  A  <  +  *  A  f  1 


(22) 


The  correlation  coefficient  is  simply  the  normalized  spacetime  covariance 
function  evaluated  at  zero  time  shift,  and  Ar is  the  spacing  between  an¬ 
tennas  located  at  x j,  y1?  Zj  and  x 2>  y2,  z 2, 


P 


( 0 , Ar  j  2 ) 


i  =  l 


cr2e 

l 


“  i  k 


A  r 


1  2 


(23) 


The  Fourier  transform  of  fiu(At,Ar)  is  the  scattering  function 

n 

S(\,kz,ky,kt)  =  jyj7flH(At,Ar)<ri  [2^At  +  k  •  A~]rfAtdAxc/AydAz  (24) 

where 

k  •  Ar  =  k  Ax  +  k  Ay  +  k  Az  (Tv) 


For  At , Ar)  of  the  form  of  Eq.  (22),  Eq.  (24)  becomes 


S  ( X ,  k  ,  k  ,  k  ) 

x  y  ’  z 


£  O'2  §(X  -  X  )  8(  k 
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Equation  (26)  shows  that  the  scattering  function  consists  ol  a  number  ol 
delta  functions. 


Let  us  define  a  new  function,  the  wave-number  dispersion  limit  ion 
H'(  k  t ,  k  ,  k  y  ,  k  { )  .  We  integrate  over-all  Doppler  shift  X  to  obtain  a  scat 
tering  function  in  wave  number  only 
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nkx,ky,kz)  =  jS(K,  k  x,  k  ^ ,  k  z)dK 


=  i  a2Hk  -  k  )Hk  ~  k  )Uk  -  k  )  .  (27) 

^  =  j  t  x  xi/vy  y  i  '  2  J  i 

Tlte  wave-number  dispersion  function  is  the  Fourier  transform  of  /?w(0,Ar;. 

It  indicates  how  the  channel  disperses  a  single  plane  wave  into  many  plane 
waves  of  varying  magnitude  and  direction.  Inany  principal  direction,  a 
slowly  converging  wave-number  dispersion  function  implies  a  rapidly  con¬ 
verging  spacial  autocorrelation  function  /?^(0,Ar).  This  implies  that  an¬ 
tennas  spaced  only  a  short  distance  apart  will  have  independent  complex 
envelopes  and  hence  that  one  need  not  increase  antenna  spacing  to  have 
effective  space- di versi ty  correlation.  Conversely,  a  narrow  wave-number 
dispersion  function  in  one  principal  direction  implies  a  wide  spacial 
autocorrelation  function  in  that  directiot  and  that  space  diversity  is 
rather  ineffective. 

The  results  of  this  section  are  applicable  only  where  there  is  no 
polarization-selective  fading  or  if  only  circularly  polarized  antennas 
are  used.  The  more  general  case,  where  time-varying  polarization  causes 
changes  in  amplitude  on  an  elliptically  polarized  antenna  is  not  treated 
here. 

C.  MEASUREMENTS 

In  this  section,  the  measurements  made  to  determine  normalized  co- 
variance  between  spaced  antennas  are  outlined.  No  effort  is  made  to 
determine  the  wave-number  dispersion  function.  Since  only  two  antennas 
were  available,  the  only  measurement  made  was  the  recording  of  correlation 
between  these  two  antennas  during  periods  of  error- rate  measurement. 

For  convenience,  it  was  decided  to  measure  the  correlation  between 
the  envelopes  of  the  transmitted  carrier  signal  as  received  on  two  receivers 
each  connected  to  its  respective  antenna.  The  two  receivers  were  Technical 
Material  Corporation  Model  GIR-9Q,  having  about  100-cps  IF  bandwidth.  The 
envelopes  were  obtained  by  monitoring  the  VG(  voltage  of  each  receiver. 

Th  c  amplitude  response  ol  the  \GC  is  411  as  i  -  1  oga  r  1  t  hm  i  c  and  is  approximately 
the  same  as  the  response  of  th*'  combination  of  H-ilQOA  111!!  and  CY-I.n? 
dt  ■scribed  in  Chapter  Ml.  Figure  78  shows  a  block  diagram  ol  the 
me as a  remen t . 
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0  -  «t?t  -  1*9 

FIG.  78  BLOCK  DIAGRAM  OF  OVER-ALL  CORRELATION  MEASUREMENT 
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1.  C0BBELAT10N  COMPIJTEB— .pHtNClPLES  Of  Opehation 

The  correlation  computer  (see  Fig.  79)  provides  an  output  that  is  an 
approximation  to  the  correlation  defined  in  Fq .  (1).  The  approximation 
consists  of  correlation  of  a  clipped  version  of  each  input  after  removal 
of  their  respective  means. 

The  input  signals  are  amplified  and  each  is  fed  into  u  l()0-.second 
integrator.  The  output  of  each  of  the  integrators  is  approximately  of 
the  form 
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X  (  t)  (It 

}o 


(28) 
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where  K  a  circuit  constant.  The  integrator  output  is  then  an  estimate 
of  the  mean  based  on  the  previous  100  seconds  of  input. 

The  output  of  each  integrator  is  subtracted  from  the  input  and  clipped 
at  a  very  low  level  by  a  two-sided  elippe  •.  The  clipper  gives  a  positive 
6  volts  if  the  signal  is  greater  than  its  and  a  negative  b  volts  if 

the  signal  is  less  than  its  mean. 

The  output  from  the  two  clippers  is  combined  in  an  '*  oxclusi  ve-o:” 
circuit.  This  circuit  gives  a  positive  output  when  both  inputs  are  of  the 
same  sign  and  a  negative  output  when  the  inputs  are  of  different  sign. 

Thus  if  input  signals  tend  to  deviate  from  their  mean  in  opposite  senses 
(p  <  0),  the  circuit  will  give  a  negative  voltage.  When  the  signals  de¬ 
viate  from  their  means  in  the  same  sense  (p  v-  0 ) . ,  the  output  will  be  posi¬ 
tive.  If  the  signals  are  uncorrelated  (p  =  O',  the  output  will  be 
positive  half  the  time  and  negative  the  other  half.  The  output  is  integrated 
with  a  time  constant o f  100  seconds  and  recorded.  An  example  of  this  record 
is  shown  lit  fig.  80  This  integrated  output  is  approximately  proportional 
to  t  he-  cor  re }  a  1 1  on  •coefficient  between  the  two  input  signals. 

2.  \  ctBAi  t  or  the  (.i  tpnrp  Cottar i  at i ov  \»*i*ho\ i -a*  t  i  o> 

■Vssume  the  input  signals  to  th*-  correlator  ate  ««*tne  turn  t  tons  di  f 
fr  ring  only  in  a  const  ant  phase  angle,  t* .  The  true  correlation  of  these 
two  signals  as  a  function  of  the  phase  wnglr  difference  t  a  l *o  called  the 
autocorrelation)  is  easily  found  Similar  >.  the  correlation  of  the 
c  tinned  verst  on  o  f  t  h  i  s  autorni  rel  at  ion  ha*  been  determined,  it  is  plotted 
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14  NOVEMBER  1964 


F!G.  80  EXAMPLE  OF  CORRELATION  COMPUTER  OUTPUT 


in  Fig.  81  along  with  the  true  autocorrelation.  The  greatest  difference 
occurs  at  t*  =  ti 74,  where  the  true  correlation  is  0.707  while  the  clipped 
approximation  is  0.5. 

For  a  Gaussian  random  process.  Van  \Teck~has  shown  that  if  rlf)  is 
the  original  process  and  jr.(i)  its  clipped  version,  then  the  normalized 
autocorrelation  function  of  the  original  process,  p  ,  in  terms  of  the 
clipped  version,  p^  is 


Thus  for  a  Gaussian  random  process,  a  deterministic  relationship  exists 
between  the  two  autocorrelation  functions. 


FIG.  81  CORRELATION  COMPUTER  APPROXIMATION  TO  TRUE  CORRELATION 
FOR  A  SINUSOID 
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I).  COMPARISON  OF  SPACE  - 1)  I  VERS  I TY  CORRELATION  AM)  TIME  SPREAD 

It  is  interesting  to  compare  the  hourly  trends  between  space-diversity 
correlation  and  time  spread.  In  Sec.  U.  3  it  was  shown  that  a  -large  spread 
in  wave  number  would  cause  antennas  spaced  only  a  short  distance  apart  tc 
have  uncorrelfited  envelopes.  If  it  is  further  assumed  that  a  situation  in 
which  modes  are  separated  in  wave  number  will  also  be  one  in  which  large 
separations  exist  in  time  delay,  then  a  large  time  spread  should  indicate 
low  correlation  and  vice  versa. 

The  apparent  validity  of  this  assumption  is  illustrated  in  Fig.  82 
for  7  Me  and  to  a  lesser  extent  in  Fig.  83  for  14  Me.  These  curves  con¬ 
sist  of  averages  over  each  5-minute  period  for  space-diversity  correlation 
9  October  to  14  November  1964  and  fur  rime  spread  from  4  to  14  November  1964. 
They  show  the  hourly  trends  of  each  parameter  and  thus  tend  to  average  out 
daily  fluctuations  from  these  trends. 

The  7-Mc  average  time  spread  and  space  correlation  shown  in  Fig.  82 
appear  to  be  almost  negat i ves  o f  each  other.  This  is  not  quite  as  apparent 
for  the  14-Mc  data.  A  possible  reason  for  this  difference  between  the  two 
frequencies  is  the  consistent  presence  of  a  greater  number  of  modes  on  the 
lower  frequency.  At  14  Me  much  of  the  uncorrelated  fadi ng  may  have  been  due 
to  differences  ia  polarization  between  the  two  ant ennas  and  eo r respond , ngly 
less  to  spread  in  wave  number.  It  should  be  noted  that  from  2200  to  2300 
the  14-Mc  time-delay  spread  dropped  to  about  lOC/rsec,  which  would  cor¬ 
respond  to  a  single  plane  wave.  The  space-diversity  correlation  did  not 
rise  as  fast  during  this  period,  nor  did  it  attain  a  com  on  of 
which  would  be  indicated  for  a  single  inode.  Investigators  such  as  Balser 
and  Smith"®  and  Hedlund  and  Edwards"'  have  noted  the  presence  of  polari¬ 
zation  fading  on  a  single  mode.  Furthermore,  with  such  a  narrow  wave- 
number  dispersion  function  as  would  be  present  with  a  single  mode,  the 
presence  of  receiver  noise  or  variation  in  antenna  patterns  would  also 
tend  to  lower  the  correlation. 

Figures  84  and  85  show  the  day-to-day  variance  of  space- di versi ty 
correlation  and  time  spread  for  both  7  and  14  Me.  The  variance  of  each 
was  computed  for  the  same  data  as  the  averages  cited  previously.  A  high 
variance  indicates  a  large  day-to-day  variation  for  the  particular 
parameter  at  the  particular  time,  while  a  small  variance  indicates  rel¬ 
atively  high  day-to-day  stability.  For  the  7-Mc  data,  it  is  clear  that 
the  change  in  variance  of  time  spread  with  time  of  day  is  opposite  to 
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VARIANCE  OF  TIME  SPREA.'  VARIANCE  OF  SPACE  CORRELATION 


FIG.  84  COMPARISON  OF  VARIANCES  OF  CORRELATION  AND  TIME-DELAY  SPREAD, 
7  Me 
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that  of  spare-diversity  correlation  and  that  the  change  in  vaiiance  of 
lu'tli  time  spi'end  and  space-diversity  correlation  with  time  of  day  is 
roughly  tlte  same  as  the  change  with  time  of  day  of  their  respective 
rages . 

It  may  be  deduced  that  a  stable  situation  as  far  as  time  spread  is 
concerned  (t.e.,  very  few  modes)  is  relatively  unstable  as  far  as  space- 
diversity  correlation  is  concerned.  If  few  modes  are  present  (stable 
variance),  a  change  in  elevation  or  azimuth  of  any  one  will  have  more 
drastic  effects  on  the  wave-number  dispersion  function  (and  hence  on  the 
space-diversity  correlation)  than  if  many  modes  were  present. 

The  tendency  of  tiie  aveiages  to  follow  the  variances  indicates  that 
day-to-day  variations  in  a  parameter  are  proportional  to  the  magnitude 
of  the  parameter  at  that  time.  This,  in  turn,  indicates  a  somewhat  con¬ 
sistent  measurement.  Again,  both  these  trends  are  somewhat  less  apparent 
at  11  Me  than  at  7  Me.  The  explanation  is  probably  the  same  as  that  for 
the  averages  discussed  earlier. 

Figure  86  shows  the  histograms  of  7-  and  14-Mc  space- di versi ty  cor¬ 
relation.  For  each  frequency,  the  number  of  occurrences  is  given  without 
regard  for  time  of  day.  The  7  Me  data  are,  on  the  average,  more  uncorre¬ 
lated  than  those  taken  at  14  Me,  although  the  1000- foot  spacing  between 
each  set  of  antennas  represented  about  7-1/2  wavelengths  ac  7  Me  and 
15  wavelengths  at  14  Me.  This  is  certainly  contradictory  to  much  engi¬ 
neering  practice,  which  assumes  that  the  greater  the  spacing  in  wavelength 
the  lower  is  the  correlation  between  antennas.  Again,  the  explanation 
probably  lies  in  the  fact  that  more  modes  were  present  in  the  7-Mc  data 
than  in  the  14-Mc . 

In  Sec.  A. 2,  it  was  pointed  out  that  the  correlation  coefficient 
between  the  envelopes  of  two  jointly  Gaussian  random  processes  should 
never  be  negative.  The  appearance  on  the  histograms  of  negative  correla¬ 
tions  may  indicate  that  5-minute  samples  were  not  long  enough  to  average 
out  the  1 ow- f requency  variations  of  the  processes.  It  is  felt  that  if  the 
sampling  period  had  been  increased  to  10  minutes  the  average  correlation 
would  have  been  negative. 
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E.  OCCURRENCE  OF  ERRORS  WITH  DUAL  DIVERSITY 

Figure  87  shows  the  binary  error  probability  lor  all  the  data  taken 
with  the  use  of  a  dual - di versi t y  system.  These  data  were  taken  from  9  to 
12  November  on  both  7  and  14  Me.  It  is  clear  that  almost  all  of  the  points 
fall  between  the  values  that  would  be  expected  if  no  correlation  existed 
(FSK  2)  and  those  that  would  be  obtained  if  complete  correlation  existed 
(FSK  1). 

On  12  November  errors  were  made  during  five  non- conti guous  30-second 
intervals.  Shown  in  Fig.  88  are  raw  data  (without  300-second  visual  in¬ 
tegration)  5  minutes  before  and  after  the  30-second  interval  in  question. 

In  each  case,  the  correlation  stayed  above  0.6  for  several  minutes  after 
the  errors  were  made.  Because  of  the  100-second  time  constant,  the  con¬ 
dition  of  the  channel  (for  space  correlation)  for  a  given  time  was  best 
obtained  by  observing  the  space- cor r el  at i on  data  for  at  least  100  seconds 
after  the  given  time.  These  results  indicate  that  errors  were  made  when 
the  correlation  coefficient  was  high.  It  should  not  be  assumed  that  the 
converse  is  true:  high  correlation  coefficient  does  not  cause  errors. 

The  most  probable  explanation  is  that  time-  or  f requency - sel ecti ve  fading 
is  more  likely  to  cause  errors  when  correlation  coefficient  is  high  than 
when  it  is  low. 
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FIG.  88  SPACE  CORRELATION  BEFORE  AND  AFTER  DETECTED  ERRORS, 
14  Me  12  NOVEMBER  1964 


VIII  CO?KXUSIO*>S 


A  mathematical  model  for  an  FSK  data  transmission  system  operat  i ng 
through  an  HF  channel  waS'developed  by  SHI  to  predict-  binary  error  rate. 

The  channel  portion  of  this  model  considers  SAY  as  well  as  time^delay 
and  Doppler- frequency  effects  introduced  by  thb  ch&niifcl.  The  system 
portion  of  the  model  incorporates  var  ious  real  istic  FSK  system  features 
including  a  finite  transition  time  between  the  modulation  .^states.  ' 

An  experiment  was  conducted  to  simu'l taneousby  measure  the  binary  error 
rate  of  an  AN/FGC- 29- type  FSK  system  antfthe  S/S,-  tiwje-deiay  spread  and  . 
Doppl er- frequency  spread  of  the  channel.  The  results  of  thSs  experiment 
are  summarized  in  Figs.  46  through:55.  The  asymptotic  bottoming  out  of . 
the  error-rate  curves  at  high  S/Ns  for  various  channel  conditions  is 
clearly  demonstrated.  :  -  ,  _  /  A 

For  normalized  time-delay  and  Doppler  spreads  less  than  0.05,  the  v 

asymptotic  error  rate  for  the  theoretical  model  is  much  i«ore  sensitive 
to  system  transition  time  than  to  changes  in  the  t ime-.de hay.  or  Doppler-  - 
f requenev  spreading  of.  the  channel-  The  model  shows  that  a  reduction  of  .? 
system  transition  time  from  200  to  y^sec  would  reduce' the  asymptotic 
error  rate  from  1  x  10  3  to  approximately  1  *  10  1  for  relatively  low 
spread  conditions  on  the  channel.  This  design  feature  should  receive 
careful  consideration  in  future  FSK  systems. 

Although  the  ranges  of  time-delay  and  Doppler  spread  observed  in  the 
experiment  were  small,  the  performance  of  the  experimental  system  showed 
a  well-defined  sensitivity  to  time-  and  f requency- se lect i ve  effects. 

Although  this  sensitivity  was  greater  than  that  predicted  by  the  model, 
it  was  reasonable  in  view  of  the  differences  between  the  detection  pro¬ 
cedure  of  the  modeled  system  and  that  of  the  experimental  system. 

Insufficient  data  were  collected  with  dua  1  - di ve r s i t y  operation  to 
permit  a  confident  comparison  with  the  theoretical  model.  However,  indi¬ 
cations  were  that  the  existence  of  correlated  space-diversity  channels 
would  produce  higher  measured  error  rates  than  predicted  for  the  inde¬ 
pendent  dua  1  -d i vers ity  model.  Measurements  of  antenna  correlation  have 
shown  appreciable  time  intervals  when  correlation  coefficients  exceeded  0,6. 
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Thr  alternative  iransaissiun  sequence  used  in  th»-  expei  icu-iit  was 
chosen  to  maximize  the  f  requeue v- s»  led  i  ve  effects.  Kith  a  more*  random 
frasswisfioii  sequence,  representative  of  normal  traffic,  it  could  be  ex¬ 
pected  that  the  irreducible  error  rate  would  be  lower  than  measured  here- 
The  109-baud  transmission  rate  used  in  this  experiment  was  slightly  higher 
than  the  75-baud  rate  used  operationally  with  this  system;  the  irreducible 
error  rate  would  be  slightly  lower  when  used  at  the  lower  rate. 

For  a  shorter  HF  propagation  path  than  the  4iC0-km  circuit  used  in 
this  experiment,  larger  time-delay  spreads  and  hence  a  higher  irreducible 
error  rate  could  be  expected.  Similarly,  on  paths  through  auroral  activity 
or  other  highly  turbulent  ionospheric  conditions,  a  higher  irreducible 
error  rate  would  be  likely.  A  longer  path  would  probably  display  smaller 
variations  in  time-delay  and  Doppler  spread. 

In  general,  the  results  of  the  experiment  indicate  the  importance  of 
the  dispersive  channel  effects  in  limiting  the  performance  o*  a  communi ca¬ 
tions  system.  For  the  range  of  spread  conditions  observed  cn  the  channel, 
the  model  yields  a  good  fit  to  the  average  performance  o!  the  experimental 
system  over  all  dispersive  states  observed.  However,  to  obtain  a  mere  ac¬ 
curate  error- rate  estimation  for  Lower  and  higher  channe 1 -spread  condi¬ 
tions,  it  may  be  necessary  to  include  in  the  model  nomnstantaneous  energy- 
quenching  and  other  filtering  characteristics  of  the  particular  experimental 
system. 

The  adequacy  of  this,  or  any,  error- rate  model  depends  on  its  intended 
application.  If  this  theoretical  FSK  model  is  to  be  used  to  compare  the 
performance  of  the  AN/FGC-29  system  with  the  performance  predicted  by  the 
theoretical  model  of  another  system  under  a  wide  range  c.f  channel  condi¬ 
tions,  it  may  be  necessary  to  model  the  FSK  system  more  closely.  The  re¬ 
sults  of  the  experiment  give  positive  indications  that  the  present  FSK 
model  can  be  modified  to  yield  an  adequate  model  of  the  AN/FGC-29  system 
for  arbitrary  dispersive  states  of  the  channel.  Consideration  should  be 
given,  however,  to  the  variation  of  critical  system  parameters  among 
AN/FGC-29  systems. 

If,  however,  the  model  is  to  be  used  to  estimate  outage  time  of  a 
circuit — based  on  some  modest  user-defined  minimum  preformance  threshold — 
then  the  model  is  very  likely  to  be  adequate,  especially  if  the  variations 
in  propagation  conditions  are  extreme.  Redundant  - 1 anguage  traffic  typi¬ 
cally  requires  binary  error  rates  of  10-3  for  intelligibility.  Probably 
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the  high-accuracy  requirement  for  data  transmission  will  In-  achieved  in 
current  HF  systems  by  the  use  of  coding.  Such  systems  * ■ I  I  still  operate 
with  channel  binary  error  rates  in  the  10' 3  range.  in  a  rase  where  one 
is  interested  only  in  predicting  whether  the  binary  error  rate  is  above 
or  belcw  that  required  for  mini'-um  performance,  the  range  of  interest  in 
predicted  error  rates  will  be  limited  to  moderate  values,  and  the  present 
system  model  will  suffice. 
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APPENDIX  A 

ESTIMATION  OF  RMS  VOLTAGE  LEVELS 
WITH  AN  ENVELOPE  DETECTOR 
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APPENDIX  A 


ESTIMATION  OF  RMS  VOLTAGE  LEVELS 
WITH  AN  ENVELOPE  DETECTOR 


Let  x(t)  be  a  narrowband  random  process  with  zero  mean.  Consider 
the  estimation  of  the  power  in  the  process  by  the  detection  scheme  shown 

in  Fig.  A- 1 . 
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FIG.  A-l  ENVELOPE  DETECTOR-FILTER  COMBINATION 


The  input  may  be  written  in  the  form 

x(t )  =  V{t)  cos  [ 2vt/ 0 e  +  (pit)] 

where  the  envelope,  F(t)  >  0  for  all  t.  A  measure  of  the  power  in  the 
input  signal,  x(t) ,  is  the  variance,  a2x(t)- 

The  output  of  the  half-wave  linear  rectifier,  yit),  is 

y(t)  -  V(  t )  cos  [2  0 1  +  0(t)l  when  cos  [  2'77/  0  f  +  <p{  t )  >  0 

“  0  otherwise 

By  a  Fourier  series  expansion,  yit)  may  be  written  as 

y(  t )  ■  V(t )  ^60  +  6j  cos  [2 nf0t  +  (p{t)]  +  6 2  cos  (2  [2 nf  Qt  +  4i  t  )J }  +  . . 

Since  z(t)  is  a  low-pass  filtered  version  of  y(t)  with  a  cutoff  frequency 
much  smaller  than  /0,  then 
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si  *)  -  U0K(t)]*/»(<) 


j 


:  bn  J  V(T)hi t  -  r)dr 

-QD 


Without  any  loss  in  generality,  both  6U  and  J  h( t  ~  T)d’  will  be  » hosen 

-CD 

equal  to  one,  since  these  are  arbitrary  gain  constants.  The  value  of 
:{t)  will  be  chosen  as  an  estimate  of  cr  (  t ) ,  the  rrns  value  of  the  input 
signal  process.  The  expected  value  of  the  estimate  is 

7(t)  =  j"  V(r  )h(t  ~  r)dr  .  ( A-  1 ) 

-00 

If  the  mean  of  V(t)  does  not  vary  appreciably  during  the  memory  of  the 
low-pass  filter,  then  Eq.  ( A- 1 )  reduces  to 

7(t)  =  Fit) 

For  a  simple  RC  filter,  a  realistic  value  of  the  filter  memory  is 
on  the  order  of  two  or  three  time  constants. 

It  will  be  shown  that  for  some  processes  the  square  root  of  the 
input  power  and  the  expected  value  of  the  input  signal  envelope  are  ap¬ 
proximately  proportional, 


F(t)  -  Kcrjt)  (A- 2) 

where  czit)  is  the  rrns  value  of  the  input  signal  process  at  time  t  and  A 
is  a  constant  dependent  upon  the  process.  The  functional  dependence  upon 
t  is  retained  to  include  nonstati onary  processes.  When  power  ratios  are 
formed,  as  in  S/N  calculations,  the  constants  will  cancel  if  the  noise 
and  signal  plus  noise  processes  have  the  same  constant.  If  the  constants 
for  the  two  processes  are  not  identical,  then  an  additional  bias  is 
introduced  into  the  result. 

To  illustrate  possible  values  of  K,  the  expected  value  of  the  random 
variable  F  may  be  obtained  for  various  assumed  distributions.  li  l  has 
a  Rayleigh  distribution,  then 

F  =  Vv/2  a  ,  i.e . ,  A  =  ^n/2 
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If  the  input  signal  is  a  sinusoid  of  fixed  amplitude  A,  the  output, 
z(t),  will  be  a  constant  value  equal  to  A.  The  output  is  thus  equal  to 
v/2  times  the  rms  value  of  the  input  signal. 

If  V  has  a  Rice  distribution,  and  <r xit),  the  rms  value  of  x(t),  is 
lU2/2)  +  cr2]'^;  then28 


where  is  a  hypergeometric  function.  If 


A 2 

—  >10 
2cr  2 


then,  by  using  the  first  two  terms  in  the  expansion  of  the  hypergeometric 
function,  the  following  result  may  be  obtained. 


In  this  case  we  do  not  have  a  linear  relationship  between  V  and  cr x,  the 
rms  value  of  x(t).  However,  under  the  conditions  of  the  approximation, 
the  relationship  is  roughly  true  with  the  constant  of  proportionality 
equal  to  v 2  • 


as 


It  is  interesting  to  note  that  a  parameter  designated  V.  and  defined 


161 


has  been  extensively  investigated  experiments i ly  for  HF  atmospheric:  noise 
by  the  National  Bureau  of  Standards.®  The  parameter  is  related  to  K  as 
given  in  Eq.  (A-21  by 

,  .  »...  (?)  ■  . 

The  NBS  work  indicates  that  the  value  of  lies  between  2  and  4  for  the 
frequency  range  7  to  15  Me  with  the  variation  occurring  between  winter 
daytime  and  summer  nighttime.  These  values  correspond  to  values  of  K 
equal  to  1. 12-  and  0.89,  respectively.  The  values  are  strictly  valio  only 
for  a  particular  year  and  the  site  near  Boulder,  Colorado;  however,  they 
give  some  indication  of  the. values  of  K  that  might  be  expected  for  atmo¬ 
spheric  noise. 

As  a  result  of  higher  values  of  K  for  signal  measurement  than  for 
noise  measurement,  calculated  S/N  assuming  equal  values  of  K  may  be  as 
much  as  4  db  high.  This  is  the  worst-case  estimate;  however,  the.  calcu¬ 
lation  may  be  2  and  possibly  3  db  high. 
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APPENDIX  B 


ESTIMATION  OF  POWER  AND  POWER  RATIOS 


It  has  been  shown  in  Appendix  A  that  if  the  time  variations  of  the 
input  envelope  are  slow  compared  to  the  filter  memory,  then  the  expected 
output  of  a  linear  detector- fi  1  ter  combination,  “(f),  i  .s  the  expected 
value  of  the  envelope,  V*  ( f ) ,  of  the  input  signal.  Further,  it  has  been 
shown  that  the  square  of  the  expected  value  of  the  envelope  is  approxi¬ 
mately  proportional  to  the  power  in  the  input  signal. 

Suppose  we  wish  to  obtain  an  estimate  of  the  average  power  in  the 
signal  over  a  time  interval  T'  ■  Assume  that  a  series  of  n  independent 

sample  values,  2 ( f  j ) ,  z(t0) . z ( t n ) ,  are  available.  The  actual  power 

(the  time  average  of  the  ensemble  a.  erage  power),  P  T>  ,  is  given  by 

1  r' 

<P>  r'  *  I  *2(t)dt 

1  0 

and  by  approximating  the  integration  by  a  summation  becomes 

1  * 

<P>  j '  ~  —  z.  a  2  ( t  (  )  .  t  t  1  n  ( 0 .  T 

"  .=  1 

where  o  "(f)  is  the  ensemble  variance  of  the  signal  process  * ( f  >  at  time  l 
For  a  stationary  process 


2 

r 


t' 

I 

0 


1 

9  f  * 


It  has  been  shown  that  for  the  cases  «>f  interest 


1  ( i  \  ’  Hr)  *  t  > 
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Define  an  error  term 


which  has  zero  mean  for  the  stationary  case.  Assuming  that  the  errors 
are  independently,  jointly  normally  distributed,  the  joint  probability 
density  function -is 


( B-  1 ) 


where  z{t_) 
K 


is  defined  as  z 

l 

known  to  be  the 


Given  (z  z  ),  the  best  estimate  of 

1  ft 

maximum- 1  i ke  1  i hood  estimate.  50  This  estimate 


<*  v^7) 


is  determined Iby  maximizing  the  likelihood  function 


for  each  ( z  1 ,  z  z  n  ) .  For  the  function  of  Eq.  (B-l),  this  is  most 

easily  achieved  by  maximizing  tlie  natural  logarithm  of  the  function.  The 
resulting  maximum- 1 i k el ihoo d  estimate  becomes 


and  the  best  estimate  of  the  power  becomes 


It  is  often  necessary  to  estimate  S/Ns  from  the  power  estimate  of 
S  +  N  and  noise.  Define  a  measure  of  the  power  as  H  where 
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1  " 

-  1  2 

n  .=  ) 


.] 


For  S  +  N  and  for  noise  only  present,  respectively,  wr-  have 


If 


S+JV 


1/2  n2 
n  5+  V  5+ N 


*;  ■  *k 


Since  the  signal  and  noise  are  independent  and  zero  mean  processes, 
the  total  power  is  equal  to  the  sum  of  the  individual  powers, 


^IrN 


aS  +  aN 


The  true  S/N  is  given  by 


S 

~N 


a 


a 


s  +l¥ 


-  1 


y ; 


*!♦*  % 


Since  we  have  insufficient  information  to  specify  the  ratio  ( K v  fci+v) 
with  any  degree  of  confidence,  we  choose  as  our  S/N  estimate 


A 

5 


.s+,v 
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The  expected  value  of  the  estimate  is 


where  cN  represents  the  variation  of  the  noise  power  estimate  about  its 
mean.  Since  the  noise  power  estimate  will  always  be  greater  than  zero, 


With  the  assumption  that  the  sum, 


is  negligible  compared  to  one,  we  have  approximately 
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Normally  the  last  term  will  be  insignificant,  since  K  is  of  the  same 
order  of  magnitude  as  KN-  The  multiplicative  constant,  {K  s.t  N/K  „)2 ,  in 
the  first  term  introduces  an  additive  bias  when  the  S/N  is  expressed  in 
decibels.  The  constant  will  normally  he  greater  than  one  and  increase 
with  S/N  (see  Appendix  A).  Tbs  resulting  estimate  may  thus  be  2  or  3  db 
high. 
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APPENDIX  C 


ERROR  ANALYSIS  OF  S/N  ESTIMATIONS 


In  the  formulation  of  the  S/N  in  decibels,  the  logarithm  of  a  power 
ratio  is  formed.  In  the  error  analysis,  theeffert  of  errors  in  the 
power  measurement  upon  the  S/N  ( expressed  ir.  decibels)  must  be  considered. 
In  the  material  that  follows,  the  error  propagation  effects  are  analyzed. 

First,  consider  the  expansion  of  the  estimate  in  a  Taylor  series 
about  the  true  value: 


x  +  e 


;  - 


n=  ]L  \y, 


In  +  1  ~  exp 


where 


x  +  e 


y  +  S  y 


(y  +  ey) 


Therefore , 


x  +  e 


y  +  e 


Inf— 


(7)  t  1 


(y  +  (.y  +  fcy)- 


If  the  error  of  the  original  measurements  may  be  bounded  in  the  following 


3  e 
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that  is,  i f  we  can  bound  both  terms  with  a  fixed  percentage  error,  then 


1  ~  exp 


:*) 


<  1  -  exp 


yt 


+  £ 


(  y  +  C  )  x 
J  y 


exp 


In  the  Taylor  series  expansion  it  has  been  convenient  t < •  use  natural  log¬ 
arithms:  an  error  bound  may  be  given  in  decibels  that  is  independent  of 

the  S/N  and  depends  only  upon  the  percentage  error,  e.  assumed  for  the 
original  power  measurements.  The  lower  bound  is 


10 


In  10 


1  -  exp 


2e 


.1  +  e  / 


and  the  upper  bound  is 


10 


In  10 


1  “  exp 


1  -  e 


These  error  bounds  are  plotted  as  a  function  of  the  percentage  error  in 
the  original  power  estimate  in  Fig.  C-l.  To  illustrate,  an  error  of 
±25  percent  in  the  power  measurement  (e  =  0.25)  will  result  in,  at  most, 
an  error  of  ±2  db  in  the  S/N. 
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APPEND TX  D 


CONFIDENCE  LIMITS  IN  THE  MEASUREMENT 
OF  LOW  ERROR  PROBABILITIES 

1.  INTRODUCTION 

The  experimental  validation  of  an  error-rate  model  for  a  communica¬ 
tions  system  requires  an  estimate  of  the  true  binary  error  rate  from  a 
finite  set  of  measurements.  Two  questions  then  arise:  How  much  confi¬ 
dence  can  be  placed  in  error-rate  estimates  made  with  limited  measured 
data?  How  much  measuring  must  be  done  to  obtain  a  given  degree  of  con¬ 
fidence  in  an  estimate? 

The  general  technique  for  measuring  the  binary  error  rate  of  a  com¬ 
munications  system  is  to  transmit  a  known  binary  sequence  over  the 
channel,  compare  the  receiver’s  decision  with  a  local  reference,  and 
compute  the  ratio  of  the  number  of  errors  to  the  number  of  transmitted 
bits  for  some  interval.  The  receiver  decides,  only  through  its  normal 
detection  scheme,  which  binary  signal  was  transmitted.  After  some  mea¬ 
suring  interval,  we  are  asked  to  give  an  estimate  of  the  true  binary 
error  rate  for  the  system  from  this  data.  Because  our  estimate  must  be 
based  on  a  finite  number  of  samples  of  a  random  phenomenon,  it  will  be 
subject  tc  some  error.  If  we  can  describe  this  error  in  terms  of  a 
probability  law,  we  can  better  assess  the  accuracy  of  our  answer.  When 
the  true  error  rate  is  small,  we  may  obtain  a  very  few  errors  in  our 
sample.  We  then  have  the  classical  problem  of  estimating  the  rate  of 
occurrence  of  rare  events  from  a  small  sample. 

2.  ERROR- RATE  PROBABILITY  LAW 

A  binary  communications  system  transmits  one  of  two  possible  signal 
(mark  or  space),  and  the  receiver  must  decide  at  the  end  ul  each  signal¬ 
ing  element  time  which  signal  was  transmitted.  These  receiver  derisions 
can  be  described  by  Bernoulli  trial*  where  the  results  are  traditionally 
termed  success  or  failure.  In  this  Appendix,  the  occurrence  of  an  error 
will  be  designated  by  success  and  the  absence  of  an  error  bv  failure. 
1’his  definition  is  consistent  with  statistical  conventions. 
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For  a  symmetrical  binary  cliatmel,  we  define  p  us  the  probability  of 
a  detected  error,  where  a  detected  error  is  defined  as  a  mark  r» reived 
when  a  space  was  transmitted  or  a  space  received  when  a  mark  was  trans¬ 
mitted.  We  desire  to  estimate  the  parameter  p  from  our  measurements. 

3.  LOW- ERROR  CONFIDENCE  INTERVALS 

From  a  Bernoulli  probability  distribution,  i f  we  take  n  independent 
samples  where  the  probability  of  error  (or  success)  is  p,  the  probability 
of  obtaining  x  errors  will  follow  a  binomial  distribution.  (The  case  when 
the  samples  are  not  independent  is  discussed  later.)  II  we  know  a  priori 
that  the  probability  of  error  is  low  (i.e.  p  <  0.1),  we  can  approximate 
this  adequately  by  a  Poisson  distribution. 

We  then  find  confidence  intervals  for  this  distribution  by  determin¬ 
ing  probabilities  that  k  j  or  less  errors  occu.  and  k ,,  or  more  errors 
occur,  given  that  the  true  error  rate  is  p: 


Pix  < 


e~np  i  np)* 


P[x  > 


^  g~nP(nP)' 

*  =k  , 


where  n  is  the  number  of  samples  taken. 

Note  that  the  product,  np,  is  also  the  expected  number  oi  errors  we 
would  observe  with  this  error  rate  and  number  ol  samples.  I  he  expected 
number  of  errors  (or  successes)  is  more  convent  ioi*u  1  !  v  defined  as  the 
parameter  of  the  Poisson  distribution. 

If  we  have  observed  a  particular  number  of  errors  in  an  experiment 
hut  do  not  know  the  value  of  the  parameter,  we  can  reverse  the  procedure 
and  determine  ranges  of  parameters  that  could  still  give  the  observed 
value  within  a  given  probability. 
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From  cumulative  Poisson  distribution  tables  we  determine  the 
fol lowing: 

(1)  The  smallest  integf.  ilue  of  the  parameter,  np  such 
that  with  k  observed  errors 


*  =  0 


(2)  The  largest  integer  value  of  the  parameter,  np2,  such 
that  with  k  observed  errors 

V-  e”"*  2 ( «p 2  )  * 

L - ; —  i  c 

x : 

x  =  h 


where  C  is  some  confidence  level  (e.g.,  0.05). 

If  we  observe  k  errors  in  a  sample  of  n,  the  best  estimator,  p,  of 
the  true  binary  error  rate  is  simply  the  ratio  k/n.  The  upper  and  lower 
confidence  limits  at  some  confidence  level,  C,  art-  p  j  and  p  y  Valuej  of 
rip  j  and  np2  for  a  confidence  level  of  5  percent  are  shown  in  Table  D-  I . 

4.  CONFIDENCE  CRITEHIA 

It  is  convenient  at  this  stage  to  seek  a  s i ng i e- va l ued  confidence 
parameter.  For  example,  we  could  def>*r  a  ratio,  ,4,  in  terms  of  the 
upper  and  lower  confidence  estimates  as 


4 


Pi 


Or  «c  i- on  Id  define  a  percentage  error  /'.  lienern  tl»e  <  on  I  i  d«-n<  <  estimates 
and  the  estimate  itself  as  follovs: 


Itl 
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Tui>  !*■  n-i 

1'PPKR  AND  LOWED  90- PERCENT  CONFIDENCE  ESTIMATES 
OF  POl SSON  PARAMETER  WITH  k  OBSERVED  SUCCESSES 


* 

UPPER  CONFIDENCE  LEVEL* 

<nj»j) 

LOWER  CONFIDENCE  LEVEL* 

RATIO 
(4  ) 

PERCENT ACE  ERROR 
l/'.’ 

0 

3.0 

-- 

-- 

-- 

1 

4.75 

0.052 

91.3 

170 

2 

6.3 

0.36 

17.5 

JO" 

3 

7.76 

0.82 

9.46 

231 

4 

9.15 

1.38 

6 .  ti  3 

191 

5 

10.5 

1.97 

5.33 

171 

6 

11.85 

2.61 

4. 54 

151 

7 

13.15 

3.29 

4.0 

141 

8 

14.43 

3.63 

3.62 

111 

9 

15.8 

4.7 

3.  3ti 

123 

10 

17. 

5.4-1 

3.12 

116 

it 

18.2 

6.18 

2.94 

109 

12 

19.5 

6.94 

2.81 

105 

13 

20.8 

7. 7 

O  7 

101 

14 

21.9 

8.46 

2.  'n 

15 

23.1 

9.25 

2.50 

92 

r. 

25.8 

10.85 

2.36 

88 

’O 

29.2 

13.27 

2.20 

80 

25 

34.9 

17.3 

2.02 

70 

30 

40.6 

22.4 

1.61 

6 1 

50 

64. 

i9.9 

1.60 

48 

70 

86.5 

58.  i 

1.49 

41 

100 

118.8 

85.6 

1.39 

33 

200 

228. 

182. 

1.25 

23  ■- 

300 

335. 

279. 

1.20 

19 

50f 

543. 

467. 

1. 16 

15 

10 

1063. 

961. 

1. 11 

10 

Confidence  terel  -  0.0$. 


Both  of  these  paraniters  are  also  tabulated  in  Table  D- 1 .  Notice  that 
these  confidence  criteria  depend  directly  on  the  number  of  errors  observed 
and  not  on  the  number  of  samples  taken. 

We  .ire  now  in  a  position  to  determine  the  amount  of  sampling  we  must 
do  to  obtain  a  given  degree  of  confidence  in  our  parameter  estimate.  If 
we  use  one  of  the  above  confidence  definitions,  our  rule  is  to  continue 
to  sample  until  the  observed  errors  are  greater  than  a  required  number. 
Note  chat  as  the  sample  size  is  always  increased  until  the  number  ol  ob¬ 
served  errors  exceeds  a  constant,  then  the  number  of  samples  required  will 
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be  roughly  inversely  proportional  to  the  error- rate  estimate.  This  is 
intuitively  satisfying  in  that  more  samples  are  required  for  lower 
error-rate  estimates. 

5 .  EXAMPLE 

Assume  we  desire  to  state  with  90-percent  confidence  that  the  dif¬ 
ference  between  the  upper  and  lower  confidence  levels  is  always  less 
than  the  estimate  itself.  From  Table  0- I  we  see  that  this  would  require 
a  minimum  of  14  observed  errors  in  our  sample.  If  the  number  of  observed 
errors  in  a  string  of  1000  transmitted  bits  were  6,  we  would  defer  making 
an  error  probability  estimate  until  we  observed,  say,  an  additional 
1000  bits.  If  the  next  1000  bits  contained  8  errors  we  would  estimate 
the  error  probability  as  14/2000  or  7  x  10-3.  And  we  could  say  with  90- 
percent  confidence  that  the  true  error  probability  lies  within  4. 2 x  10  3 
and  1. 1 x 10“2. 

6-  CORRECTIONS  FOR  NONINDEPENDENT  ERRORS 

It  has  been  observed  at  HF  that  errors  do  not  occur  independently. 
However,  without  a  detailed  knowledge  of  the  exact  statistics  of  the 
error  distributions,  we  assume  that  the  errors  occur  in  average  bursts 
of  length,  L,  and  that  the  bursts  of  errors  are  independent.  To  obtain 
confidence  limits  we  consider  that  an  “error"  is  in  reality  a  burst  of  L 
errors  so  that  the  minimum  number  of  actual  errors  required  for  a  given 
confidence  is  L  <.imes  the  previously  defined  errors. 

In  the  previous  example,  assume  that  the  average  burst  is  7  errors 
in  length.  To  observe  14  independent  bursts  of  errors,  we  would  need  to 
accumulate  slightly  over  100  errors  before  obtaining  the  degree  of  con- 
f idenc e  des i red . 


APPENDIX  E 


TIME  INTERVAL  CLASSIFICATIONS  OF  CHANNEL  ACCORDING  TO 
MEASURED  TIME-DELAY-  AND  DOPPLER- FREQUENCY -SPREAD  CONDITIONS 


APPENDIX  E 


TIME  INTERVAL  CLASSIFICATIONS  OF  CHANNEL  ACCORDING  TO 
MEASURED  TIME -DELAY-  AND  DOPPLER- FREQUENCY - SPREAD  CONDITIONS 


TABLE 

TIME-DELAY- SPREAD  RANGE 
(msec ) 

— 

DOPPLER- FREQUENCY-  SPRE  AD  RANGE 

(  cp  s  ) 

1 

0  -  0.25 

0  -  0.25 

2 

1 

0.2')  -  0.5 

3 

V 

0.3  -  l.‘i 

4 

0.25  -  0.75 

0  -  u.  25 

5 

| 

0.25  -  0.5 

6 

V 

0.5  -  I.o 

7 

0.75  -  1.5 

0  -  0. 25 

8 

1 

0.2.3  -  0.5 

9 

V 

0.5  -  1.0 

187 


V fi « j>^«.  ^.- j|R.xr-^as;.  .* 


DATA  TIME  INtERYALS  TABLE 


• 

4 

11061  0 

TO 

11- 

4 

1034(31 

4 

• 

4 

1634*31 

TO 

11- 

4 

1030 1  0 

4 

m 

4 

10381  0 

TO 

11- 

4 

1059131 

4 

m 

4 

1859131 

TO 

11- 

4 

1904131 

4 

m 

4 

1 V  0  4  l  3 1 

TO 

11- 

4 

1909(30 

4 

m 

4 

1909130 

TO 

11- 

4 

1914(30 

4 

m 

4 

1914130 

TO 
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APPENDIX  F 

TABULATED  RESULTS  OF  MEASURED  S/N  AND  BINARY 
ERROR  RATE  CLASSIFIED  BY  TIME- RELAY*  AND 
DOPPLER- FREQUENCY -SPREAD  CONDITIONS 


DESCRIPTION  OF  COLUMN  HEADINGS 


MAXIMUM  Maximum  possible  measured  S/N  in  any  30-second 
S/N  data  block  in  this  data  bracket. 

AVERAGE  Weighted  average  S/N  of  data  in  this  bracket. 

S/N 

BER  Estimated  binary  error  rate  for  data  accumulated 

in  this  S/N  bracket. 


TOTAL 


TOTAL 

ERRORS 


Total  number  of  transmitted  bits  accumulated  in 
this  bracket. 

Total  number  of  bit  errors  accumulated  in  this 
bracket . 
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TI«E  DELAY  SPREAD*  0,00  IU  1.50  MS 
DOPPLER  SPREAD*  0,00  ft)  l.ou  CPS 
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162331 

653 

31.5 

31,0 

3.369-03 

182331 

613 

32.5 

32.0 

3.326-03 

186833 

620 

33.5 

33.0 

3.149-03 

218347 

685 

34.5 

34.0 

2.979-03 

229602 

682 

35.5 

35.0 

1.269-03 

344403 

435 

36.5 

36.0 

2.659-03 

357909 

948 

37.5 

37.0 

2.339-03 

330897 

770 

30.5 

36.0 

2.0*9-03 

389423 

612 

39.5 

39.0 

1.879-03 

454702 

649 

40.5 

40.0 

1,589-03 

468206 

738 

41.5 

41.0 

1.789-03 

560499 

996 

42.5 

42.0 

1.729-03 

549244 

945 

43.5 

43.0 

2.049-03 

569503 

1162 

*4.5 

44.0 

i. 659-03 

650539 

1076 

*5.5 

45.0 

1 « 369-03 

736326 

1007 

46.5 

46,0 

1.419-03 

700061 

987 

47.5 

47.0 

1.389-03 

621615 

1135 

48,5 

46.0 

1.249*03 

745599 

1014 

49,5 

49,0 

1.209-03 

799105 

960 

50.5 

50.0 

4.179-04 

871137 

799 

51.5 

51.0 

7.419-04 

745041 

589 

52,5 

52.0 

7.759-04 

601017 

466 

53.5 

53.0 

1.009-03 

562750 

563 

54.5 

5**0 

*.219-0* 

398427 

367 

55.5 

35.0 

7,929-04 

3i:i48 

264 

56.5 

56.0 

5.469-04 

144562 

110 

57.5 

37.0 

1.129-01 

114601 

129 

59.5 

56.4 

6.959-04 

146315 

131 

6*o 

60.7 

5.719-04 

31514 

18 

TIME  DELAY  SPKEADi  O.OU  10  0.25  Mb 
DOPPLER  SPREAD!  0.00  TO  0.25  CPS 


MAXIMUM 

average 

8ER 

TOTAL 

TOTAL 

S/N 

S/N 

RITS 

ERRORS 

1A.5 

12.3 

3.299-02 

6753 

222 

16.5 

15.7 

9,779-03 

13506 

132 

17.5 

17.0 

1.649-02 

9004 

166 

16.5 

18.0 

9.569-03 

20259 

194 

19.5 

19.0 

7. 029-03 

22510 

156 

21.5 

20.6 

4.919-03 

36016 

177 

22.5 

22.0 

i. 269-03 

36016 

118 

23.5 

23.0 

4.079-03 

42769 

17A 

24.5 

24.0 

A. 669-03 

45020 

210 

25.5 

25.0 

5.379-03 

27012 

145 

26.5 

26.0 

6.4A9-03 

40518 

261 

26.5 

27.5 

2.559*03 

90040 

230 

29.5 

29.0 

5.149-03 

47271 

243 

30.5 

30.0 

2.729-03 

38267 

104 

31.5 

31.0 

4.229-03 

27012 

114 

32.5 

32.0 

3.299-03 

38267 

126 

35.5 

34.2 

1.279-03 

78785 

100 

42.5 

38.6 

4,799-04 

229602 

110 

47.5 

44.9 

4.979-04 

227351 

113 

52.5 

50.0 

4.469-04 

245359 

110 

60.5 

54.1 

6.429-04 

42769 

36 

| 

.1 

I 
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TI»£  DELAY  SPREAD!  0.00  TU  0.25  MS 
DOPPLER  SPREAD!  0,2b  TU  0.50  CPS 


MAXIMUM 

S/N 

AVERAGE 

S/N 

BER 

TOTAL 

BITS 

TOTAL 

ERRORS 

0.5 

0.0 

4.056-01 

20259 

8215 

1.5 

1.0 

2.536*01 

9004 

2278 

2.5 

2.0 

3.13P-01 

4502 

1410 

3.5 

3.0 

2.706-01 

20259 

5476 

*.5 

4.0 

3.616-01 

9004 

3428 

5.5 

5.0 

2.576-01 

15757 

4049 

6.5 

6.0 

2.356-01 

15757 

3704 

7  .5 

7.0 

1.206-01 

40510 

4066 

6.5 

6.0 

1 .646-01 

33765 

5554 

9.5 

9.0 

1.906-01 

11255 

2134 

10.5 

10.0 

7.246-02 

13506 

978 

lx. 5 

13.1 

5.016-02 

15757 

789 

15.5 

15.0 

4.206-02 

4502 

189 

17.5 

17.0 

4.936-02 

2251 

111 

16.5 

18.0 

4.076-02 

6753 

207 

24.5 

21.7 

1.506-02 

6753 

101 

32.5 

29.6 

2. 776-. 3 

56275 

156 

34.5 

33.6 

3.176-03 

36016 

114 

37.5 

36.3 

1.966-03 

72032 

141 

36.5 

38.0 

3. 006-03 

36016 

137 

39.5 

39.0 

1.666-03 

60777 

101 

40,5 

40.0 

2.146-03 

51773 

111 

41.5 

41.0 

1.796-03 

56275 

101 

43.5 

42.6 

1.126-03 

U6. ;5 

164 

44.5 

44,0 

1,636-03 

61036 

132 

45.5 

45.0 

1,266-03 

94542 

119 

47.5 

46.5 

7.456-04 

13956? 

104 

51.5 

49,« 

3,936-04 

330697 

130 

54.5 

52.9 

4,026-04 

213645 

103 

59.5 

55.9 

4,446-04 

92291 

41 
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TI«E  DELAY  sPHEADi  O.OU  10  0.2*  Mi 
DQPPLEfl  SPREAD*  0.50  fU  1 .00  CPS 


MAM  MUM 
S/M 

AVERAGE 

S/M 

B£R 

TOTAL 

BITS 

TfliAL 

ERrtORS 

19.  5 

17.6 

7.7 09-03 

20259 

156 

20.5 

20.0 

1 .00**02 

11255 

113 

21.5 

21.0 

2.199-02 

11255 

246 

23.5 

22.6 

1.049-02 

22510 

235 

25.5 

24.8 

5.269-03 

45020 

237 

26.5 

26.0 

5,551**03 

36016 

200 

26.5 

27,6 

5, 769*03 

60777 

230 

29.5 

29.0 

<i.  566*03 

29263 

134 

31.5 

30.5 

2.649-03 

67530 

178 

33,5 

32.5 

2.949-03 

56526 

172 

34.5 

34.0 

2.129-03 

47271 

100 

36,5 

35.5 

3.049-03 

123805 

376 

37.5 

37.0 

4.249-03 

51773 

166 

39.5 

3«.5 

2.059-03 

99044 

203 

41.5 

40.6 

1.1M9-03 

132809 

157 

43.5 

42.4 

1  *  379-03 

83267 

114 

45.5 

44.5 

1.099-03 

101295 

110 

52.5 

47.3 

1.039-03 

56526 

60 

31$ 


mm 


TINC  DELAY  SPRCAOl 
DOPPLER  SPREAD! 


0.25  TO  0.75  MS 
0*25  TO  O.Su  CPS 


maximu* 

AVERAGE 

S/N 

S/N 

0.5 

o.o 

1.5 

1*0 

2.5 

2.0 

3.5 

3.0 

4.3 

4.0 

3.5 

3.0 

6.5 

6.0 

7.5 

7.0 

6.5 

6.0 

9.3 

9.0 

10.5 

10.0 

17.5 

16.0 

20.5 

19.0 

24.5 

23.6 

26.5 

25.6 

26.5 

27.5 

29.5 

29.0 

30.5 

30.0 

31.5 

31.0 

32.5 

32.0 

33.5 

33.0 

34.3 

34.0 

35.5 

33.0 

36.5 

36.0 

37.5 

37.0 

36.5 

36.0 

39.5 

49.0 

•0.3 

40.0 

41.5 

41.0 

42.5 

42.0 

43.5 

43.0 

•4.3 

44.0 

43.5 

45.0 

•6.3 

44.0 

•7,5 

47.0 

46.5 

46.0 

•9.3 

49.0 

30.5 

30.0 

31.3 

51.0 

32,3 

32.0 

53.3 

53.0 

34.3 

34.0 

33.5 

53.0 

37.3 

36.4 

64.5 

39.1 

BER 

TOTAL 

BITS 

4. 476*01 

24761 

3.166*01 

9004 

3.066*01 

15757 

2.136*01 

9004 

2.406*01 

15757 

4.566-01 

22510 

3.136*01 

20259 

3.776-01 

9004 

3.416-01 

22510 

2.646*01 

16008 

4.626*02 

2251 

1.306*02 

9004 

6.166*03 

18000 

3.066-03 

22510 

4.746-03 

42769 

4.966-03 

49522 

3. 476-03 

36016 

3.306*05 

56275 

4.166-03 

56275 

2.406-03 

51773 

4.006-03 

76534 

3.636-03 

76534 

1.436-03 

123605 

3,006-03 

123605 

2.466-03 

140566 

2.136-03 

171076 

2,166-03 

175576 

1.646-03 

191335 

1.706-03 

109064 

2.256-03 

207092 

2.796-03 

177629 

2.116-03 

209343 

1.566-03 

106633 

1.476-03 

175570 

1.316-03 

230606 

1.766-03 

272371 

1.346-03 

256614 

1.346-03 

274622 

7.036-04 

256614 

1.266-03 

200339 

0*146-116 

169004 

1.146-03 

162972 

7.646*04 

1*1013 

7,316-04 

139562 

7.7 56-04 

105797 

TOTAL 

ERRORS 

11065 

2659 

4657 

1937 

3774 

6016 

6336 

3395 

7673 

4753 

149 

117 

147 

114 

160 

196 

125 

296 

234 

124 

306 

276 

177 

372 

366 

365 

379 

314 

322 

465 

497 

•41 

292 

256 

313 

479 

345 

379 

201 

252 

134 

163 

111 

iC? 

62 


264 
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TIME  DELAY  SPREAO*  0.25  TO  0.75  MS 
DOPPLER  SPREAO*  0.00  TO  0.25  CPS 


MAXIMUM 

AVCKA6E 

BER 

TOTAL 

TOTAL 

S/M 

S/N 

HITS 

ERRORS 

0.5 

0.0 

4.98P-Q1 

4502 

2240 

27. 5 

23.7 

3. 766-03 

29263 

110 

29.5 

26.6 

9.066*03 

15757 

143 

31.5 

30.5 

5.536*03 

24761 

137 

32.5 

32.0 

7.336-03 

16006 

132 

34.5 

33.6 

3.626-03 

45020 

163 

36.5 

35.6 

1.646-03 

67530 

111 

36.5 

37.5 

1.746-03 

78785 

137 

40.5 

3*.5 

1.426*03 

135060 

192 

41.5 

41.0 

1,456-03 

85536 

124 

43.5 

1.2.5 

6.666-04 

243106 

216 

44.5 

44.0 

9.916-04 

175578 

174 

45.5  - 

45.0 

6,166-04 

220596 

160 

46.5 

46.0 

1.206-03 

256614 

307 

47.5 

•  7.0 

9,246-04 

342152 

316 

•6.5 

46.0 

6.676-04 

321693 

279 

*9.5 

49.0 

9,296-04 

342152 

316 

50.5 

50.0 

6,436-04 

407431 

262 

51.5 

51*0 

7.576*04 

319642 

242 

52.5 

52.0 

7.216-04 

267669 

193 

53.5 

53.0 

1.206-03 

249661 

299 

54.5 

54.0 

7.146-04 

166574 

119 

55.5 

55.0 

9.046-04 

132609 

120 

57,5 

56.4 

9.696-04 

119303 

116 

61.5 

56.5 

6.376*04 

56526 

49 

TIME  OELAY  sPREADi  0,25  TO  0.75  MS 
DOPPLER  SPREAD!  0.50  TU  1.00  CPS 


MAXIMUM 

AVEHAGE 

8ER 

TOTAL 

TOTAL 

S/N 

S/M 

BITS 

ERKORS 

16.5 

15.0 

8,889-03 

11255 

100 

21.5 

20.0 

¥.059-03 

16008 

163 

25.5 

24.1 

4.119-03 

36016 

148 

26.5 

27.1 

2.6 19-03 

5*526 

153 

30.5 

29.6 

2.709-03 

65279 

176 

32.5 

31.4 

2.729-03 

765  34 

206 

33.5 

33.0 

2.229-03 

54024 

120 

3<l  .5 

34.0 

2.129-03 

49522 

105 

35.5 

35.0 

1.849-03 

67530 

124 

36.5 

36.0 

1. 559-03 

76534 

119 

36.5 

37.7 

2.729-03 

85536 

233 

39.5 

39,0 

2.079-03 

67530 

140 

MO. 5 

M0,0 

1.719-03 

60777 

104 

Ml. 5 

41.0 

2.119-03 

112550 

238 

M2. 5 

42.0 

2,099-03 

96793 

202 

M3. 5 

*3,0 

4,019-03 

83287 

334 

AM. 5 

44.0 

1.769-03 

81036 

143 

M5.5 

*5,0 

2.049-03 

90040 

184 

*7.5 

46. M 

1 , 3*9-03 

155319 

21* 

53.5 

49.3 

1.029-03 

101295 

103 

208 
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TIME  DELAY  SPREAD!  0.7b  10  1.50  MS 
DOPPLER  SPREAD!  0.25  TO  0,  .'  CPS 


MAXIMUM 

S/N 

AVERAGE 

S/M 

SER 

TOTAL 

BITS 

TOTAL 

ERRORS 

A1.5 

41.0 

1,666-02 

6753 

112 

44.5 

43.5 

5.446-03 

18006 

107 

45.5 

45.0 

5.616-03 

24761 

139 

46.5 

46.0 

/. 686-03 

24761 

194 

47.5 

47,0 

4. 966-03 

33765 

169 

48.5 

46.0 

6.966-03 

20259 

141 

49.5 

49.0 

4.566-03 

42769 

196 

53.5 

50  aO 

8.456-03 

38267 

132 

56,5 

54. V 

1.716-03 

15757 

27 

m 
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TIME  DELAY  SPREADl  0.75  TO  1.50  ms 
DOPPLEN  SPREAD1  0.50  TU  1.00  CPS 


MAXIMUM 

$/M 

AVEKAOE 

5/  M 

8ER 

TOTAL 

PITS 

TOTAL 

ERRORS 

46.5 

44.6 

4 .424*03 

38267 

169 

48.5 

47.4 

3. 629-03 

31514 

114 

51.5 

49.7 

4.599-03 

13506 

62 

216 
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TIME  DELAY  SPREAD!  0.00  HJ  1.50  MS 
DOPPLER  SPREAD!  0.00  TO  0.25  CPS 


MAXIMUM 

AVERAGE 

SER 

TOTAL 

TOTAL 

S/N 

S/N 

BITS 

ERRORS 

0.5 

0.0 

4.969-01 

4502 

2240 

14.5 

12.3 

3 .299-02 

6753 

222 

16.5 

15.7 

9.779-03 

13506 

132 

17.5 

17.0 

1.649-02 

9004 

166 

16.5 

16.0 

6.719-03 

22510 

196 

19.5 

19.0 

7.669-03 

27012 

207 

21.5 

20.6 

4.919-03 

36016 

177 

22.5 

22.0 

3.289-03 

36016 

116 

23.5 

23.0 

3.669-03 

45020 

174 

24.5 

24.0 

4.119-03 

51773 

213 

25.5 

25.0 

3. 499-03 

31514 

173 

26.5 

26.0 

6,109-03 

42769 

261 

27.5 

27.0 

2.269-03 

49522 

112 

26.5 

26.0 

3.139-03 

54024 

169 

29.5 

29.0 

6.459-03 

56275 

363 

30.5 

30.0 

3.709-03 

49522 

183 

31.5 

31.0 

4.259-03 

40518 

172 

32.5 

32.0 

4.589-03 

56275 

258 

34.5 

33.4 

2.949-03 

61036 

236 

36.5 

35.5 

1.089-03 

153066 

166 

36.5 

37.5 

1.169-03 

141613 

165 

40.5 

39.5 

9.919-04 

204641 

203 

41.5 

41.0 

1.269-03 

119303 

150 

42.5 

42.0 

8.509-04 

130558 

111 

43.5 

43.0 

9.2  -04 

184582 

171 

44,5 

44,0 

9.279-04 

207092 

192 

45.5 

45.0 

7,129-04 

283626 

<02 

46.5 

46.0 

1.039-03 

299383 

309 

47.5 

47.0 

8.639-04 

360419 

336 

46.5 

46.0 

7,669-04 

375917 

266 

49,5 

49.0 

8.409-04 

391674 

329 

SO.i 

50.0 

6.429-04 

450200 

307 

51,5 

51.0 

7.549-04 

364662 

275 

52.5 

52.0 

6.379-04 

321893 

205 

53.5 

53.0 

1.139-03 

276873 

313 

54.5 

54,0 

7.199-04 

171076 

123 

55,5 

55.0 

8.749-04 

137311 

120 

57,5 

56.4 

1,039-03 

123605 

127 

61,5 

56.5 

9.549-04 

60777 

56 

211 


TIME  DELAY  SPREAOl  0.00  TO  1.50  MS 
DOPPLER  SPREAD*  0.25  TO  0.50  CPf 


MAXIMUM 

AVERAGE 

BER 

TOTAL 

TOTAL 

S/N 

S/N 

BITS 

ERRORS 

0.5 

0.0 

4.268-01 

45020 

19280 

1.5 

1.0 

2.858-01 

16008 

5137 

2.5 

2.0 

3.098-01 

20259 

6267 

3.5 

3.0 

2.538-01 

29263 

7413 

4.5 

4.0 

2.918-01 

24761 

7202 

5.5 

5.0 

3.158-01 

38267 

12067 

6.5 

6.0 

2.798-01 

36016 

10040 

7.5 

7.0 

1.678-01 

49522 

6261 

6.5 

8.0 

2.358-01 

56275 

13227 

9.5 

9.0 

2.358-01 

29263 

6087 

10.5 

10. 0 

7.158-02 

15757 

1127 

13.5 

12.3 

1.138-02 

9004 

102 

14.5 

14.0 

7.758-02 

9004 

698 

15.5 

15.0 

4.208-02 

4502 

189 

17.5 

17.0 

2.418-02 

9004 

217 

18.5 

16.0 

1.928-02 

15757 

303 

20.5 

19.8 

1.248-02 

11255 

140 

24.5 

23.5 

4.668-03 

27012 

126 

26.5 

25.6 

3.388-03 

47271 

160 

27.5 

27.0 

3.558-03 

29263 

104 

28.5 

26.3 

3.628-03 

31514 

114 

29.5 

29.0 

3.648-03 

47271 

172 

30.5 

30.0 

4.928-03 

63028 

310 

31.5 

31.0 

3.788-03 

65279 

247 

32.5 

32.0 

2.858-03 

85279 

186 

33.5 

33.0 

4.108-03 

92291 

378 

34.5 

34.0 

3.318-03 

96793 

320 

35.5 

35.0 

1.358-03 

139562 

189 

36.5 

36.0 

2.868-03 

141613 

409 

37.5 

37.0 

2.458-03 

186833 

458 

36.5 

36.0 

2.428-03 

207092 

502 

39.5 

39.0 

2.038-03 

236355 

480 

40.5 

40.  ) 

1.758-03 

243108 

425 

41.5 

41.0 

2.128-03 

252112 

535 

42.5 

42.0 

2.048-03 

274622 

560 

43.5 

43.0 

2.328-03 

263367 

612 

44,5 

44.0 

2.108-03 

301634 

634 

45.5 

45.9 

1,608-03 

306136 

550 

46,5 

46.0 

1.808-03 

272371 

491 

47.5 

47.0 

1.618-03 

339901 

546 

46.5 

46.0 

1.918-03 

344403 

657 

49,5 

49.0 

1.508-03 

369164 

553 

50.5 

50.0 

1.168-03 

405160 

469 

51.5 

51.0 

7.668-04 

366913 

269 

52.5 

52.0 

V, 448-04 

270120 

255 

53,5 

53.0 

6,748-04 

261375 

246 

54.5 

54.0 

1.078-03 

227351 

244 

55.5 

55.0 

7.358-04 

195637 

144 

57,5 

56.4 

6,368-04 

175578 

It? 

64.5 

59.0 

7.778-04 

117052 

91 
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I 


i 


TIME  DELAY  SPREAD* 
DOPPLER  SPREAD* 


0.00  TU  1.50  MS 
0.50  TU  1.00  CPS 


5 


MAXIMUM  AVERAGE 

S/N  S/N 


BER 


TOTAL 

BITS 


TOTAL 

ERRORS 


J 


16.5 

19.5 

20.5 

21.5 

23.5 

25.5 

26.5 

27.5 

26.5 
2V.5 

30.5 

31.5 

32.5 

33.5 

34.5 

35.5 

36.5 

37.5 

36.5 

39.5 

40.5 

41.5 

42.5 

43.5 

44.5 

45.5 

46.5 

47.5 

49.5 

53.5 


15.3 

18.4 
20.0 
21.0 

22.7 

24.7 
26.0 
27.0 
26.0 
29.0 
30.0 
31*0 
32.0 
33.0 
34.0 
35.0 
36.0 
37.0 
36.0 
39.0 
40.0 
41.0 
42.0 
43.0 
44.0 
45.0 
46.0 
47.0 
46.4 
51.1 


y. 529-03 
7. 519-03 
V. 339-03 

I. 649-02 

7.709- 03 
5.169-03 
4.729-03 
2.899-03 
3.479-03 

J.  969-03 
2.299-03 
2.539-03 

2.709- 03 
2.799-03 

2.129- 03 
1.579-03 
3.049-03 
2.999-03 
2.119-03 
2.419-03 
1.619-03 
1,669-03 
1,909-03 

3.129- 03 
1.769-03 
1.729-03 
1.469-03 
2.509-03 
1,429-03 
i, 369-03 


15757 

150 

22510 

169 

13506 

126 

20259 

333 

33765 

260 

69761 

360 

54024 

255 

40518 

117 

60777 

211 

58526 

232 

69781 

160 

76534 

194 

65279 

176 

81036 

226 

96793 

205 

132809 

206 

135060 

411 

76534 

229 

108046 

228 

119303 

267 

119303 

192 

189064 

313 

144064 

274 

121554 

379 

141613 

250 

146566 

255 

128307 

167 

101295 

253 

103546 

147 

42769 

56 

213 


TIME  DELAY  bPREADl  0.00  TO  0.25  Mb 
DOPPLER  SPREAOi  0.00  TO  1.00  CPS 


MAXI  n^N 

AVERAGE 

BER 

TOTAL 

TOTAL 

b/.. 

S/N 

BITS 

ERRORS 

0.5 

0.0 

4,059-01 

20259 

6215 

1.5 

1.0 

2.539-01 

9004 

2278 

2.5 

2.0 

3.136-01 

4502 

1410 

3.5 

3.0 

2.709-01 

20259 

5476 

4.5 

4.0 

3.619-01 

9004 

3428 

5.5 

5.0 

2.579-01 

15757 

4049 

6.5 

6.0 

2.359-01 

15757 

3704 

7.5 

7.0 

1.209-01 

40516 

4666 

6.5 

6.0 

1.649-01 

33765 

5554 

9.5 

9.0 

1.909-01 

11255 

2134 

10.5 

10.0 

( ,249-02 

13506 

978 

11.5 

11.0 

3.299-02 

4502 

148 

14.5 

13.4 

4,799-02 

16006 

863 

15.5 

15.0 

3.129-02 

9004 

281 

17.5 

16.5 

1.369-02 

27012 

373 

18.5 

18.0 

1.319-02 

33765 

441 

19.5 

19.0 

9.749-03 

31514 

307 

20.5 

2  C  •  > 

7.779-0 3 

27012 

210 

21.5 

21.0 

1.039-02 

31514 

326 

22.5 

22.0 

2.619-03 

47271 

133 

23.5 

23.0 

7.079-03 

56275 

398 

24.5 

24,0 

4.649-03 

56275 

261 

25.5 

25.0 

5.369-03 

63026 

339 

26.5 

26.0 

5.699-03 

61036 

461 

27.5 

27.0 

2.429-03 

72032 

174 

26.5 

26.0 

3.429-03 

90040 

308 

29.5 

29.0 

4.639-03 

87789 

424 

30.5 

30.0 

2.519-03 

78785 

198 

31.5 

31.0 

3.209-03 

69781 

223 

32.5 

32.0 

3.059-03 

83287 

254 

33.5 

33.0 

3,149-03 

67530 

21? 

34.5 

34,0 

2.329-03 

78785 

183 

35.5 

35.0 

9.779-04 

123805 

121 

3*5 

36.0 

3.019-03 

119303 

359 

37. 5 

37.0 

2.359-03 

121554 

266 

36.5 

36.0 

1.709-03 

114601 

195 

39,5 

39.0 

1.619-03 

137311 

246 

40.5 

40.0 

1.359-03 

155319 

210 

41.5 

41.0 

1.199-03 

164323 

196 

42.5 

42.0 

1.329-03 

132609 

175 

43,5 

43.0 

1.009-03 

166625 

169 

4«.* 

44.0 

1,349-03 

162072 

217 

45.5 

45.0 

6.799-04 

209343 

;«4 

47,5 

46.5 

6.619-04 

256865 

171 

50.5 

49.1 

3.V19-04 

396176 

155 

52.5 

51.5 

3.619-04 

270120 

103 

5«.5 

53.4 

6,729-04 

175576 

116 

60.5 

55.9 

5.709-04 

103546 

59 
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-  ... 

I 

I 

TIME  DELAY  SPREADi  0.25 

TO  0.75  Mb 

I 

i 

DUPPLEH  SPREAD! 

0.00 

10  1.00  CPS 

MAXIMUM 

AVEKAGE 

SER 

TOTAL 

TOTAL  I 

S/N 

S/N 

BITS 

EKnORS  J 

i 

0.5 

0.0 

4.556-01 

29263 

13305  f 

1.5 

1.0 

3.166-01 

9004 

2859  f 

2.5 

2.0 

3.086-01 

15757 

4657  | 

3.5 

3.0 

2.156-01 

9004 

1937  | 

4.5 

4.0 

2.406-01 

15757 

3774  1 

5.5 

5.0 

3.566-01 

22510 

6016  ] 

l 

6.5 

6.0 

3.136-01 

20259 

6336 

7.5 

7.0 

3. 776-01 

9004 

3395 

6.5 

6.0 

3. 416-01 

22510 

7673  1 

9.5 

9.0 

2.646-01 

18006 

4753 

1  . 

10.5 

10.0 

6.626-02 

2251 

149 

1 

1 

16.5 

14.7 

6. 226-03 

13506 

in 

17.5 

17.0 

1.576-02 

6753 

106 

I 

16.5 

16.0 

6 .59M-0  3 

13506 

116 

! 

20.5 

19.6 

7.606-03 

20259 

156 

f 

23.5 

22.4 

4.92l»-03 

31514 

155 

f 

f 

25.5 

24.6 

4.168-03 

67530 

261 

i 

t 

26.5 

26.0 

3.536-03 

45020 

159 

% 

27.5 

27.0 

3.366-03 

47271 

159 

§: 

26.5 

26.0 

3.316-03 

56275 

186 

I 

2V. 5 

29.0 

4.626-03 

74283 

343 

| 

30.5 

30.0 

4.396-03 

103546 

455 

:4 

$ 

31.5 

31.0 

3,476-03 

112550 

390 

If 

32.5 

32.0 

3.539-03 

103546 

366 

33.5 

33.0 

3.146-03 

150817 

473 

! 

34.5 

34.0 

3.316-03 

150617 

499 

; 

35.5 

35.0 

1.426-03 

220598 

314 

> 

36.5 

36.0 

2.476-03 

236606 

589 

37.5 

37.0 

2.316-03 

209343 

464 

36.5 

36.0 

2.256-03 

2741 22 

617 

39,5 

39.0 

1.896-03 

317  191 

601 

i 

40,5 

40.0 

1.696-03 

312389 

526 

41.5 

41.0 

1.776-03 

387172 

664 

42.5 

42.0 

1.846-03 

414184 

763 

43.5 

43.0 

2.418-03 

3  9  3  V  2  5 

951 

44.5 

44.0 

1.638-03 

465957 

758 

i 

45.5 

45.0 

1.326-03 

497471 

656 

46.5 

46.0 

1.236-03 

526734 

647 

47.5 

47.0 

1.196-03 

641535 

761 

46.5 

46.0 

1.246-03 

634762 

767 

49.5 

49.0 

1 .139-03 

626029 

712 

50.5 

50.0 

9.266-04 

691057 

641 

51.5 

51.0 

7. 806-04 

*87511 

458 

i 

52.5 

52.0 

9,506-0* 

474961 

451 

* 

53.5 

33.0 

1.036-03 

443447 

457 

l 

54.5 

54.  C 

9.256-04 

326646 

t 

55.5 

55.0 

6.416-04 

27462? 

231 

56.5 

56.0 

6.416-04 

157570 

101 

37.5 

37.0 

1.176-03 

10U95 

1!9 

39.5 

36.3 

9.036-04 

135060 

12? 

64.5 

60.6 

3. 066-04 

29263 

9 

% 

i 
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1 

TIME  OEcAY  SPREaDi 
DOPPLER  SPREAD* 


0 .75  TU  1.50  HS 
O.OU  fU  1.00  CPS 


MAXIMUM 

AVEftAGE 

S/N 

'./N 

A1.5 

41.0 

AM. 5 

•  3.6 

45.5 

45.0 

•  6.5 

46.0 

•  7.5 

47.0 

•  6.5 

46.0 

49.5 

49.0 

SI. 5 

50.4 

56.5 

54.3 

TOTAL 

TOTAL 

BITS 

ERRORS 

9004 

118 

31514 

150 

31514 

167 

40518 

286 

54024 

257 

31514 

166 

49522 

225 

36267 

130 

22510 

6? 

6E« 


1.31P-02 
4.76P-03 
5.30P-03 
f .06P-03 
4.76P-U3 
5.27P-03 
4.54P-03 
3.A0P-03 
2.75P-03 


TIME  DELAY  SPREAD! 
OUPPlEK  SPREAD* 


DUAL  DIVERSITY-ALL  DATA 

U.OO  ?u  5.Uw  Ms 
0 . OU  lU  5.00  CHS 


MAXIMUM 

AVERAGE 

S/N 

6/ 

6.5 

6  .  J 

7. 5 

7.0 

6.5 

<1.0 

9.5 

v.o 

11.5 

11.0 

12.5 

12.0 

15.5 

14.6 

16.5 

16.0 

17.5 

17.0 

21.5 

20.1 

23,5 

22.5 

2s. 5 

26.5 

31.5 

26,6 

36. 5 

33.0 

37.5 

36.2 

39.5 

36.5 

*1.5 

40.5 

*5.5 

*3,6 

*7.5 

*6,5 

59,5 

69,6 

OCR 

total 

BITS 

6.156-02 

6753 

1 .650-01 

6753 

1 .450-01 

4004 

6.  /«*-u2 

6753 

1 . 390-02 

11255 

2.560-J2 

11255 

1.550-02 

22510 

1.5 Om-o? 

11255 

1 .960*02 

6753 

6.100-uj 

40510 

5.910-U3 

*5020 

4.040-03 

*2760 

1 , 75#-03 

1.12609 

4,660-04 

103546 

2.070-03 

157570 

4. 160-04 

1**07»4 

*. 430-04 

16*323 

1 .450-06 

63*543 

2.200-04 

*63706 

6.260-05 

'630MV 

TOTAL 

EKHQRS 


*15 

lit* 

1310 

45rt 

156 

266 

349 

169 

13# 

2*7 

266 

m 

233 

100 

426 

132 

155 

122 

102 

63 
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APPENDIX  G 


FSK  CROSS  AMBIGUITY  FUNCTIONS 


In  this  appendix  the  cross  ambiguity  functions  arc  documented  for  FSK 
transmissions  possessing  a  linear  frequency  variation  at  the  transition 
intervals.  We  require  the  function 


Ar0l0(\,T)  =  JSrU)S*010U  -  r)e-'2A,dt 


where  SQ(t)  and  Sj(()  are  the  two  ideal  FSK  transmissions,  and  S010O) 
is  the  non-ideal  FSK  transmission  depicted  in  Fig.  1.  The  signaling- 
element  duration  is  K  seconds,  the  mark  and  space  frequencies  are  -1/2# 
cps  and  1/2#  cps,  respectively;  and  the  total  transition  time  is  2L  seconds. 

For  |t|  >  L,  the  ambiguity  functions  assume  the  form: 

Ail0(KT)  "  J~f  exp{'  ‘Tr(^  +  X#  +  \r  +  (sgn  sin|;  •  j»  -  !r|  -iJJ 

+  -/r"(i?f\ir)  eXp  {l7T[(  Sgn  T)(  1  +  K*]  ~  *r  *  < *#n  t  ”  l  /AlJ  si  n  |tt(  I  ♦  X»)(  “  "  ~  | 

+  2^T7)  exp{"  ‘"[^2|r+^T^!Tl f  r- 

•  i  -  ,yf  ( i +  xi)  J  -  ^  i  - « k y£Z#J  | 

C°<^>  *  ***‘{,*[u«"  T)L ^  * M'  7  - 


*<|-.\I) 


*  ~  **<**  ““  NT*  (sgn  r )\£  *  (sgn  r-  1)V#  k,«  «M  M  ~  L) 


^ exp - . . -  “  *  «’  -  Xv  *  X’.'ijl 


For  1  r |  i  L,  the  ambiguity  functions  assume  the  form: 

4?10(^.T)  =  eXP  j“i77j~  ‘  2Xr  *  XH]  }  S"'{,"'  !1(  2/J} 


+  ?£^^"exp  (~i77  2M+  2\r  +  j 

ifd-t)  l  L  ^JJ 


ha 


-  i )  frt 

-  (2LVlf+L  +  T) 
LW 


-  <F 


(1- i)\ 


*/"  } 


+  ^^eXp^i77|2U2lf+M-2\T- AH-- 

lf(l  + 


-'.]} 


H(-i— /-  ( 2L\lf  +  L  -  r)  -  4>  ( 1  +  i  i  . 

2  LW  J  L  J  •> 


2 E 


(1-Xlf) 


exp  |^t77^;“  SVr- Mtj  sin  (l-.YH)^—  _  ^ 


•&S4 

+  I"  l7J  2*  -  2MT  -  2Vr  -  2U.  x  I  ) 

k  i  + 1 )  H  L  v  » > 


(1  -  i)(,YK  -  1) 


HA 


v  1  +  t 


(XM-t-L) 


4' 


vL  1\ 

—  (I  f  i  )(\»-  1)  J. 


4Kx) 


>  « 


e'* 


•  77  C 


*here  x  i«  a  complex  variable. 
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